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A Quantitative Study of the Expansion Method 
for Liquefying Helium 


By G. L. PICKARD * anp F. E. SIMON 


Clarendon Laboratory, Oxford 
* Now at University of British Columbia, Vancouver 


MS. received 12 December 1947 


ABSTRACT. The yield of the Expansion Method for producing liquid helium has been: 
measured for various starting conditions and the results are presented in the form of two 
diagrams. 

In addition some data are given on the volumes and specific heats of compressed helium 
gas; although these are only of a preliminary nature, they give all the data necessary for the 
design of such an expansion liquefier. 

‘The entropy diagram has been extended to higher pressures and the actual yield com- 
pared with the “ideal” yield of a piston and cylinder arrangement. It appears that the 
simple expansion method results in an efficiency which is 60°% of the ideal when starting 
from the usual temperature and pressure conditions. The reasons for this high efficiency 
are discussed. 


Se TNR OD Calon 


HE principle of the Expansion Method (Simon 1932a, b, 1935, 1936, 
1937, Simon and Ahlberg 1933) of liquefying gases is a very simple 
one. ‘The gas is compressed isothermally into a container C (figure 1) 
called the expansion chamber, which is pre-cooled in a low temperature bath A 
to the lowest temperature possible. The container is then 
isolated by evacuating B and the gas expanded slowly through 
valve V outside the apparatus. After this expansion a certain 
proportion of the gas remains as liquid inside the container. 
Such a method would be quite impossible with most gases, 
and is only feasible for gases with very low boiling points. ‘The 
reason is that at ordinary temperatures, say room temperature, 
the container has always a much higher heat capacity than 
the gas which it contains. With falling temperature, however, 
conditions change radically because, first, at constant pressure 
a given vessel contains more gas at low temperatures, and 
secondly, the specific heat of solid bodies disappears with the 
approach to absolute zero. 

Table 1 gives the heat capacities of a steel container of 
150 cm? capacity, designed (with the conventional safety factor) 
to stand up to 100 atmospheres, together with the heat capacities 
of the helium which it would contain at two temperatures. Figure 1- 
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Table 1. Heat capacities (cal/deg.) of steel container and helium 


Room temp. TOK 
Steel container 50 0-1 
Helium at 100 atm. : 2 20 
Ratio 25 0-005 


Thus while at room temperature the container has 25 times the heat capacity 
of the gas, at 10° it has only $%; this represents a change by a factor 5000 of the 
ratio between the capacities of the container of which the greater part is due to 
the falling off in the specific heat of the steel. It is still possible to liquefy 
hydrogen (Simon etal. 1935) by this method, but its main application is in the. 
liquefaction of helium. 

The method has been used extensively at the Clarendon Laboratory and also 
in a few other laboratories (for a full list of work see Simon 1940). In most 
of the experiments carried out, the appropriate low-temperature apparatus has 
been attached to the expansion chamber. ‘This method of working has many 
advantages, chiefly that the whole apparatus is surrounded throughout by liquid 
hydrogen and thus the heat losses are kept very low. An expansion chamber 
of about 100cm? capacity provides enough helium for an experiment of 10-20 
hours’ duration, and as a result, expansion chambers larger than about 100 cm? 
have seldom been used. ‘The expansion method is therefore generally regarded 
as a small scale method, but there is in fact no limitation to the size of an expansion 
liquefier. Just before the war the first liquefier of this type was built in which 
liquid helium could be siphoned over from an expansion chamber of 450 cm? into 
an exterior container (Cooke et al. 1939). Since the war another such liquefier 
has been described (Scott and Cook 1947) of a capacity of 400.cm?, and one of 
14 litre capacity is now under construction at the Clarendon Laboratory. With 
the widening interest in this method we think it worth while to publish some 
data, obtained in 1936, from which it is possible to decide upon the best starting 
conditions. Some earlier experiments on the liquid yields have been described 
by Simon and Ahlberg (1933) and a preliminary account of the present results 
was given by Simon in 1936. 

The main object of these experiments was to determine the effect of varying 
starting conditions on the yield of liquid helium, and while doing so, to obtain _ |} 
as much data as possible about the diagram of state and thermal properties of | 
helium, of which little was known at the time in the region of operation of the | 
expansion method. It should be emphasized, however, that these latter experi- 
ments have been regarded only as auxiliary and were not carried out with high 
accuracy. 


§2. GENERAL ARRANGEMENT 


For the experiments an expansion chamber (enclosed in a vacuum case) 
was fitted with an electrical resistance thermometer and heating coil and the 
helium was allowed to expand into a large reservoir of known volume. Measure- 
ments of the pressure and temperature of the gas in the expansion chamber 
and of the gas expanded were made at intervals during the expansion. The 
volume of helium remaining as liquid after expansion to the normal boiling 
point was measured by heating the cylinder and measuring the volume of gas 
obtained. A series of seven such experiments from different starting conditions 
gave sufhcient data for the construction of curves showing the relations between 
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the starting conditions and the yield of liquid obtained. The heating was 
carried out electrically so as to measure the heat capacity of the helium gas up 
to about 20°k. before it was let out into the reservoir for measurement. 


Se IOI WAC OF THE APPARATUS 


‘The expansion chamber was the hemispherically ended cylinder of Hecla/ AMF 
steel used in the previous work on hydrogen (Simon et al. 1935); it had an 
internal volume of 144cm® and weighted only 212gm. The heater and ther- 
mometer wires were wound directly on to the cylinder and attached by bakelite 
varnish, which was not subsequently hardened by baking. The cylinder was 
then covered with bright aluminium foil to reduce radiation effects. The 
thermometer coil of No. 47 s.w.c. Eureka wire was of about 500 resistance 
and the heater coil of No. 40 s.w.c. Eureka wire was of about 200 Q resistance. 
The wires were wound in longitudinal zig-zags up and down the side of the 
cylinder to reduce the risk of their being strained by distension of the cylinder 
at high pressures. 

The single tube to the cylinder was of cupro-nickel and served to support 
the cylinder within a brass vacuum case, which was immersed in a bath of liquid 
hydrogen. ‘The temperature of the bath could be reduced to the triple point 
(14°x.) by reducing the vapour pressure over the liquid. All the necessary 
connecting tubes for the helium supply etc. were of as small a bore as practicable 
in order to reduce the dead volume correction. ‘The pressure of the helium 
in the cylinder was measured by means of a small volume bent-tube manometer 
and the pressure in the reservoir by a precision bent-tube manometer calibrated 
against a mercury U-tube manometer. 

The resistance thermometer was calibrated from 10°-20°k. using the known 
vapour pressure-temperature relation for hydrogen after liquefying hydrogen 
in the cylinder; from 2°-5°K. it was calibrated against the vapour pressure 
of liquid helium, and in the intermediate region from 5°-10°K. by using the 
cylinder as the bulb of a helium gas thermometer. 

In performing the experiments, the vacuum case, precooled with liquid air, 
‘was immersed in a bath of liquid hydrogen and a small quantity of helium gas 
admitted to permit thermal exchange between the cylinder and the bath. Helium 
_ was compressed into the expansion chamber to a pressure of about 160 atmospheres 
and the temperature of the bath lowered to the desired starting temperature. 
When thermal equilibrium had been attained the exchange gas was pumped 
from the vacuum case to isolate the cylinder. ‘The heat inflow along the con- 
necting tube could be neglected during the expansion. 

The expansion was effected in stages, the gas being permitted to expand 
into the reservoir and the pressure and temperature of the gas in the cylinder 
| and in the reservoir measured at each stage. On reaching a pressure of | atm. 
| in the cylinder the valve to the reservoir was closed. ‘The specific heat of the 
| helium remaining in the cylinder was measured up to 20°k. and the volume 
| of this helium determined by releasing it into the reservoir. From this volume 
was calculated how much had been present as liquid at the boiling point after 
| correction for the helium in the form of vapour above the liquid. 

The “yield”? of liquid was then defined for practical purposes as the percentage 
of the volume of the cylinder which remained filled with liquid after the expansion, 
‘The range of starting conditions covered during the series of expansions enabled 


2'7--2 
2]/--2 


& Pressure (acm) 
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a yield diagram to be constructed showing the relationship between starting 


conditions and the yield obtained. Furthermore, from the measurements of | 


corresponding pressure, temperature and amount of gas contained in the known 
volume of the cylinder it was possible to draw a diagram of state for the gas under 
the conditions of the expansion method. 

With the above apparatus it was practicable to use pressures up to 160 atm. 
and temperatures as low as 14°x., which resulted in yields up to about 45%. 
It was found impossible to reach lower starting temperatures by pumping the 
bath alone, and therefore in order to extend the measurements a second apparatus. 
was used. In this apparatus the helium cylinder had a hydrogen chamber 
attached to it and by reducing the pressure over this hydrogen the temperature 
of the helium before expansion could be reduced to 11° or 10°xK. With this 


apparatus it was possible to obtain yields up to 80%. The measurements. | 
with the second apparatus were rather less accurate than those with the first | 
and no attempt was made to obtain any diagram of state or specific heat results. | 
The effect on the yield of the residual solid hydrogen in the hydrogen chamber | 
was calculated and also checked by experiment. This factor introduced a correc- 
tion of a few per cent to the yield obtained. The temperatures with this apparatus. | 
were measured by means of a simple gas thermometer of the Simon type | 


(Mendelssohn 1931). 


With both sets of apparatus the necessary corrections were made throughout | 


for the dead volume of the tube system and a correction was applied to allow for the 


volumes of the cylinders for thermal contraction. The correction for distension | 
of the cylinder due to the internal pressure was negligible except at the highest | 


pressures. 
~S4 RESULTS] 


The experimental data for the pressure-temperature relationships during | 


expansion are represented in figure 2, the curves with roman figures being those 


a 
o 


5S 
uv 
2 
= 
40 


Temperature(*K) 


| 


Starting Pressure (atm) 
Figure 2. Experimental pressure-temperature data Figure 3. Yield of liquid helium by the |} 
observed during expansions. expansion method. (Percentage of 


volume of expansion chamber re- 
maining filled.) 
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made with the main apparatus, while A and B represent the expansions made 
with the second apparatus. 

The most generally convenient method of plotting the yield results is that 
of figure 3 in which the yield is plotted against the starting pressure for a series 
of starting temperatures from 7° to 20-4°K. The extent to which the yield is 
increased by raising the starting pressure and lowering the starting temperature 
is clearly seen. A comparison between these results and those of Simon and 
Ahlberg (1933) was made. The curved form of the present yield isotherms is 
very marked, while Simon and Ahlberg whose experiments were of a less precise 
nature than the present ones and did not extend to such high yields, represented 
their yield isotherms as straight lines as a first approximation. A closer inspection 
of the original experimental results shows that a better representation is given 
by yield isotherms of similar form to the present ones. With this readjustment 
the agreement between the earlier results and the present ones is good. Some 
figures for the yields are given in table 2. 


‘Table 2. Yield of liquid helium obtained by the expansion method expressed 
as the percentage of the volume of the expansion chamber remaining 
filled after the expansion 


Starting temperature (° K.) 


8° Om 15 20°4° 
Starting 150 > 100 100 37) 10 
pressure 100 => 100 76 ap) 3 
(atm.) 50 82 43 5 0 


Some of the diagram of state data are plotted in figure 4 in the form of atomic 
volume isotherms at pressures up to 150 atmospheres. Values for the atomic 
volume are also given in table 3. 


Table 3. Atomic volume of helium gas (cm?) 


Pressure (atm.) 


5 10 20 30 50 100 150 

premp. K. 8 88-0 45°5 32-0 28-0 — — al 
10 137-0 64°5 38:5 SES 26:5 22:0 OES 

15 240-0 112-0 60:5 43°5 32:0 PRS) PANY, 

20 — 171-0 85-0 = 39-0 27-0 23:0 


The figures in table 3 agree with the 
data of Zelmanov (1944) within the 
accuracy of the present measurements 
(estimated as +3°/) except in the region 
of 20atm. pressure and 8°-12°K. where 
the present figures are about 6% greater 
than those of Zelmanov. 

The atomic heat of helium gas at 
constant volume was measured at various 
concentrations. These results show a 
considerable scatter in places and this 
is attributed to the varying conditions 
under which the readings were taken. t A — 45 
‘These measurements were made at the Pressure (atm) 
end of the experimental “ 


runs”? and Figure 4. Atomic volume of helium gas. 
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very often a shortage of liquid hydrogen caused the level in the bath to fall below 
the upper end of the vacuum case. In consequence the temperature drift a 
the cylinder increased and varied, due to increased conduction along the supporting 
tube, and to a reduction in the efficiency of the thermal insulation caused by 
desorption of gas inside the case. Although the accuracy is not higher than 
+3%, we think it worth while to give in table 4 the values of the smoothed curves. 
for two different concentrations as hardly any other measurements have yet 
been made in this region. 


Table 4. Atomic heat at constant volume (cal/deg.) 


Temperature (° K.) 


TE IO 15S 203 
Density 5 3307) 3-0 23c(0) Sill 
(gm.atoms/litre) 14 2s 2°8 2eOe 3-0; 


The only direct measurements in this region which are available from other 
sources are those of Eucken (1916) at 18° and 22°xk. for a concentration of 
9.3 gram atoms per litre, which give a value at both temperatures of 3-0 cal/deg. 

Some specific heat measurements were made below 7°K. but these are not 
quoted as they are considered unreliable on account of the phenomenon of bad 
heat exchange through the helium in the neighbourhood of the critical point. 
No effect of this nature took place during the expansion because then the heat 
absorption was homogeneous throughout the whole volume of gas. In view 
of the fact that temperature equilibrium during the expansion was established 
in less than 20 seconds at each stage, it was considered that the thermal exchange 
between the gas, cylinder and wires was good and that poor thermal conductivity 
of the material of the cylinder could not be responsible for the effects noted above 
when heating the cylinder and contents after the expansion. 


§5. DISCUSSION 

Figure 3, giving the yield isotherms, provides all the information needed 
for the construction of an expansion liquefier. The data are represented in a 
somewhat different form in figure 5, in which “‘ isoyields ” are plotted as functions. 
of starting pressure and temperature. The 
shaded areas indicate the regions of zero 
yield and yield above 100%. 

An approximate formula for the (per- 
centage) yield isotherms of figure 3 is 


Yield = 5-5 Paty,/'T +10(12 ~7) 


within the limits p=75 to 150 atm. and 
T=10° to 12°x., which covers the region 
of greatest practical interest. This formula 
represents the yield isotherms to an 
accuracy of +2°% yield. 

As can be seen from the diagrams, 
the yield increases considerably with  °5 15 20 
pressure. A practical limit to the it dar coh 
pressure is set by the dimensions of the 
expansion chamber and. the connecting 


Yield) 100 80 60 40 20 Wu 


Starting Pressure (atm) 


Figure 5. Yield diagram, curves of 
constant yield. 
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tube. We generally use pressures between 100 and 150atm., and as very satis- 
factory yields can be obtained with these pressures at starting temperatures of 
the order of 10°-12°, which are easily obtainable, there does not seem to be any 
point in going to higher pressures. The question whether a combination of 
low temperature and low pressure, or higher temperature and higher pressures, 
is preferable, and under what circumstances, has previously been fully discussed 
by Simon (1940), and we can, therefore, refer to this paper. * 

The form of the diagram of state, figure 4, presents no unusual features. It is 
of great interest, however, that at fairly low pressures volumes are shown which 
are considerably lower than that of the liquid at the boiling point. This is due 
to the fact that liquid helium is “blown up” by its considerable zero-point 
energy (Simon 1934) to a volume about four times that which would be predicted 
from the data of gas kinetics. The fact that at certain starting conditions the 
gas has a smaller volume than that of the liquid is also the explanation for the 
result that we can get yields of above 100°% (measured as the degree of filling 
of the vessel). 

The entropy diagram of helium has been published by Keesom (1942) for 
a range of temperatures and pressures. An extension of this range is possible 
by employing the values from the isentropic pressure-temperature curves 
obtained during the expansions made in the present experiments (figure 5). By 
this means isobars for pressures up to 160 atm. have been drawn on the entropy 
diagram; some figures are quoted in table 5. {t should be noted that the entropy 
is given for 1 gram atom and that the zero of entropy has been taken as that of 
the liquid at the boiling point (4-22°) and 1 atm. 


Table 5. Entropy of helium gas (cals/gm. atom °K.) 


Pressure (atm.) 


. 160 140 120 100 80 60 40 

Temp.°K. 7:5 = — == = —0-15* +0-05* 0-45 
10 0-20 0-30 0-40 0-50 0-75 1-10 1:65 

is 1-45 1-60 1°35 PNG, 25) 3°05 — 

20 MTL PSs) 3-20 3-60 4-00 — os 


* Slightly extrapolated. 


The ‘‘ Ideal’’ Yield 


These extended data cannot be regarded as very accurate but are sufficient 
to permit calculations to be made as to the yields which would be obtained by 
the operation of an expansion in an ideal piston and cylinder process. For 
suppose that 1 gram atom of helium is compressed into a cylinder by a_ piston 
and that after an ideal adiabatic expansion to | atm. pressure has taken place 
the helium remains in the form of a fraction x as liquid and (1 —.~) as gas at the 
normal boiling point. Then if the entropy of the compressed gas at the starting 


* In 1939 Van Itterbeek published some data on the liquid yields of a small number of expansions 
which were considerably higher than in the data given by us. As one of the authors has shown 
(Simon 1940), there are valid reasons why Van Itterbeek’s experiments might have been affected 
by erroneous determinations of the temperature. Later, Van Itterbeek (1943) published the results 
of some new determinations at conditions which resulted in a rather low yield ; they agree with our 
values within the limit of error. It should be added that we have now made more than 1000 
expansions at the Clarendon Laboratory, in most of which a rough check on the yield was taken, 
and that we have never found any disagreement with the results of our diagram. 
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conditions is S, and those of the liquid and gas at the boiling point are S, and S, 
respectively, we have, since the expansion is an isentropic process, 


S,=x S;+(1—x)S,: 
Now if we put S,;=0, then S,=4.65 entropy units and we have 


To express this as ‘‘yield”’ as defined previously it is necessary to compare 
the volumes occupied by a fraction x of 1 gm. atom of liquid with that of 1 gm.atom 
of gas at the starting conditions used. This latter would be the volume of the. 
expansion chamber in a practical liquefier. 

Since the atomic volume of the liquid is 32cm, the volume occupied by 
xgm.atoms is 32cm? and if the volume occupied by 1 gm.atom of gas at the 
starting conditions is V;, then the yield for an ideal expansion is 


V,q= 32x 100/V;, % =(32/V;)(1 —S;/4-65)100 %. 


From this expression it is found that the ideal yield (Y,,) may be greater than 
100°, in the case of good starting conditions for the reason already mentioned, 
namely that the atomic volume of the compressed gas is very considerably less 
than that of the liquid. In a practical apparatus of course the maximum yield 
is 100°%, as any surplus evaporates at once in the connecting tubes. 

It is interesting to compare the observed yields with these ideal yields for 
given starting conditions. In figure 6 are plotted the values of the ideal yield 
as calculated above and the observed 
yield, for starting temperatures of 10°, 
15°, 20°K. Naturally the observed yield 
approximates to the theoretical. yield at 
lower starting temperatures and higher 
pressures. It is, however, interesting 
to see how high a percentage of the 
ideal yield can actually be obtained with 
such a simple arrangement, about 60°, 
when starting from the usual pressures 
and temperatures. There are two 
reasons: firstly, the heat capacity of the 
containers 1s practically negligible in com- 
parison with that of the working gas— o 50 100 150 atm 
we have the advantage of working with aati ie 
“mathematical” walls: secondly, pres- Figure 6. Comparison of ideal and observed 

: yields from expansion method for 
sures as low as 100 atm. are very high feliumuier scott 
““reduced”’ pressures for helium, as 
its internal pressure is of the order of 
20 atm.; thus when working with helium 
we have, at 100atm. pressure, conditions corresponding to about 10000atm. 
pressure for nitrogen, which has an internal pressure about 100 times as great. 

The data already given enable us to calculate the amount of gas leaving the 
apparatus during an expansion, which has to be taken up by the gasometer. 
No diagrams have been prepared, but it should be mentioned that under average 
working conditions, say 100 to 150.atm. and 10-12°, about half the gas leaves 
the apparatus during the expansion. 


150, 


Yield (°%) 


un 
io) 


— — — Ideal yield of cylinder and piston. 
———— Observed yield of expansion method. 


A study of the expansion method for liquefying helium 413 


As has been indicated in earlier work (Simon and Ahlberg 1933), an improve- 
ment of the yield can be obtained by making use of the Joule-Thomson effect. 
‘Our data on the equation of state have been used to calculate the Joule-Thomson 
inversion curve which is found to be in quite good agreement with the direct 
experiments of Zelmanov (1940). A calculation shows that in general the 
yield from a simple expansion is so high that the gain which would result from 
using the Joule-Thomson effect as well would not be worth the complication 
of the apparatus and technique which would be necessary. 

For the larger expansion apparatus which we are now contemplating at the 
Clarendon Laboratory, giving a yield of one or more litres per expansion, the 
position may however be different, and we are studying this matter at the present 
moment. We shall also postpone a discussion on the efficiency of the method 
from the point of view of consumption of liquid hydrogen until these larger 
pieces of equipment are in operation. 
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ABSTRACT. Absolute values of the penetration depth A of a magnetic field into a 
‘superconductor at different temperatures can be obtained from the relative values indicated 
by measurements on colloids, if AA, the difference between Aata variable and at a fixed 
temperature, is known. Various methods of measuring AA are briefly reviewed, and a 
detailed account is given of one particular method in which A) is deduced from the changes 
of susceptibility of thin cylinders with temperature. The changes found for mercury 
cylinders were, within experimental error, inversely proportional to the radii as simple theory 
demands, and the temperature variation was consistent with that to be expected from the 
‘colloid results. The absolute value of ) at 0° K. is estimated as 7°6 x 10~® cm. In the case 
of tin the values of AX measured by the present method were in rough agreement with those 
.of other methods, but since no colloid results were available, the absolute value of A remains 
‘undetermined. The mercury results are discussed also in relation to critical field measure- 
ments on thin films. Some preliminary results bearing on the possible dependence of A 
on magnetic field are mentioned ; it is concluded that A does not vary appreciably in fields 


up to half the critical value. 
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Sie WINNIE DUE: TON 


HEN a superconductor is in a magnetic field, currents flow on its surface 
in such a way as to keep the magnetic field out of the interior of the 
superconductor. These currents require a certain depth in which to 

flow, and this is also the depth to which the applied magnetic field penetrates 
the superconductor. If the size of the superconductor is large compared with 
the penetration depth A, it behaves to a close approximation as if it had zero 
magnetic permeability or diamagnetic susceptibility 1/47. If, however, the 
size is comparable to or less than A the superconductor has a weaker diamagnetism. 

In the phenomenological theory of F. and H. London (1935) the penetration 
is described by the differential equation 


VAS a eee (1) 
and A is given in terms of the number, », of superconducting electrons per unit 
volume, their effective mass m, and their charge e, by 


M= mean.) 2) ee eee (2) 


Although the detailed form of penetration law predicted by (1) has not yet 
been proved experimentally, equation (2) shows the theoretical importance of 
A, in giving the ratio m/n for the superconducting electrons. 

The first direct experimental evidence abcut A came from measurements 
of the magnetic properties of small mercury particles in colloids (Shoenberg 1940) 
and showed that A was of order 10-° cm. and increased rapidly as the transition 
temperature was approached. Similar evidence came from measurements of 
the critical fields of thin mercury films (Appleyard e¢ a/. 1939), though the 
interpretation in this case involved certain additional assumptions about the 
thermodynamics of the destruction of superconductivity and was therefore 
less certain. ‘The colloid results, though indicating the order of magnitude 
of A, could not give a precise absolute estimate since the average particle size 
was not accurately known, and even if it were, an absolute estimate would involve 
assuming the truth of (1); they did, however, give a fairly accurate measure of 
the variation of A(T)/A(T)) with temperature 7, where Ty, is a fixed temperature. 
If it were possible to measure also A(T) —A(T,), which we shall denote by Ad, the 
combination of such measurements with those of A(F)/A(Ty) would at once give 
the absolute size of A. 


Fortunately, several methods are possible for measurement of AA, and we |} 


shall indicate three which have been tried successfully : 

(1) The diamagnetic susceptibility y of a long cylinder of radius 7 much | 
larger than A, should be slightly less than the value x 9(=1/47) for r infinite. . | 
As is obvious from elementary considerations, and also follows from the appro- 
priate solution of (1), we should have 


Xi Xop= l= 2A, 5 ee ee (3) 
and if r is of the order 10-$cm., 2A/r should be of order 2%. In principle, 
accurate measurement of y should give A absolutely, but this is very difficult 
in practice, and, as will be explained later, only relative changes of y as the 
temperature is varied can be measured with sufficient accuracy. If we denote 
x(To) —x(T) by Ax, we have in fact 


Ax/x = Ax/ x9 = 2AA/r, fm mabe ae 4 (4) 
and so accurate measurements of Ay/y will give AA absolutely. ‘This method 
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should work also with thin superconducting plates, but some preliminary experi- 
ments showed that the results were very sensitive to slight non-parallelism of 
the plate and the field, and the geometry of the measuring apparatus will have to 
be improved before reliable results on plates can be obtained. A full account 
of these experiments will be given in the present paper; a preliminary note has 
already been published (Desirant and Shoenberg 1947). 

(2) The mutual inductance of two coils wound on a superconducting core 
should decrease slightly as the temperature is lowered, owing to the decrease 
of A, and the change should be proportional to AA. This method was tried by 
Casimir (1940), but owing to a technical complication the experiment appeared 
to give much less temperature variation of AA than was to be expected from the 
colloid results. A repetition of the experiment, avoiding the technical com- 
plication, has been made by Laurmann and Shoenberg and has confirmed the 
validity of the colloid results, and has also given interesting evidence for anistropy 
of penetration depth in single crystals. This work will be fully described in 
the second paper of this series; a preliminary note has already been published 
(Laurmann and Shoenberg 1947). 

(3) Pippard (1947) has used an R.F. method which measures effectively the 
difference between A in the superconducting state, and the skin depth 6, of 
ordinary eddy currents when superconductivity has been destroyed by a magnetic 
field. The value of 8, is not accurately known, but for tin it is almost certainly 
independent of temperature, so that measurements of A—6,, at various tempera- 
tures give AA as a function of temperature; the results (which were in fact the 
first for tin) agree well with those of the other two methods. For mercury, 
however, 5, varies with temperature, so that the method cannot be used to give 
AA until a complete theory of the skin effect at high frequencies has been 
developed. 

Further work is in hand to extend and improve the colloid method of studying 
A, and also to investigate the validity of the London’s penetration law (equation (1)) 
and to look for a possible field dependence of A. This will be dealt with in later 


papers in this series. 
§2. EXPERIMENTAL DETAILS 
(i) Method of Measurement 


The method of measuring susceptibility was the same in principle as that used 
in the colloid experiments, but considerable improvements in sensitivity and 


_ accuracy were necessary before it was suitable for the present purpose. Anaccount 
of these improvements was given at the Physical Society’s conference in 1946 


(Shoenberg 1947), and only those features of the method important for the 
present purpose will be described in detail. The principle of the method is 
to displace the specimen smartly in a uniform magnetic field from the centre 
of one coil to the centre of a similar oppositely wound coil, connected in series 
with the first and a ballistic galvanometer. The small change of flux linkage 
produced by the movement, if the specimen is magnetized, then gives a galvano- 
meter deflection proportional to the magnetic moment of the specimen; two 
coils, rather than one, are used to minimize the disturbing effects of any fluctua- 
tions of current in the solenoid producing the uniform magnetizing field. 

The main improvements introduced for the present experiments were (1) to 
immerse the coil system in liquid helium (instead of liquid nitrogen as in the 
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colloid experiments) thus enabling a greater number of turns to be used without 
increasing the circuit resistance, (ii) to use photoelectric amplification of the 
galvanometer deflection, thus increasing the deflection about 100 times, and 
(iii) to use a null method in which the magnetic moment of the specimen 1s 
compensated by the opposite magnetic moment of a small current in a 5-turn 
coil surrounding it. The overall sensitivity of the apparatus was about 
3x 105mm. per ¢.G.s. unit of uncompensated magnetic moment, and the zero 
instability was of order 1 mm. in low magnetic fields and of order 1 cm. in a field 
of 300 gauss. 

In order to understand how the measurements were actually carried out we 
shall now describe the null method in greater detail. As shown in figure 1, 
the current in the 5-turn coil is taken from 
across a fixed resistance a in the solenoid 
circuit and is therefore a fraction a/(R +a +b) 
of the solenoid current, where R is the resis- 
tance in the decade box which controls the 
fraction, and 6 the resistance of the 5-turn 
coil and its leads. When the decade box is 
set so that the magnetic moment of the speci- 
men is compensated, we have 


Aai|(R+a+b)=,VRi, ...... (5) 


where y is the susceptibility of the specimen Figure1. Schematic circuit diagram. 
defined as I/H, V is its volume, k the field pa 80 volt battery supply. 
per unit current of the solenoid, zthe solenoid ma _ milliammeter. 
current, and A the area turns ofthe 5-turn coil. 2™ settee 1 

i o : galvanometer system. 
Thus independently of the Te field- th ___rheostats controlling current in s. 
strength, the value of y is inversely propor- 


s long solenoid ; 92:4 gauss/amp. 
tional toR+a-+b6. Our conditions were such R’ resistance box controlling compen- 
that R was a few hundred ohms when the ae current of eat ea 
specimen was superconducting, while a@ was R_ resistance box controlling current 


about 50hms and b about lohm. As the ' ae nee ay 

3 2 4 = . ti-t site halves 
liquid helium fell, the resistance of the leads "(22 000 turns = tS em 
and, therefore, b increased by about 0-3 ohm, symmetrically at centre of s. 


c : G 5-turn coil movable from top to 
and this change was meee for where its Eoiton Rate Gey. 
effect was more than 0:1 yor a fixed resistance (5:2 ohm). 

It will be noticed that equation (5) assumes 


implicitly that the whole magnetic moment of the specimen is due to the solenoid 
field. In practice there are two sources of magnetic moment independent of 7; 
first the earth’s vertical field is present in addition to the solenoid field, and 
secondly there may be small residual moments due either to ferromagnetic im- 
purities or to frozen-in moments arising from previous magnetization cycles of the 
superconductor. Ferromagnetic effects may be very nearly eliminated by a suit- 
able cycle of operations above the transition temperature of the superconductor; 
it was found that if a negative field of order 1000 gauss was switched on and off 
an appreciable negative remanence was obtained and this could then be eliminated 
by switching a positive field of about 300 gauss on and off. Provided no negative 
fields were used and positive fields did not exceed 300 gauss the remanence 
remained zero subsequently. Frozen-in moments were eliminated by always 
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warming the specimen above its transition temperature before a new series of 
measurements was started. Finally the earth’s vertical field was compensated 
as well as possible by passing through the solenoid an additional current which 
was not controlled by the main rheostats and did not pass through the resistance 
a cr the ammeter. This compensation was achieved by making the specimen 
superconducting and adjusting R’ (figure 1) until with no current in a the deflection 
was zero for pulling the specimen. 

The limitation to complete elimination of the various effects just discussed 
was the zero instability of the galvanometer. Since it was not possible to detect 
reliably deflections smaller than about 1mm., appreciable errors could still be 
caused where the measurement was in very low solenoid fields. Suppose that 
the residual uncompensated effects are such as would correspond to a magnetic 
‘moment of the specimen x (which by itself would give a deflection comparable 
to the zero instability) then equation (5) must be replaced by 

Aah ath yVRAxt yi | naa (6) 
and we see that such uncompensated effects can be allowed for by plotting 
1/(R +a +6) against 1/2 and extrapolating to zero 1/7. This was in fact always 
done, and as will be seen from the examples shown in figure 2, a linear plot is 
indeed obtained; the value of 1/(R +a+5) for zero 1/1 is proportional to the true 
value of y, and the slope of the graph is a measure of x, the uncompensated extra- 
neous effect. Confirmation of the validity of this procedure was obtained by 
deliberately increasing x (by deliberate uncompensation of the earth’s field) 
and as can be seen from figure 2 a steeper graph is obtained but still passing 
through the same 960 
point on the axis of 
1/((R+a+6). In the 
actual measurements we 
always made x as small 
as possible, and it is sat- 
isfactory that the linear 
plots were all nearly 
horizontal and about as 
often with positive as 
with negative slopes. 

The exact setting of 
R to give zero galvano- 
meter deflection proved 
rather difficult owing to 
the appearance of 


“double kicks’”’ close to 0 10 20 30 40 50 
the balance position Reciprocal of Solenoid Current (amp™') 
These ‘double kicks’”’ Figure 2. Illustrating extrapolation method of deducing true value 
; of 1/(R+a+d). 
arise from the fact that (a) S5:2:16°K., (b) S5:2:16°%., but earth’s vertical field 
the galvanometer has a giohtly uncompensated, (c) $5: 3-649°K.,  (d) HS : 4-044" x., 
finite period and so does (e) H5: 4:127° x. 

i i ism | i ing has to 

rovide completel A correction for paramagnetism of the specimen mounting 

‘ P 5 cs ov, be added to each intercept on the 1/(R+a+d) axis. This is 15-2 
Bee etive) IntegTation ‘ for (a) and (b), 9:0 for (c), 9-9 for (d) and 9:6 for (e). The turning 
| the E.M.F.s. developed 11 gown of the curves (d) and (e) corresponds to the falling off of the 
the coil. Even when the susceptibility as the field approaches the critical field. 


(I/R +a+b) x10° (ohm-!) 
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net moment of the specimen and its surrounding coil is zero, this zero moment 
is obtained effectively as the sum of slightly positive and slightly negative 
contributions (due to the detailed geometry of the arrangement), and so the 
galvanometer gets impulses first one way and then the other. 

The exact “profile” of the double kick corresponding to the balance point 
is difficult to assess, and it was found that the best accuracy could be obtained by 
interpolation, noting the galvanometer deflections for several values of R definitely 
below and above the balance value, eA 
and plotting the deflections against 
1/(R+a+6). To reduce the effect of 
zero instability several readings were 
taken for each value of R and an aver- 
age taken. A typical graph is shown 
in figure 3, and it will be seen that the 
value of 1/(R +a +6) corresponding to 
zero deflection can be assessed with an 
accuracy of a few parts in a thousand. 
The precision increases of course as 
the magnetic moment, and hence the 
uncompensated galvanometer deflec- 
tions, become larger, though a limit 


Galvanometer detlection (mm) 


360 420 


380 ; 400 
|/R+a+b x 10° (ohm) 


Figure 3. Illustrating interpolation of 1/(R+a+)) 


= reached When the zero instability to zero galvanometer deflection. H5, 4-044° k. 
begins to grow with the applied field. (a) solenoid current 20 ma., (b) 50 ma. The 
As can be seen from equation (6) interpolated values appear as two of the points 


our procedure finally gives a value of ae de 


xVk/Aa, but since k and A could be determined with an accuracy only of order 
1% and V only with an accuracy of order 5°% (see below), the method is not suitable 
to give absolute values of y with sufficient accuracy to measure the small deviation 
from yx, (see equation (3)). Moreover, equation (6) requires slight correction, 
if, as was the case in most of our measurements, the specimen is appreciably 
longer than the 5-turn coil; the sensitivity of the double coil method falls off 
if the specimen does not move exactly from the centre of one coil to the centre 
of the other, and since evidently this cannot be true for all points of a long specimen, 
the specimen and compensating coil produce different effects per unit magnetic 
moment if they occupy different lengths. This correction can be estimated by 
auxiliary experiments and amounts to nearly 10% for the longest specimens. 
Fortunately all these uncertainties are irrelevant if we are interested only in 
relative changes of y with temperature, for Ay/y does not involve the calibration 
constants at all, provided that the temperature changes do not alter any of the 
constants. ‘The only possibility of a change of calibration constants with tem- 
perature arose from possible slight movement of the Dewar vessel in which the 
coils were mounted as the helium vapour pressure was changed; such movement 
would have altered the relative positions of specimen and coils and hence the 
calibration constants, but auxiliary measurements showed that this effect was 
negligible. ‘Thus the values of Ay/x could be measured with an accuracy of 
a few parts in a thousand. 

Since the measured value of y includes the magnetic effects of the specimen 
mounting (mostly glass and distrene), it was essential to determine these separately, 
in case they too varied with temperature. This was done by measuring ina 
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field high enough to destroy superconductivity completely, and the residual 
effect was found to be paramagnetic, with an effective susceptibility (i.e. per 
unit volume of the superconductor) of the order of a few per cent of that due 
to the superconductor. It did indeed vary with temperature, roughly according 
to Curie’s law, so that this correction was quite important; it has been made in 
all the results quoted below. 


(ii) Temperature Measurement 


All temperatures quoted were deduced from the vapour pressure of the 
helium bath by means of the 1932 Leiden scale. 


(iu) The Specimens 


‘The specimens were prepared by filling fine capillaries with spectroscopically 
pure mercury or tin. In order to obtain enough total volume a large number of 
these filled capillaries were bundled together so that they were parallel to each 
other; and separated from each other by about 10 radii. In one case (H2) 
the necessary volume was obtained by winding a long length of capillary into a 
spiral with the turns well separated. In the case of the spiral (H2) and of one 
of the bundles (S5) the field was perpendicular to the length of the superconductor. 
It should be noted that with this arrangement y, is 1/27 instead of 1/47, but 
equation (3) is unaffected. ‘The gain in sensitivity is offset by the fact that fields 
only half as large can be used before superconductivity begins to be destroyed. 
Full magnetization curves of H2 and S5 were taken for the purpose of studying’ 
the intermediate state; the results of these measurements are discussed elsewhere 
| (Désirant and Shoenberg 1948). In the case of the other two specimens (H5 
and H6) the field was parallel to the cylinder lengths. Details of the specimens 
are summarized in table 1. The radii quoted were worked out from the values 
of xVk/Aa deduced from the susceptibility measurements after making due 
| allowance for the various corrections mentioned above, and assuming an approxi- 
- mate value of A in using equation (3) for y. ‘This was considered more accurate 
than direct measurement, since measurements with a travelling microscope were 
difficult, and weighing was also difficult in view of the large numbers of capillaries 
and the small masses involved. ‘The radii deduced in this way did in fact agree 
' with those estimated by travelling microscope well within the limits of error 
of the latter method. The microscope readings were useful in indicating the 
order of magnitude of the variation of radii between the individual capillaries ; 
the estimated mean deviation was 7°% of the radius for H5 and H6, and about 
20% for S5. Strictly speaking the radius deduced from V is (7?)* while that 
required in deducing AA from the experiment is 7?/r, which is slightly different 
if the mean deviation of 7 is taken into account. Assuming a Gaussian distribution 
of the radii, it is easily shown that the difference is negligible for H5 and H6, 
but amounts to 3% for S5. The corrected value of r for 55 is shown in table 2. 


§3. EXPERIMENTAL RESULTS 


For each specimen measurements were made at several temperatures from 
close to the transition temperature (T',) down to about 2:1° k., and from the changes 
of 1/(R+a+6) which is proportional to x, values of A(T) — A(2-1°K.) = AA were 
deduced. The exact rather than the approximate form of equation (4) was 
cused for this purpose, using our final estimate of A(2-1°K.) to correct x to x 
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Table 1 
4d/A(2-1) 
algo ius 7 me See 1/(R+a+b) Ax “4 
Specimen ere Description IPP TR, I( 105 eae 10 ee 
H2 22°8 Pyrex capillary 2:14 303-7-+0°5 
eel filled with 3-691 302:1+0°5 6-143 052m 
with mercury. 3-968 300°6-+40°5 11-943 1°25 
About 6:4 mg. 4-046 300:6+41-0 11:96 1-83 
HS.10 763 4-085 296:9+1°5 25:5+8 2:34 
H5 15°55 Bundle of 98 2:16 396-4+0°5 
pyrex capillaries. 3-928 392:7+0°5 7242 1-03 
Mean length 4-000 391-0-+0°5 10°5+2 1:45 
18-3 mm. 4-088 387:-4+1:0 17:5+3 2°38. 
About 19 mg. 
PASM eOQS. 172 396:9+0°5 
(2-14) (396:8+0°5) 
3-699 393-8+1:0 5943 0:52 
4-046 390-°8+0°5 11:6+2 1:82 
4-046 391-:8+1-0 9-743 1-82 
4-127 384:611:0 \ 23-8431 3-97 
—15-0f —3 
(2:14) (394-8+0°5) 
3-19 394-1-40°5 
4-044 388:°5+0°5 12-342 1:82 
4-127 379-5+0°5 29-8-+-2 3-97 
4-045 266:7+0°5 
H6 9-03 Bundle of 198 4-129 255-8--5:0. Und eee 
pyrex capillaries. ==l05) J) 
Mean length 
18-6 mm. (2-14) 277:5+0°5 
About 13 mg. 3:19 276:6+0°5 
HS.11 023 4-000 270:141:0 11-:8-43 1:45 
4-048 269:4+1:0 12:9-+3 1°87 
4-085 265:0+0°5 20:0+2 2:33 
~ (4-129) Se acane 30:1+2 4-04 
—1-0 —i1 


Notes to table 1. 


(i) The different series of values represent experiments on different occasions ; the slightly 
different values of 1/(R-++a-+-b) at 2-1° k. in different series are due to slight differences in the exact. 
geometry of the specimen relative to the measuring apparatus. The values of 1/(R-++a+D) for differ- 
ent specimens are not exactly proportional to the masses, since different coils c were used on different 
occasions. 


(ii) Of the two readings at 4:046° k. on HS, the first was taken before the low temperature. 
reading and the second after. 


(iii) The first of the readings on H5 at 4-127° x. has such large limits of error that it has been 
ignored in the calculations. 


(iv) The first series of readings on H6 was accidentally cut short, but the difference between the 
4-045° x. and 4:129° k. readings has been transferred to the next series to give an estimate of 4) at 
A 29S KK: 


Table 2 
Ax 
: oe, Aa Goan Oe 
Specimen saber tt Description Tex, VU ee ) Be gs cm. x 10° (Laurmann 
(Pippard) and 
Shoenberg) 
S5 17:8 209 pyrex 2:16 974-0+1 
capillaries 3-009 971-442 2:3+2 1:2 2°4 
in 5 layers, 3559 965-141 8-1+1 6:2 DAS) 
field applied 35985 9606)-2e 2 DO ae 2 78 9-2 
transversely. 3-649 958-612 14:0+2 9-9 12:3 
Mean length 3-667 949 +8 22:°7-+7 11-2 13-9 
7-9 mm. 
About 11 mg. 


JM.12 966 
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by means of equation (3)*. The results of the various series of measurements 

for mercury and tin are collected in tables 1 and 2; the errors indicated as + are 

rough estimates of probable error based on the individual extrapolations (similar 

to those of figure 2) which yield the appropriate values of 1/(R+a+b). Ina 

few cases it was not possible to take accurate readings at 2:1°K., and in such cases 
an estimated value of 1/(R +a +5) at 2-1°. has been putin the table for the purpose 
of calculating AA. ‘These estimates (shown in brackets in table 1) were made 
from the value at the nearest temperature used, on the basis of the colloid results 
assuming an approximate absolute value of A, or by interpolation, and the error 
involved is small, since the A~T curve is nearly flat over a wide range (e.g. below 
3-7°k. for mercury). 

The usual order of experiment was to measure at successively lower tempera- 

tures, and it was necessary to make sure that the observed changes of 1/(R +a +) 
were genuinely due to changes of penetration depth rather than to some stray 
effect which was associated with the passage of time (e.g. falling of the liquid 
| helium level in the Dewar vessel). For this purpose, in one experiment (on H5) 
‘the measurements at 4:048°K. were repeated after those at 2:1°K. had been 
|taken; as can be seen from table 1, no significant difference was found, thus 
| confirming that the observed changes were indeed only due to the temperature 
) changes. ; 
In order to derive an absolute value of A for mercury our results must now 
be compared with those from the colloid measurements (Shoenberg 1940). 
)From the latter we have derived values of AA/A(2-1°k.) at each temperature 
jused in the present experiments, and these are shown in the last column of table 1. 
i The comparison assumes} that 7, for our specimens was the same (4:17°K.) 
as for the colloid specimen, but actually this assumption was not very carefully 
i verified, since its importance was realized only later. From such few measure- 
) ments of critical field as were made, however, it is unlikely that 7, differs from 
4-17°«. by more than 0-:01°x. The simplest method of deriving (2:1°K.) is 
ito plot the observed values of AA from the present experiments against those of 
}AA/A(2:1°K.) from the colloid results. This has been done in figure 4, and it 
twill be seen that within the limits of experimental error the relation is a linear 
Jone. This confirms the validity of equation (4), since there is no systematic 
i difference between the points for specimens of different radii, and also the consis- 
itency of the present results with the colloid measurements. 

The slope of the best straight line through the points of figure 4 gives A(2:1° kK.) 
in absolute measure, and a least squares calculation using weights inversely 
proportional to the errors shown in table 1 gave A(2:1°K.)=7-9 x 10-®cm. with 
ja standard deviation of 0:3 x 10-$cm. Extrapolation to absolute zero by means 
lof the colloid results then gives A(0°K.)=7-6x10-%cm. ‘This estimate may 
twell be in error by a greater amount than indicated by the small standard deviation, 
since the calculation cannot, of course, allow for possible systematic errors in 
‘the colloid results, for possible differences in T, between the present specimens 
jand the colloid (which might introduce a systematic error) and for unsuspected 
systematic errors in the present experiments. In figure 5 the values of AA are 


* In the case of tin, A(2-1° K.) was assumed to be 5X 10~* cm. for this purpose, but since the 
tcorrection is only about $%, the uncertainty in this assumption is unimportant. 

+ It is probable that for small changes of T,, A remains the same function of T.—T, sont dias 
different for two specimens, the observed values of T must be adjusted in order to compare values 


J 
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shown plotted against AT=T, — T, (assuming T,=417°.) together with a 
graph of the colloid results assuming A(2:1°K.) =7-9 x 10° em. The broken-line 
curve comes from the results of Appleyard et al. (1939) which are discussed below. 
The results for tin can- 
not be treated in the same 
way as those for mercury, 
since no colloid measure- 
ments are as yet available. 
They may, however, be 
compared with those of 
Pippard (1947) and with 
“those of Laurmann and 
Shoenberg (1947). In mak- 
ing such a comparison due 
allowance must be made for 
the fact that S5 had a 
transition temperature of 2 | 
T,=3:738°. while the ADJA (21 °K) 
specimens of the other Figure 4. Relation between 4A from present work and 4A/A 
two experiments had. 2 deduced from colloid experiments. © H2, /\ Hs, + HO. 
T.=3-711°K., for? each | \ l | 
temperature of measure- xf er | cE 


Ad «10° cm 


| i 


4 


ment of S5 we therefore 
show in table 2 the value of 
AX obtained in each of the 
other two experiments at 
a temperature 0-028°K. 
lower. It will be seen that 
except at the highest tem- oS r 4 
perature (where owing to = ss Sas 
the limited range of fields : 
available the accuracy was 0 ai aa 1 
rather poor), the results ON Te 

agree fairly well with those Figure 5. Relation between 4A and AT (assuming To=4:17° Ke 
of Laurmann and Shoen- © AZ os a H6. Full curve deduced from colloids 


berg, but are a little higher assuming A(2-1° K.)=7:9 X 10-° cm. Broken curve from 
than those of Pippard thin film measurements of critical field explained in text. 


It is possible that this discrepancy is associated with anisotropic effects, but it 
will be convenient to defer a detailed discussion of this point to paper IT. 


§4. RELATION BETWEEN PENETRATION DEPTH 
AN DEG RTMEGAI Er ED ae 

As has already been mentioned, information about penetration depth may be 
obtained also from the measurements of critical fields of thin mercury films by | 
Appleyard et al. (1939). The interpretation is complicated by the necessity | 
of allowing for a possible difference between the surface tensions «, and a, at 
a boundary between an insulator and the normal and superconducting phases 
respectively. If the film thickness is not too small, Ginsburg (1945) has shown 
thermodynamically that the critical field h of a film of thickness 27 is given by 


hiH;=1-H(\-E Be (7) 
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where B = 8n(a, — os)/H2 
and H, is the critical field of the bulk metal. Equation (7) is valid, however, 
only if (A+8)<r, while most of the thin films of Appleyard et al. had r.SA+P. 
For such thin films the thermodynamic theory can be developed only by assuming 
a detailed form of penetration law, and Ginsburg in fact assumed the exact 
validity of the Londons’ law (1); he assumed further that the concept of a surface 
tension still applies even when r<A+f, which would seem to require some 
justification. On the basis of these assumptions Ginsburg calculated from 
the thin film data a curve of A against T’ which disagrees seriously with that of 
figure 5. The lack of agreement suggests that Ginsburg’s assumptions are 
probably unjustified. 

The analysis of their results made by Appleyard et al. suggested that within 
experimental accuracy h/H, was a function of a single temperature dependent 
parameter divided by 7. ‘They identified this parameter with the penetration 
_ depth A, but in view of the complications just discussed, it is safer to suppose 
: that it has some other meaning and we shall denote it by A’; for thick films, for 
instance, A’ should bethesameasA+f. Ifweassumed’(2:1°K.) = 10-8 x 10-8 cm.*, 
we find that the curve of AN =d’'(T) —A'(2-1°K.) against AT (the broken line of 
figure 5) agrees fairly well with the colloid curve of AA against AT. This suggests 
that A’—A is fairly independent of temperature, and if, speculatively, we identify 
W’—-X with f, this result might be taken to mean that f is about 3 x 10-® cm. at high 
and low temperatures, rising to about 4 x 10-®cm. at 3-7°K. 


$5, POSSIBLE FIELD. DEPENDENCE OF A 


An interesting possibility, first 
suggested by Ginsburg (1947), is that 
the penetration depth may depend on 
the strength of the applied magnetic . 30 
field. Most of our measurements were 
made at low fields, and up to about $H, 
there was certainly no appreciable effect 
of this kind. A few measurements were, 
however, madeat higher fields and we did 
find that y began to fall off appreciably ; 
some typical curves of $rAy/y against 
H/H, for H5 and H6 at 4-05°x. and io 
_2-16° kK. are shown in figure 6, and similar 
effects were found in longitudinal tin 
specimens. If these changes of y are 
interpreted as due to changes of X with 


2rAX/X x 10° cm 


1 
2 


field, the values of 4vAy/y would ike . He 

of A(H)—A(0); since this 

et those ( ) ( )s : : Figure 6. Relation between $v4y/yand H/He. 
quantity is greater for the thinner speci- "AC Fys4.0462K, > XOEES 2142iK- 
men H6 than for H5, this does not seem © H5 4:047° k., on another occasion, 
to be a satisfactory interpretation, though Z\ H6 4045" x. 


it should be emphasized that these measurements were not very accurate, and 
it is just possible that the differences between H5 and H6 are due to experi- 
mental error. Probably the observed effects are at least partly due to some 


* Appleyard et al. gave N(2:1° K.)=12 x 10°42 x 10~° based on the values of h!/H, for their 
thickest films, so our assumption is not inconsistent with their results. 
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slight destruction of superconductivity starting at irregularities or at the 
ends of the cylinders. A similar effect was found by Laurmann and Shoenberg 
in their repetition of the Casimir experiment, and there the apparent 
increase of X was always associated with the appearance of resistive losses. 
Although evidently this question requires further investigation, we may say 
from the present experiments that there is no appreciable increase of A up to 
1Hf,, and that for higher fields an increase cannot be excluded, though there is 
no proof of its existence. 

It may be noted that if in fact A does increase with field, this would invalidate 
the argument on which equation (7) is based, and would further complicate the 
interpretation of critical field measurements, such as those of Appleyard et al. 
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ABSTRACT. The factors are discussed which control the initial rate of increase of 
current in the arc-like discharge obtained when a condenser of several microfarads is 
discharged through a small tube containing rarefied gas. It is suggested that at higher 
pressures ionization probably proceeds “‘ thermally ’’, but that at lower pressures ionization 
takes place by direct electron collision, the cathode spot having a negligible retarding effect 
except possibly at very low pressures. The current density of emission from the cathode 
spot is at all times enormous, appearing sometimes greater than 10° amperes per cm? and. 
very seldom less than 10° amperes per cm? even for a discharge prolonged until it has the 
appearance of a normal arc. It is suggested, therefore, that the value (Druyvesteyn and 
Penning 1940) of the emission current density of about 104 amperes per cm? usually accepted 
for normal vacuum arcs is in error (Froome 1946, 1947). The rapid movement of the 
cathode spot and rate of growth of current exclude the possibility of the cathode emission. 


being derived thermionically, while the very high current densities observed favour field 
emission. 
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SINGER OD Oi Calon: 


LTHOUGH the transient arc obtained by condenser discharge through a 
tube filled with gas at reduced pressure is widely used for high speed 
photography (fcr list of references see Henry 1944), a systematic study of 

the factors determining the growth of current seems to have been neglected. 
The discharge will be termed for convenience an “‘arc”’ since it needs, in common 
with a normal arc, a cathode spot for its maintenance, although it differs consider- 
ably in other respects from a normal arc, especially during its initial stages. In 
this work a condenser of a few microfarads, charged to some 360 volts, is connected 
directly across the tube, the discharge being started by applying a brief high- 
tension pulse to an electrode wrapped round the outside of the tube, adjacent the 
cathode. The discharge tube usually has an electrode spacing of about 18 mm. 
The paper consists of two parts : firstly the study of current rise when the impe- 
dance of the leads between condenser and tube is a minimum except for a small 
inductance in series with a small resistance included in the cathode lead for 
measuring purposes; secondly the study of the cathode spot by a simple Kerr cell 
technique, where sometimes a fairly large resistance has been added to the leads 
in order that the transition to a normal arc can be studied. 


02) THE RATE OF CURRENT RISE 


Discussion will be limited to studies of times of the order of 10~’ second, since 
the effects of possible cathode vaporization and of impurities arising from this 
seem negligible in such times. It is not proposed, in this section, to trace the 
current rise above 100 amperes or so, although the current may rise to over 
1000 amperes in a microsecond if the external impedances are kept very low. 

To measure such times a cathode-ray tube is used on a repetitive system so 
that a fair sensitivity of Y-shift is obtained. ‘This should not be less than 1 cm, 
shift for 30 v. applied; otherwise too much of the condenser voltage will be taken 
up in the measuring circuit. ‘The trace is photographed with an F3-5 lens. 
With the tubes and condenser voltages used there is a variable delay of the order of 
75 to 1 microsecond, occasionally as much as 100 microseconds, between the 
application of the triggering pulse and the start of the arc. Further, if a cathode 
spot does not form during this time, a glow discharge may strike, but not an arc. 
The probability of the formation of an unwanted glow discharge can be reduced by 
making the electrode area small, and the cathode of metal of low work function, 
or of mercury; this does not upset the discharge if an arc forms, for in this case 
ionization follows a narrow intense channel. Because of this delay, an indepen- 
dent linear time base is useless for measuring durations of 10~’sec. or less unless 
the photographic writing speed of the cathode-ray tube is such that the current at 
any time in the discharge tube can be measured ona single trace. ‘This demands 
a cathode-ray tube operating at an anode voltage of at least 5000, and a faster lens, 
implying too great a sacrifice of deflection sensitivity. Thus the method adopted 
was to allow the voltage (or a fraction of it) across the discharge tube to supply the 
X-deflection, and to record the rate of growth of current, and the current, on the 
Y-shift. If a double-beam tube is used these can be recorded simultaneously, 
but if a single beam tube is used, then two exposures must be superimposed. 

This method was found to be very satisfactory, for once an arc starts, the relation 
between the potential difference across the discharge tube and the current or rate of 
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growth of current in it is always the same, mde DEES 
so that the trace obtained, even after several 
thousand superimposed flashes, is quite Trigger, 
sharp. The cathode-ray tube beam is ™' 
turned to full brilliancy for a time a little 

longer than the duration of the discharge 

plus the variable delay already mentioned. 

Figure 1 shows a typical discharge tube Tem 
and the measuring circuit employed for 
these current rise experiments. JL is a 
small inductance (about 10-7 henry) in 
series with a small resistance R (about 7000.0 
1 ohm) in the cathode lead of the discharge ¥ to 360v. vc. ~ 
tube. The junction of these, E, is earthed, Blea 
so that the potential across L or R can be Higuue = 
applied to the Y plate of the cathode-ray tube, and thus dz/dt and 7 can be | 
found if the values of L and R are known. Hence the area under a graph of dt/dz | 
against 7 gives a time scale, since 


= [ (at/di) di. 
1 


~ to Pumps etc. 


to Y- plate 


2 - 20uF 


‘The minimum detectable current is about 1 amp. and for convenience we shall 
measure time from this point. This method will resolve 10~$ second. 

Figure 2 shows the recurrent single-sweep linear time-base and pulse generator, 
the time base portion being used in the second half of this investigation. It 
provides the triggering pulse to the discharge tube, as well as switching on the 
cathode-ray tube beam to full mtensity for a time equal to that of the linear sweep. 
This can be varied between 3 and 10000 sec. in 8 steps. It is used on a VCR.97 
with anode voltage of 2500, the most sensitive plates being used for Y deflection. 
The linear sweep is used for the set up of figure 1. 

In figure 2, the thyratron V1, is used to trigger the whole sequence of events. 
It can be used as a slow relaxation oscillator, or triggered singly by a pulse applied 
through C2, or synchronized to a sub-multiple of the mains frequency by closing 
Sl. This latter alternative is used when a large number of recurrent traces is 
needed, R8, R2 being adjusted for a repetition speed of about 10 persec. Single 
triggering is used for the Kerr cell shutter investigations described later. When ! 
V1 has “fired”’, it triggers the thyratron V3, which discharges one of the time-base 
condensers C15—C22 at constant current through the beam tetrode V2, thus 
providing a single high-speed linear sweep for supplying the cathode-ray tube 
X-deflection if required. V2 should have a screen potential of about 400 v. in 
order that this may be large compared with the potential drop across V3 when in 
the conducting state. After V3 had fired, it was found that it would not recover 
until one of the tetrode screen decoupling condensers C6—C8 had been discharged. 
It was also found that for a GTIC thyratron the resistance R21 and small condensers 
C13 were needed to assist recovery. These should be as small as possible. | 
It is thus apparent that the tetrode screen decoupling condenser C6—C8 must be 
of such a value as not to discharge appreciably during the linear part of the sweep— 
but must do so before V1 produces a further triggering pulse to V3. 

During the time of linear discharge a small negative square-wave pulse appears | 
across R22; this is amplified by V4, and the resulting positive square-wave pulse 
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fed to the cathode-ray tube grid through C31, thus turning this tube to full brilliance 
for the duration of the linear discharge. It can also be used to trigger the Kerr cell 
circuit described later. The thyratron V5 can be triggered after V3 through the 
variable delay network D2 by this positive pulse from V4, or it can be triggered 
before V3, the triggering of V3 then being delayed by the variable lag D1, depend- 
‘ ing on the position of the switch $2. When V5 fires, the condenser C33 is dis- 
charged through the primary of a small induction coil, thus producing the high 
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tension pulse which finally triggers the discharge tube under investigation. ‘The 
delay networks D1, D2, are needed to ensure that the discharge tube arc starts 
whilst the cathode-ray tube beam is at full brilliance. The 8-position switches 
$3, S4, S5, D2 are ganged. 

Figures 3, 4, 5 show some experimental results for a tube of 18-5-mm. electrode 
spacing, connected across a condenser of 2pF., initially charged to 360v., at 
various pressures and natures of gas filling. 

It is seen that at about 2 mm. Hg pressure there is a marked difference in the 
rates of current rise for air, helium and argon filling. With helium the current 
rises to 5(/amp. in about 10-7 sec., whereas with air it takes about 23 times as long. 
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Even so, the speed with which the current increases is very great for both, suggest- 
ing that whatever is the emission process operating at the cathode, it can build up 
extremely rapidly. This emitting area or “‘ cathode spot” will be discussed next, 
before attempting theoretical estimates of the current rise to be expected. It has 
not been found that the nature of the cathode has a significant effect on the current 
rise, provided measurements (such as the above) are made before extensive 
vaporization takes place from it. These results of figures 3-5 have been obtained 
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with a sodium cathode with enough lead added to render it sufficiently hard for 
machining on a lathe. Before use the protective layer of oxide which forms on 
the surface is removed, the cathode quickly inserted into the discharge tube, and 


the discharge run for some thousand flashes so that the cathode spot can further 
clean the cathode under it. 


§3. THE BEHAVIOUR AND SIZE OF THE CATHODE SPOT 
A study of the cathode spot during a transient condenser discharge is important 
for two reasons. Firstly, the nature of the cathode spot is liable to influence the 
rate of growth of current. Secondly, if the arc is prolonged by addition of a 
resistance of several ohms in the lead between condenser and discharge tube, until 
it has the characteristics of a normal vacuum arc, then the small size of the tube 
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‘enables the cathode spot to be studied under much better conditions than in a 
tube designed for continuous running at comparable currents. 

A Kerr cell shutter has been used to take low power photomicrographs of the 
discharge-tube cathode at any stage during the condenser discharge. These are 
taken singly, on ciné film, and subsequently enlarged until the final overall 
magnification is about 20 diameters. The ciné film technique is of value when 
short exposures of moving objects are needed (Froome and Jarrett 1943, 1944). 

Figure 6 shows the circuit 
diagram of the Kerr cell 41+ 
shutter, which is triggered 
via a delay circuit by the 
positive pulse from the amp- 
lifier valve V4. ‘The optical 
arrangement for photograph- 
ing a tube with a mercury 
cathode is shown in figure 7. 
In this particular tube the 
distance from the top of the 
‘cathode to the optical flat eae 
covering the tube was about pak fs 
25 mm. of 

When the thyratron V6 Figure 6. 
(figure 6) is triggered, the 5-F. 
condenser is discharged into ue ee, Pa 


the primary of induction coil 2. cell _ p) Eyepiece . 

The resulting high tension de- ue gel, a oe erry [See 
veloped across the secondary cee 

then sparks over the spark gap (Objective 
S1, opening the Kerr cell =i 
shutter until S2 sparks over, 

when the shutter closes once 

more. ‘Thespark gaps S1, S2, 

must be arranged so that light 

from S1 adequately illumi- Figure 7. 
nates S2. Both Sl and S2 
can be varied by screw adjustments. By judiciously varying Sl and S2, 
--exposures can be obtained ranging from about 75 to 40 [ASEC. When it is 
_ desired to use this apparatus, the cathode-ray tube X-sweep is provided by the 
linear time-base already described, the length of exposure, time of exposure after 
initiation of discharge tube arc, and arc current being observed on the screen, 
either visually or by photography. If exposures less than about 3 psec. are 
required, only visual observation can be used, since there is a variable delay of at 
least 1 psec. between the application of the triggering pulse to the discharge tube 
and the development of the arc, even when operating under optimum conditions. 
‘There is a similar degree of uncertainty in the sparking over of the gap 51 of figure 6. 
If longer exposures are used (as when studying prolonged arcs) the cathode-ray 
tube can be photographed by repeating the whole process a sufficient number of 
times. In this case, of course, the film of the Kerr cell shutter is exposed for one 
-condenser discharge only. There is also a variation in the exposures given by the 
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Kerr cell shutter for any given setting of the gaps S1 and S2, but the exposure: 
actually obtained on the ciné film can always be observed sufficiently accurately on 
the cathode-ray tube screen. In practice, exposures have usually been made as.- 
short as possible consistent with adequate exposure of the cine film, except when 
it was desired to show movement of the cathode spot during a fairly long exposure. 
Some results for various cathodes will now be discussed ; these are all characteristic 
of a large number of similar photographs. It was not found that the pressure of 
the gas filling of the discharge tube had any significant effect on the appearance of 
the cathode spot; in all cases in this investigation the tubes have therefore been run 
at conditions of optimum ease of triggering. 


Mercury cathode 


Figure 8(a) shows.a photomicrograph of the mercury cathode during the 
whole of one discharge of 20uF. initially at 300v. The external impedances. 
between condenser and discharge tube were made a minimum, the L and R of 
figure 1 being reduced to 3 x 10-8 H. and 0-05 © respectively, and the total residual 
inductance of the condenser and leads about 10-7. Under these conditions the 
discharge lasts some 7 usec. and the current reaches a peak of 1400 amp. about 
half-way through this time. It is found that the cathode emitting area always: 
produces a similar, but not identical, photograph. The ‘‘spot”’ usually forms on 
the mercury surface adjacent to the tube wall, and expands outwards covering the 
area shown. For different discharges this area is approximately constant in size, 
but variesinshape. Figure 8 (6) shows a Kerr cell picture of an identical discharge 
with an exposure of about 10-’ second, taken with tube current at the maximum.. 
It is now apparent that the area exposed in figure 8 (a) is caused by much smaller,, 
rapidly moving emitting areas,.the current density of emission from these being, 
for this stage of the discharge, greater than 10°amp/cm?. Reducing the maximum 
current by reducing the initial condenser voltage merely seems to reduce the area: 
of emission. Before discussing 
this further, we will examine the 
behaviour of the “spot”? when 
the arc is prolonged to about 
90 usec. (thus reducing the peak 
current) by adding a resistance 
of 1-6 Q between condenser and 
discharge tube. Figure 10 shows 
the tube current and voltage drop 
plotted against time for such a as 
discharge. It also shows the Meee 
positions A, B, C of the Kerr cell exposures corresponding to the photographs in: 
figure 9 (a), (d), (c), respectively. 

Figure 9(a) shows the cathode spot for the first 5usec. of the discharge. 
Figure 9 (b) shows a longer exposure of 7 usec. taken about 18 sec. from the start, 
_and figure 9(c) shows an 8 ysec. exposure taken about 50sec. from the start, | 
when the arc has apparently become a “normal arc”. Figure 9(d) shows the 
cathode for the whole of one flash. It is easy to see how this picture is built up | 
from (a)-(c). In figure 9(a) the “spot” appears similar (though on a smaller: | 
scale) to that of the minimum impedance discharge of figure 8(a), probably | 
consisting of a very small, nebulous, rapidly moving area. In figure 9(b), the: |f 
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‘““spot’’ has meved further over the mercury surface away from the tube wall, is 
spreading out, and developing into a number of separate unit-like spots. The 
fine streaks shown indicate these moving radially outwards from the discharge 
starting point. In (c), the unit-spot formation is complete, the current being 
carried by some 15-20 spots each presumably emitting about 2amp. From this 
point the motion of the spots is random; they remain in a group, but at relatively 
large distances apart. The current density of emission is probably at least 
oe amp/cm®. In this case the average diameter of the spots appears to be about 
z mm. After this spot-formation is complete, no further change takes place; 
the number of spots merely falls with falling current. This suggests that the 
usually accepted value of the cathode spot emission current density for mercury 
arcs is at least an order low. Further photographs, with the Kerr cell shutter, of 
50amp. and 100amp. normal arcs, jjsec. after their initiation, confirm this 
(Froome 1946). There can be little doubt that an arc is ““normal’’ after this 
interval. 


Sodium cathode 


Figure 11 (a) shows the appearance of a freshly scraped sodium cathode for 
the whole of one minimum impedance discharge, under electrical conditions 
identical with those obtaining for the discharge shown in figure 8. This photo- 
graph shows an unusual, but not rare, occurrence in that the are struck simul. 
taneously in three cathode “spots”. Figure 11 (6) shows an exposure of 10~‘ sec. 
taken at the maximum of an identical discharge with the arc current at 1400 amp. 
Again it is found that the emitting area at any instant is only a fraction of the whole 
area traversed and that the emission current density is in excess of 10° amp/cm? at 
this stage. 

The result of prolonging the discharge by the addition of a series resistance of 
7 Q between condenser and tube is seen in figure 12: figure 12(a) shows the 
appearance of the cathode for the whole of one flash and figure 12(b) shows a 
2 psec. Kerr cell exposure of the cathode, taken 120 sec. from the time of arc 
initiation, with the arc current at 12amp. ‘The unit-spot formation illustrated for 
a mercury cathode in figure 9, also takes place for a sodium cathode; the number of 
emitting spots has dwindled to three (figure 12(b)). Hence, a current density 
well in excess of 2 x 10°amp/cm? is found, although the arc appears “normal”? at 
this stage, having a potential drop of some 20 v. and a negative characteristic. 


Copper cathode 


An important difference was noted between copper cathodes and mercury and 
sodium cathodes : for mercury or sodium the cathode “spots” were brilliantly 
visible during the very high current minimum impedance discharges and the early 
stages of prolonged discharges, whereas the copper cathode “spot” was not, and 
in fact was often apparently non-existent. It was found, however, that a clean, 
polished cathode was marked after one such flash in a manner not unlike the 
photographs of figures 8(a) or 11(a). This suggested that the cathode emission 
process was not significantly different, but that at the estimated velocity of the 
-emitting area for mercury or sodium (about 10° cm/sec. for a 1400-amp. “spot”’) 
the copper cathode was not vaporized and thus the “ spot”’ was much less brilliant. 
By making the conventional supposition that one-tenth the arc current was carried 
by positive ions, and assuming a cathode fall (see later) of about 20 v., calculation 
shows that a mercury cathode could be easily vaporized in the time of about 
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10-7 sec. that the cathode spot takes to move on to fresh territory, while copper — | 
would not even be heated to red heat. It has been seen, however, that for a 
prolonged discharge with mercury or sodium cathodes, the unit-spots formed after 
the arc has become “normal” move more slowly (of the order of 10% cm/sec.) and 
hence it might be expected that for a similar discharge with copper cathode these 
unit-spots should appear once the copper in their neighbourhood became 
incandescent or vaporized. ‘This was in fact found: after some 30 psec. very 
brilliant unit spots became visible. ‘They moved more slowly than for mercury 
and sodium, and were considerably brighter in appearance, presumably due to 
incandescence of the copper. ‘Their size was again very small, indicating emission 
current densities of at least 2 x 10®amp/cm?, their numbers falling with decreasing 
current. Each spot carried about 3--5 amp. 

Figure 13 shows these spots with an exposure of 4 usec. with the arc current at 
J amp., 150 psec. after the start of a discharge prolonged as for the sodium cathode 
spot photographs of figure 12. An approximate treatment for estimating the 
cathode temperature rise under the spot is as follows : 

Assume that the temperature of the surface, initially 0), is instantaneously 
raised to 6,. Then @,, at a distance x from this surface is given by an equation of 
the form 

= “y 
aot = S| exp(—y*)dy, 


where y=x/2,/(ht) and h=k/ps, k is the thermal conductivity of the cathode 
material, p its density, s its specific heat, and ¢t the time. 

Most of the heat will be concentrated in a distance x given by 
$x(kt/ps) #1. Thus for copper x=2/(0-9t/9 x 0:1)=2x10-3cm. and for 
mercury #=24/(0-02¢/14 x 0-03) =1-4<10*.em., if. =10~ sec:, the-time taken 
for the spot to move to fresh territory. We now assume that all the energy 
supplied to the cathode from the time the spot forms on it, will be used in 
heating the volume of material of depth x under it. 

The input energy to this region is about 2 x 106 w/cm? if the cathode fall is 20 v. 
and the positive ion current is one-tenth the electron current, the electron current 
density being about 10®amp/cm?. | 

Thus the temperature reached by the copper in this time is found to be 75° c. 
The temperature that would be reached by the mercury if there were no boiling 
iA ce 

Conclusions on cathode spots 


(1) At all times the emission current density is enormous, ranging between 105 
to wellin excess of 10°amp/cm?. (ii) The initial rapid movement of the emitting 
area precludes all possibility of the emission being derived thermionically. 
(11) Intense luminosity of the cathode spot is a secondary effect, only caused if the 
positive ions falling on the cathode (under the “spot’’) are able to vaporize it or 
render it incandescent. It seems probable that the mottled effects of figures 8 (a) 
and 11 (a) are due to an unevenness in the velocity of the emitting area, as though it 
were constantly slowing down and gathering speed, the cathode “‘ spot” appearing . 
brighter when moving more slowly. (iv) Previous measurements of cathode spot 
current density based upon its markings on the cathode over a long period of time 
are probably unreliable. It may be that Slepian’s (1942) “low current density 
arc” had a very high cathode current density carried by a large number of small, 
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rapidly moving areas. (v) The very high cathode current density encountered for 
both high current arcs and arcs prolonged to normality with mercury, sodium, or 
copper cathodes, favour Langmuir’s original theory (1923) that normal arcs with 
cathodes of these metals have their cathode emission derived by field emission. 
Current densities of the order of those herein described, make this theory more 
plausible than hitherto. 


$4, THEORETICAL IMPLICATIONS 


Discussion in this paper is limited to the simplest calculation of rate of growth 
of current, time is measured from the point when the current has reached | amp. 
When this has occurred there is a conducting path of positive ions and electrons 
bridging anode and cathode, but the potential difference applied across the tube is 
still substantially that across the condenser. Inside the discharge tube this 
potential difference will be formed of three parts, that due to cathode fall, that due 
to anode fall, and that due to the necessity of maintaining an electric field inside the 
plasma or ionized gas between anode and cathode. It is well known that the 
cathode fall for normal arcs is roughly equal to or less than the ionization potential 
of the gas filling and this assumption will be made for the arcs examined here, 
for increasing or decreasing the current does not alter the cathode fall of normal 
arcs appreciably. If, as seems most likely, field emission is the process by which 
electrons leave the cathode at the cathode spot, this process does not demand a 
particularly high cathode fall, and the current rise (as far as this process is con- 
cerned) will be equal to the speed at which the necessary positive ion space charge 
sheath can be built up at the cathode. ‘The anode fall arises from the need for a 
rate of ionization slightly greater just near the anode thaninthe plasma. ‘The same 
density of ionization is needed at the anode as in the plasma, but there is a slow 
drift of positive ions away from the anode, so that the anode fall is needed to make 
good such losses. Again, if we assume an anode fall equal to the gas ionization 
potential, adequate ionization can be produced. On these hypotheses we see 
that some 20-50 v. of the initial overall voltage of approximately 360 are absorbed 
by the combined anode and cathode falls. 

We shall assume that the remainder is applied to the plasma and calculate the 
rate of growth of current, on this hypothesis, from the rate of plasma ionization to 
be expected if ionization takes place by direct electron collision. In general, 
discussion is restricted to times of the order of 10’ sec., for vapour from the cathode 
and positive ions moves a negligible distance in this time. Recombination of ions 
is also neglected. If it is also assumed that the ionized gas between anode and 
cathode is a genuine“ plasma’, i.e. that it consists of equal densities of positive ions 
and electrons, the calculation is simple, for then space charge effects may be 
neglected. Potential measurements by means of probes placed along the tube 
indicate that at any instant the plasma field is approximately equal to the tube 
potential difference divided by the electrode spacing (table 1). 

If we then assume that in the plasma ionization initially takes place only by 
direct electron collision, the equations of continuity are (Thomson 1933) 
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n,=density of electrons, positive ions, 


u_, u, =drift velocity of electrons, positive ions, 
% = ionization efficiency (i.e. the number of ion pairs produced per electron 
per cm. of advance). 
x refers to distance from the cathode. 
0, 0 
Lie ie 5 te ee — 9, (2at+) =, Scare (8) 
and On_jOt=an us, 2) ee oe (9) 


or, if the discharge path is not of unit cross-section, 0N/dt=«NU, where U now © 


refers to electron drift velocity, and N is the actual number of ions present per unit 


t 
length. Hence, NE NeeD (| ade), 
0 
where N, is the initial number of ions present. But NgU =z, where z=current | 
and q=electron charge, whence 
2 
i= NogU exp (| xUdt). vee (10) 
J0 
Table 1 
Helium Argon 
‘Gas pressure (mm. Hg) WS 8 O°5 8 3 0°5 
Initial condenser voltage 310 310 310 310 310 310 
Xp (v/cm.) SOR 37-0 ee 46-1 41-5 34-6 
d (cm.) 3:0 3e0) 3-0 3:0 3-0 3-0 
Xd tity) 111 69-3 138 125 104 
Tube pedal Gigi aGe) 155 155 108 193 180 135 


Xp 1s the gradient in the plasma at (di/dt)max ; dis the electrode spacing. 


Noter.—The plasma potential gradient is obtained from the main differences in potential 
indicated by four identical “ floating ’’ probes equally spaced along the body of the 
tube. ‘The tube used had a mercury cathode. The circuit impedances were 
slightly different from those used in the measurements of figures 3, 4, 5. 


In conditions of no space charge the electric field and current in the tube at 
any time are independent of the position in the tube, i.e. X, U, and « are indepen- 
dent of x. 

If in equation (10) we measure time from the moment when z=1 amp., we get 


l=NoUgon | so a eer (11) 
t 
whence i=(U/U,) exp (J, aU dt) ; 


and U= KX where K is the electron mobility in an electric field ¥. The equation 
we have to solve is thus 


i=(X/X,)exp ([ avian). ang) | 


Unfortunately it cannot be assumed that X, and hence «, are constant with | 


time since the inductive drop of voltage in the leads and the resistive drop cause x 
to change rapidly with time. 

If the capacity of the condenser is assumed to be large, so that a negligible | 
charge flows out of it during the early stages of the discharge in which we are most 
interested, we can put 


X=d-NV,—2V;—L(dildt) aay 


| 
| 
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where d is the electrode spacing, Vy the initial applied voltage, V, the ionizing 
potential of the gas (assumed cathode and anode fall), L the total inductance and 
R the total resistance of the leads. 

a also varies with X/P where P is the residual gas pressure. For various values 
of X/P, « has been fairly accurately evaluated by several workers (Druyvesteyn and 
Penning 1940), so that equation (12) can be solved by a laborious numerical 
integration. In obtaining numerical results one is hampered by uncertainty in 
the values for the electron mobility K found by other workers. It appears to have 
been measured only at values of X/P lower than those used here, and although it is 
fairly constant over the measured range, some error may be introduced by extra- 
polation; this is, however, not likely to be great enough to alter seriously the 
conclusions reached. 

In figures 3—5 it is seen that the current rises at rates expected on the above 
hypotheses over a limited pressure range. At higher pressures (e.g. 5-5 mm. air) 
the measured rises are more rapid than expected : some additional form of ioni- 
zation must be operative, the most probable being some “ thermal’’ mechanism. 
This appears to have appreciable effect at values of X/P lower than about 10 for 
helium, 30 for argon, and 40 for air. Ionization by direct collision will account 
for the measured current growths for values of X/P probably up to at least 150 for 
these gas fillings. It must be emphasized that these results are obtained for 
pressure ranges of 2 to 40 mm. Hg for helium, 2 to 6mm. for air, and 2 to 30mm. 
for argon. <A substantial triggering pulse is needed at the higher pressure end of 
these measurements, but once the arc has started, this pulse has never been found 
to have effect upon the current-time characteristics of the resulting discharge. 
Within these limits the cathode spot appears to have no effect. 
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ABSTRACT. Bands of Ky have been photographed in absorption in the region \ 4510— 
\ 3940 a. in the first order of a 21-ft. concave grating with a dispersion of about 1:28 A- 
permm. The bands measured fall into two systems, one lying between A 4510 and A 4220 a. 
and the other between A 4160 and 43940 a. They can be represented by the following 
equations : 

d4510-A 4220 a. : v=22970-0-+ 60-60(u’) —0-20(u’?) —92-64(u’’) +0°354(u’”?), 

\ 4160-A 3940 a. : v=24627-7-+ 61:-60(u’) —0-90(u’?) + 0:0010(u’*) —0-00030(u’*) 

—92-64(u’’)+0-354(u’”), 
where u=v-+3. 
The upper states of the two systems are considered to dissociate into 4?S+5 *P and 

42S +3 2D atoms respectively. 


§1. INTRODUCTION 

OTASSIUM is known to possess a number of systems of bands in the infra-red, 
Pp visible and the ultra-violet region. ‘The infra-red bands were first photo- 
graphed by McLennan and Ainslie (1923) and their vibrational analysis has 
been carried out by Ritschl and Villars (1928) and by Crane and Christy (1930). 
The vibrational structure of the red system has been studied by Fredrickson and 
Watson (1927), and by Crane and Christy (1930) in absorption and by Loomis. 
and Nusbaum (1932) who used the technique of magnetic rotation, while the 
rotational structure of this system has been studied by Loomis (1931). A 
vibrational analysis of the blue bands has been made by Weizel and Kulp (1930) 
from the data of Walter and Barratt (1928) and by Yamamoto (1929). Several 
systems of bands have been reported in the ultra-violet region by Chakraborti 
(1936), Yoshinaga (1937) and Sinha (1945, 1945 a). In general the analyses of 
the ultra-violet bands provided by the various authors do not appear satisfactory 
in as much as the systems look fragmentary. A somewhat similar situation 
existed in case of the ultra-violet bands of Na,, but it has now been shown that 
the bands which were analysed into five different systems by Weizel and Kulp 
(1930) constitute in reality two systems only (Sinha 1947). The present investi- 
gation has therefore been undertaken with a view to improving the existing 
data and obtaining a more satisfactory analysis of the bands of K, in the blue 
and ultra-violet region. The work has so far been carried out only in the region 
A4500-A3900 a., although bands have been reported in previous works to wave- 

lengths as low as A2900 a. Work is in progress, however, in the latter region. 


§2. EXPERIMENTAL 

The bands have been studied in absorption. The absorption vessel con- 
sisted of a cylindrical steel tube, 65 cm. long and 2:5 cm. in internal diameter, 
provided with water-cooled glass or quartz windows. The central portion of 
the tube (30 cm. in length) could be electrically heated from outside to tem- 
peratures up to 1000° c. Potassium, freed from the oil in which it was stored 
by scraping off the outer layer of the oxide and washing two to three times in 
ether containing a few per cent. of alcohol, was put into a steel tube which was 
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then pushed to the centre of the absorption tube. The absorption tube was 
evacuated with a Cenco hyvac pump and filled with nitrogen at a pressure of 
a few cm. of mercury. The presence of nitrogen prevents rapid distillation of 
the metal to cooler parts when the tube is heated. A pointolite lamp was used 
to give continuous radiation. When used in proper orientation, the continuum 
given by it was free from band structure. 

Preliminary investigations to determine the conditions under which the 
bands appeared most satisfactorily in different regions were first investigated 
with a low dispersion spectrograph. It was found that with a given amount 
of nitrogen gas inside the absorption tube, a higher temperature was generally 
needed to produce the systems the smaller the wavelength. Thus when the 
total pressure of nitrogen and potassium vapour in the tube, as indicated by 
a manometer connected to it was 10-15 cm. Hg, the blue bands appeared satis- 
factorily at about 550° c. (within the system a slightly higher temperature was 
needed to obtain the long wavelength end and a slightly lower temperature to 
obtain the short wavelength end satisfactorily), and the bands between d 4160 
and A 3940 a. at about 700° c. Using these values of temperature and pressure, 
the spectrum was next photographed in the first order of a 21-ft. concave grating 
(Eagle mounting), having a dispersion of about 1:28 a. per mm. ‘The iron arc 
spectrum has been used for comparison. 


§3. APPEARANCE OF THE BANDS AND MEASUREMENTS 
The bands appear in two distinct regions: (i) 44510-A 4220 a. and (ii) A4160- 
A3940 a. ‘They are all degraded to the red. ‘he heads appear much better 


Table 1. .K, bands: A4510-A4220 a. 


Aair ” vobs— Veale Aair 


Bile vobs— Veale , 
(a.) Int. vv (cm=) (a.) Intense) ae) (nt) 
4505 :3 1 2,10 1 4343-5 10 1,0 0 
4500-0 Ps Lg 3 4338-1 5 Sipil 2 
4495-1 3 0,8 0 4332-3 i 2,0 0 
4487-8 D 2,9 3 4327-1 4 4,1 22 
4482°8 2D; 1,8 1 4320-9 7 3,0 2 
4477-7 4 0,7 0 4316-2 4 Dal 1 
4470-6 3 2,8 D 4310°0 7 4,0 0 
4465-6 3 ee =i 4304°8 4 6,1 3 
4460-3 5 0,6 ='ll 4299-0 8 5,0 1 
4453-7 2 Dei On 4294-5 5 Tell yy 
4448 -2 3 1,6 0) 4288-4 8 6,0 1 
4442-6 6 0,5 0° 4284-0 4 8,1 Ue 
4436°5 4 2,6 =i 4277-6 6 7,0 1 
4430°8 4 iW) —1 42736 6 Opi 1 
4425°5 6 0,4 —2 4269°8 3 ti 1 
4419-3 4 Dts) =1 4267-2 4. 8,0 il. 
4413-3 3 1,4 0) 42636 4 10,1 a) 
4407-7 5 0,3 —{ 4259-5 2 22 1 
4401°8 4 2,4 0 4256°8 3} 9,0 WR 
4395-9 4 1,3 0 4253-0 2) Nite 2 
4390-2 5 0,2 —1 4249-6 3 1? 1 
4378-4 6 1,2 0 4243-4 3 iil —1 
a 4372-9 7 0,1 —2 4236-7 2 1A) 1 
4367°8 5 De 1 4229°8 | 1 1552) 1 
4361-0 6 Ne 0 4223-0 1 14,1 3 
4355-1 8 0),0 —1 
4349-7 0 AI 0 | 
NN 
29 
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marked in the former region than in the latter. 
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Lines due to rotational structure 


appear over the entire region, but the branches cannot be identified at all owing 


to overlapping and closeness of structure. 


Measurements of the heads of the bands are given in tables 1 and 2. Intensities 


given are visual estimates on a sca 


le of 1-10; they have been estimated in different 


regions (separated by horizontal lines in the tables) from different spectrograms. 


Table 2. K, bands: \4165-A 3940 a. 
Aair A vobs— Veale Aair ” vobs— Veale 
(a.) Int Vv’, v (cm) (a.) Int. v'.v (cm-}) 

4164-7 1 5,10 =a 4063°5 4 a2 = 
4161-0 1 = _ 4059-9 4 5,3 3 
4159-1 1 4,9 = 4059-0 4 7,4 ae 
4155-4 2 6,10 0 4057-9 5 ZA er 
4153-2 D 3,8 =i 4054-8 4 aye 3 
4149-4 2 5,9 =o? 4033-5 6 3,0 24 
4146-0 3 7,10 — 2) 4031-1 4 54 =3 
41443 3 4,8 2h) 4029-3 2 aes — 
4138-2 3 37), 0 4024-9 6 4,0 ae 
4137°5 3 8,10 4 4022-8 4 6,1 = 
4134-6 4 5,8 1 4016-3 5 5,0 aD 
4133-1 4 2,6 0 4014-4 5 TA 0 
4131-1 4 7,9 6 4008-0 5 6,0 =< 
4127-6 5 15 1 4004-9 3 = — 
4125-6 4 6,8 2 3999-6 5 7,0 0 
4123-6 5 3,6 o-$ | 3995-1 4 14,3 ae 
4122-7 6 0,4 i 3992-6 4 120 = 
4121-9 4 = = 3991-6 5 8,0 1 
4119-8 5 5,7 0 10,1 1 
4115-1 3 = = 3989-6 3 15,3 —4 
4113-8 4 4,6 3 3986-2 3 13,2 2 
4112-8 8 1,4 0 3984-7 5 11,1 1 
4110°5 4 6,7 4 3984-1 4 9,0 2 
4108-6 6 3,5 —] 3980-8 4 a = 
4107:3 7 0,3 2 3979-8 4 14.2 Zz 
4104-7 4 5,6 3 3978-9 4 = ae 
4103-0 6 oe 0 3978-0 5 12,1 0 
4101-7 4 = = 3976-7 4 10,0 —3 
4099-1 5 4,5 0 3975-8 3 = = 
4097-4 7 1,3 1 3974:3 3 15,5 0 

6,6 = 
4095-7 4 = = 3972:1 3 13,1 1 
4093-9 5 3,4 2) 3970-3 3 11,0 0 
4092-3 8 0,2 ==) 3969-6 3 16,2 —6 
4090-2 4 5,5 0 3966-0 3 14,1 —2 
4087°5 6 23 2 3964-6 3 12,0 —6 
4084-7 4 4,4 —2 3963-2 2 172 0 
4082-7 10 i =i] 3960-6 3 15,1 0 
4081-4 5 6,5 D 3957°2 3 13,0 1 
4078-2 6 33 1 3955:3 D 16,1 =F 
4075:5 5 5,4 0 3951-6 1 14,0 2 
4073-0 4 XD) =] 3946-3 1 15,0 —5 
4069-1 2 4,3 2 3941-3 1 — _ 
4067-0 8 Te 1 | 

6,4 a 
4065-1 3 Se i | 


§4. VIBRATIONAL ANALYSIS AND MOLECULAR CONSTANTS 


All the bands measured fall into two systems, one lying between 44510 and 


A4220 a. and the other between A4160 and 43940 a. We will, for convenience, | | 
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Blue and ultra-violet bands of K, 


call them the blue and the first 


The assign- 


ultra-violet systems respectively. 


are given in tables 1 an 
and 4 respectively. 


ments of vibrational quantum numbers 
arrangements in v’, v” schemes in tables 3 
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As mentioned befcre, the blue bands were analysed by Yamamoto 
The analysis of the former ap 
present one. 


and by Weizel and Kulp (1930). 


satisfactory and is in good agreement with the 
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All the blue bands measured in the present case can be represented by the 

equation 
v = 22970-0 + 60-60u’ —0-20u'? —92-64u" +0-354u"2 ©... (1) 
where u=v+H. 

The terms in wu” are those from the work of Loomis and Nusbaum. (1932). 
‘The values of vo4.-Vveaie given in table 1 indicate fairly good agreement between 
the observed positions of the band heads and those calculated from equation (1). 

The heat of dissociation for the upper state calculated from the above formula 
by applying Birge and Sponer’s extrapolation rule is 4560 cm. This enables 
us to calculate vation, the energy of the dissociated atoms, given by the expression, 
Vatom = %o,0 + D’—D”, where the symbols used have the usual significance. 
Now v9, 9=22955 cm, D’=4560 cm and D’=4130 cm“, and therefore 
Yatom = 23385 cm}. Further, from the study of line spectra data for potassium, 
we have 42S—52S=21025 cm, 42S—32D=21540 cm and 42S—52P 
=24710 cm-?. Comparison of these with the above value of Vatom shows 
that the upper state of the blue system dissociates into 42S and 52P atoms of K. 


(b) The first ultra-violet system 

Bands in this region have previously been reported by Yoshinaga, who 
regards them as belonging to two different systems, one lying in the region 
A4160-A 4050 a. and the other in the region 44040-A 3960 a.; his systems are 
fragmentary and comparison of his measurements with table 4 indicates that all 
the bands which he measured in this region could be included in a single system. 
It appears that the presence of the 4045-6 a. line of the P series of K which 
obscures some bands between A4040 and 44050 a. led him to consider the bands 
as belonging to two different systems. ‘The region of overlap of the bands by 
the broadened P series line is indicated by asterisks in table 4. 

The bands measured in the present investigation can be represented by the 
following equation: 

v =24627-7 + 61-60u' —0-90u’? + 0:0010u'? — 0-00030u4 
— 92-640" +0:3540, eeu. (2) 

where u=v++4 as before. 

The terms in wu” have been re- 
tained, as before, from the work of 
Loomis and Nusbaum. Valuesof 50 
Vons-Veale Given. in table 2 indicate 
fairly good agreement between the ,, 
observed frequencies and those jz’ 
calculated from equation (2). (om-") 

The value of D’ cannot in this 
case be obtained from Birge and 
Sponer’s extrapolation rule. It — 20 
has therefore been determined from 
figure 1, which gives a plot of 7” c 
Ai.€. Vo, 0, V1, 0» Ya, - ++). against AT” 
(1€..V9 6 —V1,0)¥1,0 > Y2,0°* ): 2 max 
the value of 7’ corresponding to MS aT 253. 255 257 
AT’ =0, is-equal to 25590 cm, T’x 10° (cm*') 
and since v=24612cm™, D’ is Figure 1. 
approximately 980 cm7!, ATT’ curve for A4160-A3940 bands of Kg. 
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The energy of the dissociated atoms can be determined from the relation, 
Vatom= 2 max—D" and is equal to 21455 cm™}. Comparing this value of 
Vatom With the separations of the excited atomic K levels from its ground 4*S 
level, given above, it may he concluded that the dissociation products for the 
excited state of this system -f K, bands are 42S and 32D atoms of K. 

Vibrational constants for the five states of K, are summarized in table 5. 
The constants v, We, XeWey Veer Zee are those in the expression (3) and are 
given in cm-, 

V=Ve+W ll, —XWeU2+ Yew? — Swot, ~— as aaees (3) 
where u=v+H. 


Table 5. Vibrational constants for the ground and the excited states of K, 


5 Dissociation | 
States Ve We Xe we Vewe ewe D products 
139+ oe 92-64 | 0-354 3 = 4135 42gr4a2g | 
Say t 11683-6 | 69-09 | 0-153 as =e 5510 42gi42p | 
UIT y 15378-0 | 75-00 | 0-388 | 0-00437 | 0-00018 | 1780 428142P 
1Ey+(?) | 22970-0 | 60-60 | 0-20 ae a 4560 429152P  } 
94627-7 | 61:60 | 0-90 0-0010 0-00030 980 429+ 32D | 


§5. CONCLUSION 

The ground state of K, is 12,*, due to two atoms of potassium in the normal 
42,5 state. The next two states are 1%,,* and !I1, states formed from 42S and 
42P atoms. ‘Transitions between these and the ground state are responsible 
for the infra-red and the red systems of bands of K,. ‘The next higher molecular 
states revealed in absorption experiment are those involved in the blue and the 
ultra-violet systems studied in the present investigation. ‘The upper state of 
the blue system has been shown to be due to 42S +5?P atoms and from analogy 
with the states due to 42S +4?P atoms, it may be concluded that this state is of 
™,* type. The blue system should thus be due to 1X,*<'X,* transition and 
the bands should be single headed and should only have two branches. The 
general appearance of the bands suggests that they are single headed, but con- 
firmation is impossible since the dispersion used was insufficient to enable 
resolution of two separate heads. An attempt was made to employ sufficient 
dispersion so as to see clearly the branches possessed by the bands. ‘This has 
not, however, succeeded. ‘The upper state of the next system, i.e. the first 
ultra-violet system, is regarded as due to 42S +3?D atoms of potassium. It is, 
however, also equally probable that it may be due to 428+52S atoms. At 
wavelengths lower than A3900 a. several systems of bands have been proposed 
by different workers (referred to earlier), but only one system appears extensive 
(A3900-A 3700 a.). It is not, however, possible to make a correct estimate of 
D’ even for this system, and thereby to determine the states of the atoms into 
which its upper state dissociates. If however, we suppose that this state dis- 
sociates into 42S +5?P atoms, the value of D’ required for this purpose would 
be about 2000 cm=, and the state would be "I, in character. A comparison 
with the value of dissociation energy of 4II, state from 428 +42P atoms given in 
table 5 indicates that the above value of D’ is not unlikely. Investigations are 


in progress to clarify the position with regard to the bands of this system and of 
other systems at lower wavelengths. 
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Ultra-Violet Bands of Li, 


By S. P. SINHA 
Imperial College, London 


MS. received 1 August 1947 


ABSTRACT. The ultra-violet bands of Li, have been photographed in absorption 
on a quartz spectrograph, using lithium vapour in argon at a pressure of half an atmosphere, 
at temperatures between 900 and 1100°c. Bands appear between 3500 and A 3100 a. 
The band heads, which have been measured between A 3450 and A 3130 4., can be repre- 


sented by the equation : 
v=30658:5+231-5u’ —1-5u’?—351-6u""+2-6u'”. 
where u=v-+3. 
The upper state of this system of bands is considered to dissociate into 2°,S-+-3° 
atoms of lithium. 


Se INa RODWETTON 


HERE are two systems of Li, bands in the visible region: (1) A7700-A 6550 a. 

and (ii) A5600-A4500 a. These have been studied by Wurm (1929), 
Harvey and Jenkins (1930), Loomis and Nusbaum (1931) and others, 

and are due to 1X, +<+¥,* and 4II,<12,* transitions. Li, has also some bands 
in the ultra-violet region between A 3500 and A 3100 a. which have been considered 
to belong to four different systems (Vance and Huffman 1935). These systems, 
however, appear fragmentary. Moreover, the corresponding bands of Na,— 
the  3600-A3200 a. system—formerly regarded as due to several systems of 
bands have now been shown to constitute a single system (Sinha 1947). It has 
therefore been considered necessary to re-photograph the lithium ultra-violet 


bands and attempt a fresh analysis. 
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§2. EXPERIMENTAL 


_The details of the absorption tube used in the present investigation have 
been described earlier in connection with the ultra-violet bands of K, (Sinha 
1948). Lithium, freed from the oil in which it was stored, was placed at the 
centre of the tube and heated electrically from outside to temperatures up to 
1100°c. The source of continuous radiation was a 6-v. 18-amp. tungsten 
lamp which gave a fairly intense continuum down to 43100 a. A quartz spectro- 
graph (E,) was used to photograph the bands. 

In the case of sodium and potassium, experience had shown that the presence 
of a certain amount of a gas like nitrogen, which did not react chemically with 
the metal, was almost a necessity to obtain a satisfactory photograph of the ultra- 
violet bands, and nitrogen was invariably used. ‘The main function of the gas 
is probably to prevent a rapid distillation of the metal to cooler parts of the 
absorption tube. With lithium, however, it was found that as soon as the metal 
was heated, all nitrogen introduced inside the chamber was absorbed by it and 
the pressure was reduced to its initial value. Under these conditions, although 
the visible bands appeared, there was no more than a trace of the ultra-violet 
bands, nor was there any spectroscopic evidence of the formation of a lithium— 
nitrogen compound. Nitrogen had, therefore, to be discarded and argon used. 
This was not absorbed by hot lithium and in its presence the Li, bands appeared 
quite satisfactorily in the ultra-violet region. When the pressure, as read on 
a manometer connected to the absorption chamber, was about 40 cm. Hg, the 
bands appeared well developed at temperatures between 900 and 1100° c. 
(estimated from the colour of the furnace). 


§3. APPEARANCE OF FHE SPECTRUM-AND MEASUREMENTS 


The bands, which appear between A3500 and A3100 a., are all degraded 
to the red. Liuines due to rotational structure can be seen throughout the region, 
but it is not possible to draw any conclusion with regard to the number of 
branches they have, unless the bands can be photographed at a higher dispersion 
without much overlapping. It is equally difficult to say whether they are single- 
headed; for even if they were double-headed, the heads would lie too close 
together to be seen distinctly resolved at the dispersion employed. The iron 
arc has been used for comparison. Wavelengths of the band heads are given 
in table 1; intensities given are visual estimates on a scale of 1-10. 


§4. VIBRATIONAL ANALYSIS AND MOLECULAR CONSTANTS 


All the band-heads measured are found to belong to a single system. ‘The 
vibrational quantum numbers assigned to them are given in table 1, and table 


2 gives their arrangement in a Deslandres scheme. The heads can be represented 
by the equation : 


v=30658-5 + 231-5(u’) — 1-5(u'?) —351-6(u") + 2-6(u"”), 
where u=v+}. 


The terms in uw’ are those of Loomis and Nusbaum (1931). Values of 
Vops — Veale Given in table 1 indicate that the representation of the band heads 
by the above equation is satisfactory. 

The value of the energy of dissociation for the upper state, calculated from 
Birge and Sponer’s rule, is about 9000 cm-!; the value thus calculated can be 
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Table 1. Li, ultra-violet bands 


Nair (A.) Int. vv" Vobs— Veale Aair (A.) Int. v',v” — vobs— Veale 
3442-8 2 1,5 1 3266-3 2 Dis 0 
3431-2 4 0,4 =3 3253-1 10 1,0 3 
3416-3 1 235 0 3229-6 2 2,0 2 
3404-4 4 1,4 =3 3219-9 1 41 0 
3392-1 6 0,3 =) 3206-9 6 3,0 4 
3378-5 4 24 a3 3197-2 1 5,1 3 
3365-8 3 163 = 3184-2 6 4.0 1 
3353-6 10 0,2 4 3175-3 2 6,1 4 
3340-5 2 2,3 =) 3162-0 5 5,0 3 
3327°8 2 (W 1 3154-5 ) FA 3 
3315-6 9 0,1 0 3140-8 4 6,0 =4 
3290-5 4 Ae 1 3133-8 2 81 1 
Bo77'6 6 0,0 2 


. 


Table 2, Arrangement of Li, ultra-violet bands in Deslandres’ scheme 


‘Soret 1 2 3 4 pp tae 
v diff. 
ORES 05 01M Ss O12 ie 2981029472 Se 9186 
230 © ZOOM ce 23 220m 229 230 
1 ROS Go SOR 2 BOG A Bey Ss wees S OURS 
Pepys. G2) 22 Seon 0) 225 ae 272 Ciena 225 224°8 
2 30955 % 30607 29927 29590 NN 29263 
Zt Om Se) ay 219 
3 31174 
222 222 
4 | 31396 S 31048 
220 e 220 220 
5 BAG 1Gme eo 1268 
2a 213 ATROR 
6 31830 2 31481 oe 
oe 211 211 
7 31692 
208 208 
8 31900 
Mean diff. 348°5 341°5 338°5 336°7 327 


accepted only with great reserve, since the calculation involves a large range of 
extrapolation. 


$5 1SOTOPE ERFRPECT 

Lithium has two isotopes, ®Li and *Li and their relative abundance is in 
the approximate ratio 1 to 12. It would therefore appear that although it may 
not be possible to see bands due to ®Liy, strong bands of ®Li7Li may possibly be 
‘seen under favourable conditions. The separation between the heads of bands 
due to *Li’Li and *Li, may theoretically be calculated from the expression 

vi—v=(p—1)[w,'(u’) -—a,"(u")], 

where p=1/(u/p')=1-041 approximately, and the superscript z refers to the 
‘67. i’7Li molecule. Calculation from the above expression indicates that, of the 
bands of 7Li, molecule which have appeared strongly on the plate, the 0,4 0,3 0,2 
and 0,1 bands will have the components due to ®Li’Li on their long-wave side, 
i.e. on the side on which their own rotational fine structure will extend, and it 
will therefore be difficult to detect them. For the remaining strong bands, 
i.e. the 1,0 3,0 and 4,0 bands, the isotope components will lie on the short-wave 
side and may possibly be detected. Faint heads separated from the. strong 
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heads of 3,0 and 4,0 bands by 24 and 35 cm~ have been detected on the plate, | 
the corresponding calculated values being 25 and 34 cm respectively. 


$6. DISSOCIATION SP RO Dil Cars 


The energy of the dissociated atoms from the excited state of the ultra-violet 
system can be calculated from the expression: Vatom=Vo,o+D’—D", where 
the terms used have the usual significance. Since vp 9=30500 cm’, D= | 
9000 cm-! and D’’ =9200 cm! (Loomis and Nusbaum 1931), we have vatom= 
30300 cm-! Also from Li line spectrum data, we have 22S — 32S =27 205 cm™, 
22.5 —32P=30926 cm! and 22S —32D=31283 cm A comparison of these: | 
with the above value of vajom suggests that the excited state dissociates into: | 
27S +3?P atoms. | 


§7. DISCUSSION 


Altogether three systems of bands due to Li, are now known, two of which | 
are in the visible and one in the ultra-violet. 'The visible systems have been 
definitely established to be due to transitions between the ground 14,+(22S + 24S” 
atoms) and the excited 'X,*+ and 11, (22S +2?P atoms) states. Theoretical 
considerations indicate that the next higher excited states to which transitions. 
from the ground 1X,+ state are permitted will be the following: (1) 1X,* state 
from 22S+32S atoms, (ii) '2,+ and 11, states from 22S+3?P atoms and _ 
(iii) 4X,+ and 41, states from 27S +3?2D atoms. Since, from simple considera-- 
tions of Heitler-London theory, all these states would be stable, we should. 
observe five systems of bands in the region under present investigation. 
A somewhat similar discrepancy between the number of systems expected from 
considerations of Heitler-London theory and the number observed has been. 
noted also in Na, and K,. For each of these molecules three systems of bands. | 
have been found, two of which are believed to be due to 1X, *, 'I1,(m?S +(m+1)?P) 
<-1X,*(m?S + mS) and one due to transition from the ground 1X,+ state to an 
excited state dissociating into m%S-+m?2D atoms (m=3 for Nay, and for K,). 
Of these three systems, that towards the long-wave side, i.e. the A3600-A 3200 a. 
system of Na, and the A4500-A4200 a. system of K,, appears more intensely 
developed than the other two, and also appears at a lower temperature than is. 
needed for the others. From analogy with these, it is considered probable: 
that the Li, 43500-A3100 a. system corresponds to the Na, A3600-A3200 a. | 
or to the K, 4500-A4200 a. system, and its upper state dissociates into- 
2°S+3?P atoms. The systems of Li, bands corresponding to the other two: 
systems of Nay or K, bands will require higher temperatures for their development __ 
and cannot be observed with the present steel absorption tube. 

A few additional band-heads have been measured by Vance and Huffman. 
which have led them to consider that there are three more systems of bands in | 
the region investigated. Present observations, however, provide no evidence. 
of the existence of these bands. Since Vance and Huffman have not published 
the intensities of the bands which they measured, it is difficult to account for this. | 
discrepancy. 
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Ultra-Violet Bands of NaK 
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ABSTRACT. Bands due to NaK, in the region A 3820 to 44080 a., have been photo- 
_graphed in absorption in the first order of a 21-ft. concave grating and their vibrational 
analysis carried out. ‘The constants for this system are (in cm!) 
-Ve=25 201, we —95°85, xew. =0:94, D’=2320 
Oni ASO wen oes—= (40a Die 5025. 
The system is regarded as arising from transitions from the ground 12+ state to an. 
excited state dissociating into Na 3 2S+K 3 2D atoms. 


§1. INTRODUCTION 


‘HERE are three systems of bands of NaK in the infra-red and the visible 
| regions: (1) A9150-A7200 a., (11) A6150-A5600 a., and (iit) A5300- 
A4700 a. ‘They have been studied by Ritschl and Villars (1928) and 
by Loomis and Arvin (1933, 34). Some bands in the ultra-violet region, extending 
from A4020 to A3540 a., attributed to NaK have been reported by Walter and 
Barratt (1928). ‘Two different analyses of the ultra-violet bands have been 
proposed; one by Uchida (1929) and the other by. Weizel and Kulp (1930). 
Since the data given by Walter and Barratt were obtained on an instrument 
of rather low dispersion and the two analyses of these data are widely divergent, 
it has been considered desirable to re-photograph the bands with a higher dis- 
persion and attempt a fresh analysis. 


§2. EXPERIMENTAL 


The bands have been photographed in absorption. ‘lhe details of the 
apparatus used are the same as for the ultra-violet bands of K, (Sinha 1948). 
Preliminary investigations were made on a small-dispersion instrument to study 
the conditions under which the NaK bands appeared satisfactorily in different 
regions. ‘The observation of these bands is more difficult than for Na, or K, 
because they may be masked by bands due to these molecules if they have very 
strong absorption bands in any particular region. The difhculty may be over- 
come to a certain extent by adjusting the relative proportion of sodium and 
potassium in the absorption chamber. For example, in attempting to photo- 
graph NaK bands in a region likely to be overlapped by Na, bands, only a small 
proportion of sodium should be used. Thus, in the green and the yellow 
regions, both Na, and NaK are known to give bands; the former had to be 
suppressed in order to obtain NaK bands only and the trial experiments showed 
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that this was very effectively done by photographing the absorption spectrum 
of a sample of potassium in which the only sodium was that present as impurity. 


‘The bands appeared very satisfactorily at a temperature of about 500 to 600° Cal 
when the pressure inside the absorption chamber, obtained by introducing | 


nitrogen, was about 5 cm. Hg. Although the D lines appeared quite strongly 
in absorption, no trace of Na, bands was found. To photograph the ultra- 
violet bands of NaK, since overlapping by K, bands was expected, only a small 
proportion of potassium was used in the mixture. The bands appeared at 
a temperature of about 650° c. when the pressure of nitrogen inside the tube 
was about 15 cm. Hg. ‘The method would be more difficult in a region where 


both Na, and K, bands were present, but the relative intensities of the bands | 


as a function of the relative concentrations of sodium and potassium could be | 


studied as a means of deciding which belonged to Na, which to K, and which | 


to Nak. 
The bands were then photographed on high-dispersion instruments. For 
the visible bands, a tungsten lamp (6 V., 18 amp.) was used as a source of continuum 


and photographs were taken on a glass Littrow spectrograph (dispersion about | 


8 a. per mm. at A6000 and 4 a. per mm. at A4500). The ultra-violet bands 
were photographed in the first order of a 21-ft. concave grating (dispersion 
1-3 A. per mm.) using a 500 c.p. pointolite lamp to give the continuum. 


§3. APPEARANCE OF THE SPECTRUM AND MEASUREMENTS 
(a) Visible systems 

As noted above the visible systems have been studied by a number of workers, 
but since no satisfactory photograph of these bands has been published nor 
any satisfactory description given, it is considered worth while to publish the 
spectrograms which might prove useful for purposes of identification. 


The visible band-systems of NaK, taken on a glass Littrow spectrograph, 


are reproduced in the plate facing p. 452. For the yellow system (strip (a)) between 


6150 and A5600, the temperature was 500° c. and the pressure (by introducing | 


nitrogen) 5 cm. Hg; the bands are all degraded to the red, and the heads can be 
clearly detected. ‘The absorption seems stronger between 5800 and A5600 than 


at longer wavelengths. The A5300-A4700 system (strips (b) and (c)) appear at 
about 600° c.; the bands are all degraded to the red. Rotational structure. 
can be seen over the entire region, but the dispersion employed is small for its’ 


study. 


Measurements and assignments of vibrational quantum numbers are given. 
in tables 1 and 2; they agree with those due to Loomis and Arvin (1934). The 


intensities given are visual estimates on a scale of 10. 


Table 1. NaK bands between 45650 and 45970 


Aair (A.) Int. vv" Aair (A.) Int. Vv’, vu" Aair (A.) Int. vu’ 
5640°8 4 14,0 5745-4 9 7,0 5935-1 3 0,1 
5650:4 5 13,0 5763-4 7 6,0 5954-2 3 he, 
5665°3 6 12,0 5783-0 * 6 5,0 5973-0 2 253 
5676-0 3 18,2 5802°6 7 4,0 5997-0 3 eS 
5680°8 8 11,0 5823°5 6 3,0 6022-1 2 0,3 
5696-2 8 10,0 5845-2 7 2,0 6040°5 2 1,4 
Byitslc7/ 9 9,0 5869-5 * 5 1,0 6066 +2 3 0,4 
5728-2 10 8,0 5912-0* 4 iil 


* Bands not measured by Loomis and Arvin (1934). 
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Table 2. NaK bands between 4730 and A5260 


Aair (A.) Int. vv" Aair (A.) Int. U,v" Aair (A.) Int. Ue 
4702-7 4 16,0 4847-4 10 7,0 5032:4 4+ Mes: 
4717-8 4 15,0 4866-0 8 6,0 5052°6 4 13 
4733-1 + 14,0 4876°8 10 Teal 5063°5 3 24 
4748°5 4 13,0 4894-5 8 6,1 5083 °8 4 1,4 
4764-1 5 12,0 4913-5 8 551 5105-4 3 0,4 
4779-7 6 11,0 4932-2 8 4,1 5115-0 2 15 
4796-2 * 7 10,0 4943-2 5 54 5136°3 3 0,5 
4813-0 8 9,0 4962-4 i 4,2 5167°5 2 0,6 
4819-3 7 122 4982-2 7 By 5198-9 1 0,7 
4830°:3 10 8,0 5001-8 6 DD 5230°5 1 0,8 
4836°5 8 de 5011-9 5 3.3) 5261°8 1 0,9 


* Band not measured by Loomis and Arvin (1934). 
(b) Ultra-violet system 


Photographs taken on a small quartz instrument indicate that at wavelengths 
lower than A3800 a. only bands of Na, occur, and no additional bands attributable 
to NaK can be detected. There is, however, a patch of what appears at this 


, dispersion to be continuous absorption between A3820 and A4000 a.;. this 


is considered to be due to NaK. It cannot belong to Na, which is known not 
to have any bands in this region, nor can it be attributed to K,, for although K, 
has bands in this region, they would not appear without the stronger system 
of K, between A4500 and A4200, and this is absent from the spectrum. The 
spectrum was finally photographed in the first order of a 21-ft. concave grating. 


_ The strong bands lie in the region A3870-A 4010 a., while fainter ones extend to 


; A3820 a. and to 44080 a. ‘The bands are all degraded to the red. 


The wavelengths measured are given in table 3. The iron lines were used 
for comparison. Intensities given are visual estimates on a scale of 10. 

Walter and Barratt’s measurements include bands down to A3540 a. 
attributed to NaK; their measurements indicate that there are two strong 
patches of bands, one between A 3550 and 43650 a. and the other between A 3820 
and 43990 a. Of these, only the latter has been observed in the present in- 
vestigation and attributed to NaK; in the former region, no bands except those 
of Na, have been observed. It is possible that the strong patch of bands between 


23550 and A3650 noted by Walter and Barratt as NaK bands might in fact 


belong to Na,. This suggestion receives support from the fact that all their 
bands in this region can be fitted in the vibrational scheme of the first ultra- 
violet system of Na, bands and they lie over the Franck-Condon parabola. ‘Table 4 
gives. the classification of these bands. Weizel and Kulp’s arrangement does 
not include these bands at all. Uchida has included them in his arrangement 
but the intensity distribution, as given by him, does not appear satisfactory. . 
§4. VIBRATIONAL ANALYSIS AND DETERMINATION 
OF MOLECULAR CONSTANTS 

The quantum numbers assigned to the ultra-violet bands are given in tables 3 

and 5. Nearly all the bands measured can be represented by the equation 
y=25 201 +95-85(u’) —0-94(u’?2) —123-29(u") + 0-40(u), 

where u=v-+4, and the terms in wv” are those of Loomis and Arvin (1934); t 
differences between the observed values and those calculated from the ae 
equation are given in table 3. 
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Table 3. NaK bands in absorption: 4080-A 3820. 


Aair (A.) Int. Din Oe vobs— Veale Aair (A.) Int. Uv Vobs— Veale 
4079-60 1 4,9 0) : 4,0 3 
74-48 1 3'8 3 Saliba! 10: hod aC 
62:91 1 0,5 =% 06°31 4 17,8 35) 
60:61 1 4,8 == jl 03-83 4 125 =) 
55-43 1 3,7 2 01-85 4 9,3 =a3 
00:29 4 14,6 sag 
51-49 1 2,6 —2 
aA 2 15 ss 3898-39 6 1G ==(© 
41-08 2 4,7 1 96-33 5 11,4 == 
37:05 2 3,6 2 94-26 8 8,2 1 
32:24 2 2,5 —2 93-05 4 oe —3 | 
; 10,3 0 
28-12 2 1,4 4 ae Feet ste 3 | 
23.52 1 0,3 —2 ere : 
17-20 2 3,5 0 a db fe 
13-02 2 2'4 3 BS 2 oz a: 
08-49 2 13 0 85-80 8 12,4 = 
83-85 3 — = 
04-00 2 0,2 —2 6,0 1 
3997-85 2 304 1 Bate 4 ae 4 
OSE 62 1 2,3 =) 81-48 4 14,5 2 
88-63 3 iy =5 
84-33 4 0,1 == 80-13 6 16,6 3 
79-19 4 11,3 =f 
78°71 4 333 2 76:21 8 8,1 ="! 
74:12 4 Dap 0) USS 6 13,4 ==) 
69-47 4 iil 1 ; JS 10,2 = 2 
ee 5% 0,0 0 Pee vet cie 1 
4,3 0 
59-90 6 32 0 oe Sf 7,0 0 
70-05 Pe ee , 
54:86 6 Ai 0) 68-90 2 1957 3 
51-36 4 O53) as 0) 67-80 3 2538 aS 
-45-80 4 42 = 66:17 2 — — 
15,9 1 9,1 3 
43-12 2 ay : 63-96 4 144 
34 4 
= niet ES, Sy, 61-96 3 16,5 a 
59-94 3 di 0 
ESO) 2) 6,3 2 58-26 4 8,0 —3 
35-68 4 2,0 0) 56-96 4 13,3 = 5 
34-31 4 — — 54-00 3 15,4 —6 
31-90 3 52 4 
8,4 =H 3852-91 3 10,1 0 
27-00 4 4,1 2 51-79 3 1755 4 
49-93 1 — — 
24-85 5 Ops 6 48-73 2 — a 
3,0 5 : 90°.) es 
21-19 4 soe : 46-73 2 see ae 
; F62 25; 
aie 4 Ng = 43-83 2 16,4 2 
18-47 4 14,7 5 42-36 2 11,1 as 
16-21 3 16,8 =H 40-60 2 — = 
39-08 1 — a 
13-59 3 5 == 5) 36-04 2 53 0 
13-15 4 8,3 1 ; 
- 12-40 2 — — 12.2 
10-55 4 13,6 0 S404 1 Nie z 
08-80 4 — a 32-08 1 ae i! 
28°55 1 — =e 
25-81 1 18,4 = 
The heat of dissociation D’ can be calzulated from Birge and Sponer’s extra- 
polation rule. Ks value for the upper state is 2320 cm-1. The value for the | 
ground state, D”, as given by Loomis and Arvin is 5025 cm. With. 
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violet system of Na, 


reported by Walter and Barratt 
st ultra 


, arranged here as part of the fir 
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‘Table 4. Bands between A3550 and 3650 
as due to NaK 
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Voor Zo alo: Cina the energy of the dissociated atoms in the excited state, given | 
by Vatom= Yo,9 + D’ — D’, is equal to 22510cm™'. Since K4 25 —52S9=21025cm™| 
and K42S —32D=21530 cm~, and states of Na or other states of K are still 
farther apart than these, we may reasonably conclude from the above value of 
Vatom that the dissociated atoms in the excited state are Na3 2S and K5?S or 3 2Di 
It is satisfactory to note that the value of v4, calculated from the observed | 
data is higher than the value required from considerations of atomic levels. 
This can be understood, because the value of D’ from the Birge-Sponer) 
extrapolation is generally higher than the true value of heat of dissociation’ 
for such molecules (Gaydon 1946). It is, however, difficult to decide whether 
the K atom after dissociation in the upper level is in 32D or 52S state. If we | 
assume that the above values of vy and D” are correct, then in order that the K) 
atom should be in 32D or 52S state, the value of D’ should be 1345 em or) 
835 cm! respectively. Of these, in view of the above calculated value of D’,, 
2320 cm, the former seems more plausible. Hence, it may be concluded. that 
the products of dissociation are Na3 7S + K3?D atoms. 


§5. DISCUSSION 


The ground state of the NaK molecule will be due to Na3 2S and K42S atoms, | 
and will be a !X+ state. Excited states will be obtained if we consider either 
the sodium or the potassium atom to be raised to higher energy levels. ‘The 
energy levels of Na and K beginning from the lowest atomic level are in order as 
follows: Na32S or K42S, K42P, Na32P, K52S, K32D, Na42S.... Ordinarily 
we expect that in molecules like Na,, K, or NaK, the lowest molecular state 
dissociates into atoms in their lowest states, and the next higher molecular state 
into atoms in the next higher excited states and so on. Hence the molecular: 
states resulting from normal and excited atoms of Na and K will lie roughly 
in the following°order: 14+, 32+ (from Na32,S+K42S); 45+, 32+, 41, II 
(from Na37S+K4?P); 727,22, IT, $1 (from-Na3*P+K42S); 12°)°>) (irom 
Na3*S4+K52S) and 12*, 22>, 411, #11, tA, #A»(from Na3 24S = K32D een 
In transitions from the ground state, which is 12+, we are concerned only wit 
the tu* and 'I states.. Comparing the latter with similar states in Na, K 
or known state of NaK, we conclude that the 12+ will be lower than the 1II state. 
We should, therefore, have the following systems of bands beginning from the 
long wavelength side: = | 


Systems 1 and 2: 45+, 411 (Na3 2.9 +K42P)<1>+ 
Systems 3 and 4: 12+, 1Il (Na32P+K42S)<4z+ 
System 5: ape (Na3 2S + K52S)<ta+ 
Systems 6 and 7: 12*, 1IT (Na3 2S + K32D)<ta+ 


Cc 


Of these, systems 1, 2 and 4 are the infra-red, the orange and the green system 
studied by Loomis and Arvin (1934). System 3 has not been observed and ist 
believed to be overlapped by strong Na, bands. It does not seem possible to 
say whether system 3 would lie at wavelengths longer than those of system 2 or 
at shorter wavelengths. Loomis and Arvin seem inclined to think that th 
excited *X* state of system 3 is fairly deep and that the bands of this system would, 
therefore, lie in the red region where they are obscured by red bands of Nay. 

‘The present system of bands can correspond either to system 5 or to system 6. 
The energy states of K5?.S or K3?D are rather too close together to permi 


PROC. PHYS. SOC. VOL. 60, PT. 5 (s. P. SINHA) 


NaK bands: (a) yellow system, (6) and (c) green system 


ane eee 
1 ® 2 
— en 


To face paye 452 


Ultra-violet bands of Nak 453 


of an unambiguous decision with regard to the dissociation products of the upper 
state of this system. If we, however, consider the probable error involved 
in D’ we arrive at the conclusion (as shown before) that the dissociated atoms 
are Na3*S+K 3D. This will imply that the bands corresponding to system 5, 
which will lie at wavelengths longer than those for system 6, are missing. If we 
consider that like other 1+ states, the excited 42+ state involved in system 5 
is fairly deep so that the bands of this system will lie somewhere in the green 
region, it is possible that these NaK bands may not be observed at all, having 
been obscured by strong bands of the green system of Nay. 
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An Extension of Kapitza’s Theory of 5-radiation 
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Communicated by N. Feather ; MS. received 20 August 1947 


ABSTRACT. The Kapitza theory of 5-radiation is modified by the introduction of a 
parameter b specifying the radius of the initial temperature distribution. It is shown that a 
value for } inferred from the columnar ionization theory of Jaffé leads to estimates of the 
number of 5-electrons, most probable energy of the 6-electrons, and the duration of the 


~ emission process, which are in satisfactory agreement with experimental results for platinum. 


It is also shown that this satisfactory agreement holds only for quite a small range of values 


of b. 


Sie NER OD WiC ON, 
APITZA’S theory (1923) of the 6-radiation from solid bodies bombarded 
k by «-particles supposes that the energy lost by an «-particle traversing 
a solid produces a local heating of the body, that the diffusion of this 
heat through the body can be described in terms of macroscopic constants 
(e.g. the thermal conductivity and specific heat), and that those parts of the surface 


which are heated emit electrons according to Richardson’s thermionic equation. 


These thermionic electrons are supposed to constitute the 5-radiation. 
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It is assumed by Kapitza that the initial heat is localized along a mathematical | 
line coincident with the track of the «-particle. If we denote distance from this | 
line by 7, the initial temperature distribution is thus 


(T),-9 = ©, r=0; 
(T)o=0,  7>0. 


The temperature distribution at time ¢, obtained from the equation of heat 
conduction (assuming purely radial transfer), is expressed in terms of 7, the | 
specific heat c, the thermal conductivity k, the density d, and the energy lost | 
per centimetre by the «-particle Q. In principle the macroscopic constants 
c, k and d should not be used to describe phenomena localized in a space small 
compared with the spacing of the atoms in the solid; also it would be preferable — 
to assume that the initial temperature distribution occupies a finite volume. | 
Such a modification of the analysis becomes obviously necessary if, in Kapitza’s — 
treatment, we write Richardson’s equation in the form N= A’T?e-® instead 

of N=AT?*e-?'7; for otherwise, with this change, the number v of 6-electrons | 
per «-particle becomes 


y= —(A’Q?/8nked)Ei(0) = 0. 


The argument of the exponential integral function is proportional to the negative | 
reciprocal of the maximum temperature occurring at r=0; for v to be finite 
this temperature must be finite, and hence the initial distribution must occupy 
a finite volume. 


$25 THE TEMPERATURE, DISTRIBU ON 


In Jaffé’s theory (1913) of the columnar ionization produced in fluids by] 
o-particles, it is assumed that the initial ion density has a cylindrically symmetrical 
distribution proportional to exp(—r?/b?). ‘This assumption has been justified ]} 
on theoretical grounds by Lea and Kara-Michailova (1940) who show that 6} 
is independent of the energy of the ionizing particle. It is found experimentally 
that the distribution parameter bis, to a first approximation, inversely proportional 
to the density of the fluid, and independent of its nature. We shall assume 
that the initial temperature distribution produced by the passage of an «-particle 
through a solid * is also of the form exp( —r?/b?), and, further, that an estimate | 
of the magnitude of b can be obtained from the proportionality relation which 
applies for fluids. It thus appears that for platinum (see below) we may take 
b=7 x 10-$cm, or 7 4., a value which is perhaps large enough to allow the macro-| 
scopic constants to have some relevance. (It should be observed, however, 
that in the application of the columnar ionization theory to fluids it is found) 
that the numerical results obtained are very sensitive to variation of 6 and the} 
practice is to choose a value of 6 which leads to results in agreement with experi- | 
ment. In the discussion which follows it will appear that the result of the: 
calculation of the number of 8-electrons is similarly very sensitive to variations: 
of 5.) | 


* There must obviously be a certain, though extremely short, time lag between the passage of the 
a-particle and the establishment of a physically recognizable temperature distribution ; itis perhap 


in the recognition of this point that the method of the present paper differs significantly from that of 
Kapitza. | 
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For the initial temperature distribution we now write 
(Teo = ot (Ojmedb*)iexp( —77/b72), (1) 


where 7) is the temperature of the target before the entry of the z-particle. It 
follows that, at time ¢, 


O r2 
T=T) +o om XP — aes tee 
0 T red(4ht/ed + 8°) exp Aki lcd + R): 2) 
$3. THE NUMBER OF 3-ELECTRONS 
If the «-particle enters the solid normal to the surface, and we assume that 


(2) gives the temperature distribution on the surface, the total 5-emission is 
obtained by combining (2) with 


i | po Ania. dupe aren Psy (3) 


and evaluating the integral. Following Kapitza we change this to an integral 
over surface elements dw in the (r?, ¢) plane for which Tis constant. The analysis 
is considerably simplified if we drop the constant 7) of equations (2) and (3), 
and consider the range of temperature variation as being from 0 to Ty, (say) 
rather than from Ty to Ty) +Tax. If Ty) represents room temperature, say 
300°. and 7;,,, 1s of the order of 104°K. (see below), no significant error is 
introduced by this approximation. We then find that 


O2 cdb! 
ee -\ a “ serra 


and hence 
ep Ones 17 e-$IT 
w= ea —7cdb Te ar. eevee (4) 
ey bent —0,7—0, we have T= 7... = 7, = O/mcdb", so that 
A’ OB? (7m ( T. T e 
= eke —$/T a ne —$/T oe —— @ e ee ee 
ye 8 {2 e ae har (5) 


Tf 6= —4/T, 0, = —$/Tm (5) becomes 


== A'dQb* °m 1 Or é’ 
enn ["(G-B)oe 


Pes A'¢Q6? 1 On : 1 \e’™ 
giving v= ae — 3) Bil) a (1 3 ied! Tapa en Boe" (6) 


For the «-rays from polonium falling at the beginning of their range on a 
platinum target, O~2x10-!°cal/cm. For Pt, du Bridge gives A’=1-02 x 108 
electrons/sec. cm? deg”, and ¢=7-28x10!°K. ‘Taking c=0-05cal/gm., corre- 
sponding to a temperature in the neighbourhood of 2x10?°K., we find 
T,,=1-:21x104°K. From (6), therefore, we obtain v=11 electrons/«-particle. 

Experimental values for v lie between about 0:5 and 15; a value between 


5 and 10 seems most probable. 


§4. THE ENERGY OF THE 65-ELECTRONS 
As indicated by Kapitza, (4) may be regarded as specifying the “ distribution 
in temperature” of the emitted electrons. ‘The equation may be rewritten as 


/ 2 i £ 
aA ELSA EEN aah evar a a} 


ee eae Te 


4 


456 R. M. Sillitto 


which clearly vanishes when 7'=0 or 7,,. dv must therefore have a maximuny} 
at some intermediate value of 7, which can be found by differentiating the 
expression with respect to 7, and equating to zero. There results the equation} 
T? +467? 4128 2617 =0). ee (7) | 

whose solution may be called the “ most probable temperature”’ of the emitted] 
electrons, and written T,. The ratio T,/Z, is found to depend on 6, and tc} 
vary between 0-8 and 0-9 as b varies from 6 to 9a. With b=7a. we finc} 
T)/Tm= 0-87, and hence T,=1-05 x 10*°K. If we denote the most probable 
energy of the electrons by Ey we should have in this case £,=0-90ev., anc} 
should therefore expect the energy spectrum of the secondary electrons to show 
a peak in the neighbourhood of 0:9 ev., while the mean energy should be in the 
neighbourhood of 13ev. In practice the majority of the secondary electrons 
appear to have energies in the range 0-2 ev. 


§5. THE DURATION OF THE 5-EMISSION PROCESS | 

From equation (2), it follows that the time taken for the temperature excess 

(Le, LT) at the centre of the distribution to fall to 1/n of its value at t=9 is 

a=(cdb?/4k)in~l). ~ Aa > 9 eee (8) 

If we set 2 —1=100 we should obtain an effective upper limit to the duration} 

of the secondary emission process; numerically this gives 
Tao OSE. 


Work on electron multipliers (Greenblatt and Miller 1947) shows that th 
secondary emission takes place in a time less than 5 x 1071" sec. 


§6. THE DISTRIBUTION PARAMETER 6b 


The results obtained thus far, using a value of 6 inferred from the columna' 
ionization theory, are in satisfactory agreement with experiment. It is a matte} 
of some interest to investigate the dependence of the results on the magnitud] 
of 6. ‘Table 1 illustrates the variation. 


Table 1 
b (a.) 6-0 6:5 7-0 Wes) 8-0 8°5 O20) 
v 32:0 ZAEO te) 4-9 23 0:99 0:35) 
Ep (ev.) (oils TOD, 0-90 0-80 0-71 0-64 0:56. ' 
T (we Sec.) 0:57 0:67 0-78 0:89 1:01 1-14 1-28 | 


v appears to vary rather rapidly with 6; the best value for b does indeed appez 
to be about 7a. The separation of Pt atoms in the lattice is approximateli 
2:5.; the uncertainty allowed in choosing a value of b to fit the experimentd 
data appears therefore to be considerably less than the lattice spacing. 
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F, in a nuclear process, gamma rays are emitted with an energy exceeding 
one million electron volts, electron—positron pairs may be created in the 
field of the emitting nucleus. The process is exactly analogous to internal 

conversion with the creation of a pair replacing the ejection of an atomic electron. 

It is well known that nuclear transitions differ in character (as do atomic 
transitions) and that they can be classified as electric and magnetic dipole, 
quadrupole etc. transitions with different selection rules for angular momentum 
and parity applying to each type. We shall show that the statistics of the angles 
between the directions of the electron and positron emitted in internal pair 
creation provides a method of distinguishing between the different types of 
transition without requiring the comparison between counts of gamma quanta 
and electrons. Magnetic dipole radiation is indistinguishable from electric 
quadrupole radiation by this method and they are in any case likely to appear 
mixed in a nuclear transition. (The selection rules that go with them are not 
mutually exclusive). 

This method offers an alternative to the well known procedure in which 
the ratios of the conversion coefficients from the K, L... etc. shells are used. 
At higher energies where the method here proposed becomes applicable, such a 
procedure may become difficult as then the conversion electrons from different 
shells differ relatively little in energy. 

External pair creation can occur only in the presence of a field of force 
(usually the electrostatic field of a nucleus) which takes up excess momentum; 
the spontaneous transformation of a photon into a pair is not compatible with 
conservation of both energy and momentum. Hence the probability of the 
production of external pairs is proportional to the square of the atomic number Z. 
In internal pair creation, on the other hand, the photon is absorbed so close to 
the emitting nucleus that in the distance it has travelled before creating a pair, 
the uncertainty principle does not allow its momentum to be determined exactly, 
or in other words, the neighbourhood of the emitting nucleus, acting as a source 
of radiation, still makes itself felt, even apart from the nuclear electrostatic field. 
As a result, internal pair creation is as important for light as for heavy elements 
(cf. Jaeger and Hulme 1935) and for the former one may assume the components 
of the emitted pair to be free. 

Such calculations were carried out by Oppenheimer and Nedelski (1933) 
and Rose and Uhlenbeck (1933). (Exact calculations carried out later by Jaeger 
and Hulme (1935) justified the approximations made in the earlier calculations 
at least qualitatively.) Rose and Uhlenbeck give results for the number of pairs 
for which the positive and negative electrons have energies E’ and E” respectively 
and are emitted in directions forming an angle @ with each other. ‘This number 
is given as a function of H’ and @ for each value of the total energy in the case 
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of electric dipole and quadrupole transitions. For the magnetic dipole transitions 

only the total number of pairs regardless of angle is given as a function of E ‘Hl 

I have carried out a similar calculation for the angular distribution in the case’ 

of magnetic dipole transitions. The number of pairs for which 6 lies within, 

an interval d@ and in which E’ lies between €’ mc? and (¢’ + de’) mc? is found to be: | 
P(<’, @)de’ sin 6d 

(<’<” +2) —3y*(e’e” +1)? 


S aoe : d p' pti —4y2(e’e” +1)] \ de’ sin dé, 
ae eter ue ey Spr i | 
“A A VoD 2 Poe J | 
where a==, 1p'l=(<?-1), |p'l=(?—1) | 


p’.p’=|p'||p"|cos8,  y=hy/me’, al 
hv being the disintegration energy and «”=y—e’ is the energy of the positron) 
in units of mc?. &’ and E” include the rest energy. 

As it is difficult to determine experimentally the total number of gamma) 
quanta emitted by an assembly of excited nuclei, it is not easy to measure the} 
absolute probability of internal pair production per emitted gamma quantum 
(Jaeger and Hulme 1935). Another type of experimental procedure would 
be to measure the number of pairs as a function of the angle 6 between the! 
two components of the pair. A com- 
parison with the calculated results, 
discussed in this paper, will then show 
the multipolarity of the gamma ray. 
In order to facilitate this comparison 
we give graphs for the three cases, 
magnetic and electric dipole and 
electric quadrupole transitions. 

Figure 1 shows, for two given 
total energies, the number of pairs 
produced per disintegration, in equal |, 
intervals of cos@, as a function of Pxs 
cos@. Equal intervals of cos@ corres- 10 
pond to constant solid angle and are 
convenient for counting experiments 8 
in which the solid angle subtended 
by each counter will usually be inde- © 
pendent of the relative position of the 
two counters. _ 

Figures 2 and 3 show the number 
of pairs produced, per disintegration, 
against cos@, for given electron— 


20 


Y¥=3 (right hand scale) 


positron energies. We have shown | os oF -05 -| 

. . cos 
only two representative CASES ee ge Figure 1. Overall angular distribution of pairs. 
which the electron carries half the 1, 2, Electric dipole. 
energy, the other in which one carries 4, 6, Magnetic dipole. 
most. It is easy to see from these two 3, 5, Electric quadrupole. 
cases the qualitative behaviour The absolute values are given for the electric 


dipole case and the other two. normalized to 


other energy distributions. Goincide aire eceiana! 
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: at } (right hand scale) 


213} (ef hand scale) 


8 
6 
4 : } ek handscale) 
2 
0 


| 05 fe) -05 = | 05 O cos? -05 -I 


cos®@ 
The probability for a given energy 
distribution with y=5. 


The probability for a given energy Figure 3. 
distribution with y=3. 
1, 2, Electric dipole. 
5, 6, Magnetic dipole. 
3, 4, Electric quadrupole. 
The absolute values are given for the electric 
dipole case and the other two normalized to co- 


incide with it at cos 0=0. 


Figure 2. 


1, 4, Electric dipole. 
3, 5, Magnetic dipole. 
2, 6, Electric quadrupole. 
The absolute values are given for the electric 
dipole case and the other two normalized to coincide 
with it at cos 6=0. 


7b 

Lastly we show, in figure 4, the 
average value of cos@ plotted against 
various energy distributions among 
the components of the pair and for 
two values of the incident photon 
energy. 

The following conclusions may 
be drawn. We have obtained results 
symmetrical in the electron and posi- 
tron by neglecting the electrostatic 
field of the nucleus. ‘The magni- 
tude of the effect is of the order 
0-5 —1:5.10-* pairs per nuclear dis- 


integration. ‘The angular distributions 
are qualitatively similar. [et us call 
Ky and Kg respectively the ratios of 
the number of pairs per disintegration 
at cos@=1 to those at cos@= —1 in 


Peo ee isc a a8 
KSA ge 3 
Figure4+. Theaverage of cos @ for y=3 and 5. 


1, 3, Magnetic dipole or electric quadrupole. 
2, 4, Electric dipole. 
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the case of electric dipole and quadrupole (or magnetic dipole) transitions respec- 
tively and let K=Kg/Kp. (We see from the curves that electric quadrupole and 
magnetic dipole transitions are hardly distinguishable by this method.) 

In figure 1, for y=5, K&8 and for y=3, K&15. In figures 2 and 3, for 
very asymmetric energy distributions between the components of the pair, 
we find K=2. When the electron and positron share the available energy 
equally, we find that K =~ 5.103 for y equal to both 3 and 5. 

In figure 4 we find that the average angle between the pair, for y=3, is the 
same for both electric quadrupole and magnetic dipole radiation. In both 
these cases the average angle, for each distribution of energy among the com- 
ponents of the pair, is about 20°, larger than in the case of electric dipole transitions. 
The same applies for y=5 when all the average angles are somewhat smaller. 

The method may be used, together with selection rules such as parity and 
angular momentum, to interpret nuclear spectra by assigning angular momentum 
quantum numbers to the various levels. 
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ABSTRACT. Experimental details are given for the determination of self-weakening of | 


B-rays in thick samples. The self-weakening correction factor in the case of twelve 
radio-elements is found over a considerable region to depend logarithmically on the mass- 
thickness of the sample. The self-weakening half-thicknesses 7} is as a rule 2:2 times 
larger than the external weakening half-thickness 4. The logarithmic dependence of the 
selt-weakening correction factor on mass-thickness and the value of the ratio 73/ €4 are 
derived in a simplified theoretical treatment. 


————— eee eee 


* The work described in this paper was done for the Department of Scientific and Industrial 


Research, and the National Research Council, Canada, during 1942 and 1943. It was reported at | 


that time, but could not be released for publication due to security regulations. 
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si) INTRODUCTION 


HE measurement of the specific activity, ie. of the number of disintegra- 
tions per unit time and per unit mass, of a given B-ray emitting sample is 
frequently involved in the exploration of neutron densities with de- 

‘tectors, and in the determination of cross-sections for nuclear reactions leading 
to the formation of B-emitters. Such measurements are usually carried out by 
placing the sample close to a f-ray sensitive device (Geiger counter or 
electroscope), which, however, fails to respond to such f-particles as may be 
-absorbed or scattered away during their passage through the sample itself. This 
“ self-weakening ”’ has to be taken into account as a correction when interpreting 
the response of the measuring device. 

The necessity for this correction can sometimes be circumvented. If the 
sample is thin enough, self-weakening may be neglected altogether, but the 
recourse to thicker samples is offen necessary in order to get enough intensity. 
Many kinds of comparative measurements, in arbitrary units, can be carried 
‘out by using samples of exactly identical thickness, so that the self-weakening 
‘corrections cancel out. ‘The preparation of samples of exactly prescribed 
thickness is, however, inconvenient, especially if the material is short-lived 
-and needs chemical treatment before counting. Moreover, this device is of 
no avail if different 8-emitters are to be compared. 

It might be thought that the ‘“‘true”’ specific activity (i.e. that of a sample 
of vanishing thickness) could be deduced from the observed specific activity 
(of a thick sample) by applying rules on f-ray absorption, using values of absorp- 
tion coefficients predetermined in suitably chosen external screens. But absorp- 
tion and scattering of B-rays are so complex, and the sample-and-counter 
geometry is so little amenable to mathematical treatment, that a purely empirical 
approach seems desirable. Broda, Guéron and Kowarski (1942) calibrated 
the observed specific activity 2 for various mass-thicknesses of several S-emitters 
(UX,, U239, Mn56 and 1128) in samples of equal area and in their particular 
counting geometry. A linear relationship between the mass-thickness and 
the logarithm of the specific activity was found to hold in a wide range, including 
in particular the domain of small thicknesses. With the help of such calibration 
-curves the true specific activity, i.e. the specific activity as measured with very 
thin samples, can be deduced from the value m measured with thick samples. 
The same authors listed values of the ‘‘ mass-thickness for self-weakening to 
one-half”’ (7;) i.e. the mass thickness at which the f-ray intensity per unit weight 
-of a given sample is reduced to one-half. Grummitt, Guéron and Wilkinson 
(1943) investigated the same elements in a different counting geometry. ‘They 
-confirmed the general conclusions drawn by Broda, Guéron and Kowarski, 
but their values of 7; were slightly different, perhaps because of the changes 
in the geometry. Grummitt, Guéron and Wilkinson also extended the method 
to Th233, Bal39, Bal40, La140, and 1I131. Results of other authors 
(Hendricks, Bryner, Thomas and Ivie 1943, Yankwich, Rollefson and Norris 
1946, and Cook, private communication) on $35, C14 and Ms Th2 will be 
-shown to be amenable to a similar treatment. While the method, therefore, 
is widely applicable, the numerical values quoted in this paper are valid only 
under the given experimental conditions. 
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§2. EXPERIMENTAL PROCEDURES 
(i) Chemical 


(a) Mn56 and 1128. ‘These activities were obtained by exposure to slow 
neutrons of manganese oxide (empirical formula MnOj,,;) and lead iodide.. 


The activity of these compounds was measured without chemical treatment 
after irradiation. 


(6) U239. Uranium as sodium uranyl acetate was irradiated with slow 
neutrons. In order to reduce the natural B-ray activity, and to remove fission 
products, chemical treatment was necessary. Excess sodium acetate was added 
with stirring to a cold solution containing 10°% uranium by weight and 25% 
acetic acid by volume. The crystalline precipitate of sodium uranyl acetate 
was washed on the filter with saturated sodium acetate solution, a small amount 
of water, and finally with acetone. It was then dissolved in the minimum 
amount of nitric acid, the concentrations of uranium and acetic acid were adjusted 
and the uranium reprecipitated. Two sodium uranyl acetate precipitations 
before irradiation reduced the UX activity almost to zero; after irradiation, 
a further two precipitations reduced the fission product activity and newly 
grown UX activity to less than 1% of the initial U239 count. ‘The samples 
were ready for counting 10-15 minutes after the end of irradiation. 


(c) Th233. Thorium also had to be freed from its decay products. The 
procedure was to add to a solution of thorium nitrate before irradiation barium 
and lanthanum carriers (for radium and actinium) and to precipitate their 
carbonates with excess ammonium carbonate. ‘The filtrate, containing the soluble 
thorium carbonate complex, was acidified with nitric acid, and after addition 
of more barium salt barium*was precipitated again with sulphuric acid and 
removed by filtration. Now lanthanum, bismuth, lead and barium ions were 
added as hold-back carriers, and the thorium was precipitated at 60° c. with 
excess 40°, H,O, as thorium peroxynitrate. This was redissolved in the 
minimum quantity of boiling concentrated HNO, with the addition of a little 
HCl. ‘The peroxynitrate precipitation, again in presence of a hold-back carrier, |} 
was repeated, and the precipitate ignited to ThO,. The residual f-activity |} 


of non-irradiated ThO, was about 40 counts per minute per gramme (counter 
efficiency about 25%). 


(d) Bal40 and Lal40. These activities were separated from the mixture |] 
of fission products formed in the slow neutron irradiation of uranium oxide 
or uranyl nitrate. The details of the chemical procedure (Grummitt, Guéron, 


Wilkinson and Yaffe, 1947) follow in general the fuming nitric acid method 
of Hahn and Strassmann (1942). 


(1) Preparation of samples and counting 


At first cylindrical aluminium counters of wall thickness 40 mg/cm? were | 
used. ‘The counters were filled to 4cm.Hg with argon and to 2cm. Hg with | 
alcohol vapour. he finely powdered activated materials were spread evenly | 
on rectangular aluminium trays of area 2-1cm?. A few drops of an acetone | 
solution of collodion (~0-5 mg/cm? collodion) were added as binder, the sources | 
dried by gentle heating and covered with a single layer of “sellotape”. The | 
sources were disposed around the central portion of the counter tube. This | 
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method allows efficient utilization of weak sources due to the large area exposed. 
It is, however, an inconvenient and inaccurate arrangement for work with 
external absorbers. 

In the later measurements (Grummitt, Guéron and Wilkinson 1943) a 
“bell jar” type of B-ray counter with a mica window of area 7-1 cm? and mass 
thickness 5 mg/cm? was used. It was filled to 9 cm. Hg with argon, and to 0-5 cm. 
Hg with alcohol. Circular aluminium trays of area 7-1cm? were supported 
1cm. below the window by means of a sample holder placed in a slotted frame. 
Samples were prepared as previously, but were not covered with “sellotape”’. 

In no case were the samples counted without mechanical or chemical mixing 
after irradiation. This is important as the effect of unequal activation at different 
depths, due to the absorption of the effective neutrons, particularly resonance 
neutrons, has to be excluded. 


(iii) Measurements 


The measured activities of the samples were corrected for variations in the 
counter efficiency with the help of uranium standard. In the cases of U 239 
and ‘Th 233 an additional correction of the measured f-ray counts was necessary 
due to the growth of the natural f-ray activity in the samples after chemical 
purification. Correction was made with U239 by counting the UX activity 
some time after the decay of the 23-5 minutes period with extrapolation back 
to the time of counting of U239. With Th233 a non-irradiated reference 
sample was used; since the fission products of Th were too weak to interfere, 
chemical treatment after irradiation was not necessary. 


The results were expressed in ———— il 
| 


terms of n, the apparent specific 
activity of the material, i.e. as the 
number of counts per minute per 
gramme of active material at an 
arbitrary reference time. ‘The log- 
arithm of m was then plotted against 
the mass thickness M of the sample. 
In figure 1 typical plots for U239 
are given. As the neutron flux 
varied between the sets of experi- 
ments, the results of each set were 
normalized at a point in the middle | 
of the mass-thickness range studied. ie Sevan Dic aoe ate ale 
Extrapolation of this composite MASS THICKNESS  mg/em* 

curve to mass thickness zero gives a Figure 1. __ Self-weakening of U239 B-particles 
value mp, the true specific activity stood tine nan yUaeeite: 

of the sample. In figure 2 the . 
self-weakening calibration curves, i.e. the plots of log (9/n) against M, are given 
for the bodies investigated. The self-weakening correction for any thickticss 
is read directly from these curves. For all the radiations measured log ny/n 
depended linearly on M. In table 1 the values of the self-weakening half- 
thicknesses 7;, as read from the calibration curves, are collected. Further, 
the table gives the weakening half-thicknesses «, as determined by measurement 
of the B-ray absorption in aluminium foils placed midway between sample and 


APPARENT SPECIFIC ACTIVITY 


464 E. Broda et al. 
Table 
Radioactive nucleus Active material Nk Eh n3/€4 
(mg/cm’) — (mg/cm?) 
UXs Sodium uranyl 243° 135 1:81 
acetate DS 
U 239 (23:5 min.) Sodium uranyl 1132 540 Pri 
acetate 2 
Mn 56 (2:54 h.) MnO}j.47 458? 171 2-74 
p22 
1 128 (25 min.) PbI, 266? 128° 2:08 
2O1e 
Th 233 (23 min.) Tho, 134 56° 2:4 
La 140 (40 h.) LaF; 156? 76° 2:06 
Ba 140 (12-7 d.) BaSO, 918 416 223 
5:35 (87-1 d.) BaSO, 6°9¢ 3-14¢ Dee 
C 14 (4700 y.) BaCO, 6-04 34 2:0 
1 131 (8-0 d.) Cul 568 22° 2°54 
602 De 2:4 
Ba 139 (86 min.) BaSO, 3232 1475 22 
Ba(NOs)s 2952 
MsTh 2 (6:13 h.) BaBeF, 145f 


* Broda, Guéron and Kowarski. 
> Grummitt, Guéron and Wilkinson. 


¢ Hendricks, Bryner, Thomas and Ivie 1943. 


counter. It is to be noted that the 
external absorption for UX,, Mn 56 
and 1128 in the present geometrical 
arrangement does not follow an 
exponential relationship over even 
a limited portion of the range. 

Figure 2 and table 1 also give 
self-weakening data calculated from 
experiments by other authors to 
conform with the logm/n against 
M presentation as used above. It 
is noteworthy that even with p- 
radiations as soft as those of S35 
and C14 the logarithmic relation 
holds over a considerable range of 
mass thicknesses. 

In table 1, column 5, the ratio 
of the self-weakening half-thickness 
and the. weakening half-thickness 
(with external screens) is listed. 
With a few exceptions, notably UX, 
and Mn 56, the values of this ratio 
fall close to the average ,/€, = 2.2. 


$3.) DISCUSSION 


The experimental results show. 
that the relationship between the 
self-weakening correction factor and 


@Yankwich, Rollefson and Norris 1946. 
¢ Broda 1943, 1944. 
f Cook, private communication. 
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Mn56 in manganese oxide (MnQ}.,7). 
Ba139 in barium sulphate. 

1128 in lead iodide. 

UX, in sodium uranyl] acetate. 

La140 in lanthanum fluoride. 

Th233 in thorium dioxide. 

U239 in sodium uranyl acetate. 

S35 in barium sulphate (upper scale). 
C14 in barium carbonate (upper scale). 


The self-weakening correction factor m/n vs | 
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the mass-thickness of the active sample can be considered as logarithmic over 
a considerable mass range. 

The relation can be deduced by a simplified theoretical treatment with the 
following assumptions : 

(a) ‘The measurements are made in a parallel beam geometry. 

(6) ‘The weakening of f-rays as measured by external screens obeys an 
exponential law over the initial portion of the curve, and the weakening is inde- 
pendent of the nature of the screen, i.e. 2 =m,e~"” where x is the mass thickness 
of the absorber, » the mass absorption coefficient, and m, and m the numbers of 
incident and transmitted particles. The external half-thickness €, 18 0-693 /. 

(c) The back-scattering of B-rays by the sample and by the backing material 
can be considered as identical. This assumption is justifiable if backing material 
of atomic number similar to that of the active material is used. 

Then the contribution to the total f-activity of a layer of mass thickness 
dm at a depth m in a thick sample of active material of specific activity 1 is 


GAN sant 
The total activity of the sample of mass thickness M is, therefore, 
-M 
A= | nye" dm =(no/p)(1 —e *™). 
Jo 
The apparent specific activity ” of the sample is 
n= A/M=(no/uM)(1 —e7"™) 
and the ratio of the apparent and the true specific activity is 
Wins Cols) a pM p?M? peMeB 


a hs AG 16 Coie er ad 


In figure 3 the expression loguM/(1—e~"™) is plotted as a function of M. 
Thus it appears that on the above assumptions the self-weakening correction 
curve ,/n against M is approximately logarithmic over the initial region. The 
value of uM, which reduces the activity function (1) to 4, is 1-53, ie. un, =1-53 
and 4, = 2-21 «,. 

Serious divergence from the linear 
relationship occurs at values of 
pM>~2, ie. at values M>~3 «. 
The experimental correction curve 
for C14 in barium carbonate shows 
divergence from linearity at mass 
thicknesses greater than 8 mg/cm’, ie 
ie. approximately 2-7 times the 4 
external absorption half-thickness for 
the f-rays of C14. However, in the 
case of S35 no deviation is evident 


upito dW =4-5 «;. 
Feather (1943, 1944) has sub- eee a a 
i self- u 
mitted the whole problem of a Figure 3. The theoretical self-weakening 
weakening of f-rays to a caretu relation as given by a simplified 


theoretical treatment, and has given treatment. 
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explicitly the conditions under which a logarithmic relation, as found in our 
work, may be expected to hold. 

It is to be expected from the well-known complexities of f-ray scattering 
and absorption that the numerical values will vary with the geometrical conditions 
of measurement. However, the relation »,=2:2¢, is approximately obeyed 
in the majority of the experiments listed in table 1. Such discrepancies as 
exist are not accounted for by the differences between the average atomic number 
of the active material and the external absorber (Al), nor is there so far any 
obvious correlation with the energies of the radiation. Data on radiations 
from more radioelements are necessary before drawing any conclusion. 
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ABSTRACT. 'The energies of the 6 rays of thorium C’’D have been measured using a | 
magnetic spectrograph with semi-circular focusing. Corrections have been made for | 
contamination of the source by ThC and ThB deposited by the processes of B-recoil and | 
aggregate recoil. The absolute intensities of the stronger lines have been measured and the | 


shape of the continuous spectrum has been compared with theory using three recently 
proposed energy level schemes for the product nucleus Pb?°. 


These level schemes fail to explain the observed shape, partly because they omit the | 
P y-ray (859 kev) which can now be assigned to ThC’D. Evidence for the 8-excitation 


of energy levels at about 4:5 Mev is discussed. 


S INTRODUCTION 


EVERAL level schemes have been proposed for the nucleus ThD. 
Among the more recent are those of Oppenheimer (1936), Arnoult (1939) 
and Itoh and Watase (1941). It is of interest to compare the f-ray spectra 

predicted by these schemes with the f-ray spectrum of ThC’D found by 


experiment. 
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The short half life of 3-1 min. of ThC” indicates that the B-disintegration 
belongs to the class of allowed transitions. The spectrum is the sum of a number 
of partial spectra having different end points and intensities and each exciting 
a different level of the product nucleus ThD (Pb”8). All the stronger partial 
spectra should be allowed transitions and should have the shape predicted by 
the Fermi theory. 

The shape of the experimental spectrum found in the present work deviates 
from those predicted by the schemes hitherto proposed for the energy levels 


and y-ray intensities of the product nucleus Pb2°°. Much of this deviation may 


be due to the omission from the schemes of the 859 kev y ray whose assignment 


of the ThC’D disintegration has been confirmed by our measurements. 


A surplus of low energy f rays suggests that one or more levels of Pb?’ at about 
4-5 MeV receive f-excitation. 

Caution is required in the interpretation of the observations because of the 
presence of ThB and C nuclei carried over in the preparation of the ThC’ 
sources by «-recoil. Corrections for this contamination have been subtracted. 
The thin mounting used for the source should not have distorted the spectrum 
above 100 kev. 

- Alichanian and Zavelsk1 (1937) found relatively more f rays at both high and 
low energies than were found in the present work. 


$2, THE EXPERIMENTAL, METHOD 


The sources of 'ThC” were deposited on the thin aluminium leaf A (figure 1) 
by «-recoil in vacuo from a source of thorium (B+C+C”). The latter, which 
may be called the activat- 
ing source, was with- 
drawn when counting. 

The short half life of 

ThC’’ necessitated the jue eri i 
use of many recoil 
sources, whose intensities 


were correlated by Late : 
measuring their y-ray © sb 3, 


activities with an inde- 
pendent y-ray counter in 
a fixed position. 

The deposit of ThC” 
was found to be con- 
taminated with thorium 
Band C, which also emit 
B rays. This contamina- 
tion was due to the pro- 
cesses of aggregate recoil 
and f-recoil. In the 
former, anatom of ThC” 
recoiling strongly from 
the emission of the « ray 
fi the, ThC.C” decay, 


carries with it a random Figure 1. The magnetic spectrograph. 


| 

| 

| Figure lb. Details 
| of the activator 
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sample of neighbouring atoms from the surface of the activating source. In the 
second process a small fraction of ThC atoms leaves the surface of the activating | 
source owing to the weak recoil momentum received from the emission of the | 
decay electron and neutrino from ThB. | 

The aggregate recoil effect was the more serious, and its magnitude was 
followed in each day’s run by repeated measurements of the peak of the very | 
strong F B-ray line (1386 Hp) of thorium B.C. ‘The magnitude of this effect | 
showed wide variations on different days. 


§3. THE MAGNETIC SPECTROGRAPH 

During the preparation of the ThC” sources recoil particles may be deposited | 
on other surfaces in the spectrograph in addition to the source mounting A. 
By closing the shutter S (figure 1) during the time when the activating source 
was near to A, recoil particles were prevented from settling near the main defining | 
slit S, or from passing through it. If S was left open during the process of | 
activation the background counting rate of the 6-counter showed erratic variations. 
These were probably due to charged f-recoil atoms of low momentum deviated | 
through the slits by the magnetic field. These atoms of ThC contribute to the | 
count at zero magnetic field when they disintegrate near the counter. 

The baffle system shown in figure 1 was designed to exclude from the counter |] 
B rays emitted from these deposits on the surfaces of the spectrograph and also 
scattered B rays. rays coming from atoms on surfaces other than A and lying 
on the source side of the main slit S, were excluded from the counter slit S, 
by the second defining slit S,._ The surfaces C and D did not, however, lie in 
the zone of exclusion defined by the slits S,, S,, and S;._ 6 rays from the surface D 
were excluded. by adding baffles B, and B, of the type used by Neary (1940). 
B rays from the surface C and from the edges and leit sides of B, and B, were not 
excluded, but scattering from them was reduced by inclining them so that only 
their edges were exposed to the source. ‘This should also keep them free from |f} 
a-recoil particles, which, owing to their great momentum, travel in straight lines. |]} 
in the magnetic field. The bafHle B; was designed to keep secondary aggregate 
recoil particles (coming from A when the shutter S was open) away from a region 
where the thickness of lead screening the f-counter was small. . 

The aluminium leaf A was 9 mm. wide and 3-5 cm. long and had a weight |f] 
of 0-13 mg/cm?. It was inclined at 22° to the central ray of the focused beam |} 
of radius 9 cm. ‘The depth and the greatest width of the beam at S, were both | 
lcm. A high vacuum was maintained by using a 2-stage oil diffusion pump. 

The plane pole-faces of the magnet were rectangular in form of sides 19 x 28 in. 
and were 5-9 cm. apart. ‘Their edges are shown by the broken line in figure 1. 
‘The magnet showed the phenomenon of magnetic viscosity and, in order to get 
a reproducible relation between magnetizing current and field strength, readings 
were always taken on the ascending branch of the same hysteresis cycle which || 
had a constant maximum field of 1400 oersteds. The field was calibrated by | 
using a flip-coil and ballistic galvanometer. To get an absolute scale of Hp we | 
assigned to the field settings at the centroids of the f-ray lines the values of Hp 
found for the lines by Ellis (1932). This procedure, discussed by Lawson and || 
Tyler (1940), is correct if the curvature of the continuous spectrum is small. | 
(The shape of the F-line of thorium B.C is shown in figure 2 (a). The width 
of the line at half its height is about 1:35°% of H. The arithmetic mean or centroid |} 
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of the line shape lies about 0-2°% of Hp to the left of the middle of the two points 
half-way up opposite sides of the line.) 


§4. THE s- AND y-RAY COUNTERS 

The f and y Geiger-Miiller tubes (G and M in figure 1) were filled with 
mixtures of argon and alcohol vapour. The cathode of the B-counter was of 
brass 2 cm. long and 1 cm. diameter. The anode was a tungsten wire 0-1 mm. 
in diameter. ‘The reservoir R ensured constancy of the gas mixture of 6 cm. 
of argon with 2 cm. of alcohol. ‘The slit S,; was 2-4 mm. wide and 1 cm. high 
and was covered with mica. 

The use of the y rays to correlate the sources requires that the ratio of the 
efficiencies of the two counters must remain the same over an appreciable time. 
A source of error is introduced by the finite slopes of the plateaux of the character- 
istic curves of the two counters... The f-counter had a plateau about 60 volts 
long and a slope of about 0-15% per volt at 930 volts. The y-counter had a slope 
of 0-4 per volt at 1500 volts. The voltages were stabilized by small neon tubes 
and were kept as constant as possible with respect to the lower end of the 
plateaux. 

On each day the weak ends of the spectrum were counted first, the middle . 
parts being left until the source was weaker. In this way a day’s run roughly 
covered the whole spectrum so that day-to-day changes in relative efficiency 
should introduce no systematic distortion. ‘This method of sandwiching the 
measurements should make random the error due to small changes in voltage. 

The dependence of the efficiency of the B-counter on the energy of the 6 rays 
was found by measuring the spectrum of thorium (B+ C+C”) in equilibrium, 
at different pressures in the counter. ‘The efficiency at high energies (935 kev) 
did not begin to fall until the pressure of argon was below 6 cm. with 2 cm. of alcohol. 
It was concluded that at least 98°, of the 6 rays of ThC’’D were counted at all 
energies. 

The loss of counts due to the finite resolving time of the counting systems 


_ was considered to be negligible. Each counter fed into a Kipp-relay equalizing 
circuit followed by a pentode scale-of-eight circuit due to O. R. Frisch. The 


maximum counting rates were about 1200 per minute. 

The absorption of slow f rays in the f-counter window was estimated by 
measuring the spectrum of Th(B+C+C”) and comparing it with the results 
of Flammersfeld (1939), who used a much thinner window. ‘The surface density 
of the mica was about 2 mg/cm?, 6 rays begin to penetrate at 600 Hp, while the 


- intensities of the B and Bb lines at 653 and 678 Hp were respectively 19° and 30% 


of Flammersfeld’s values. At 1000 Hp the difference Snes the two 


_ spectra was only about 10%, a loss which is much less than that found, for example, 


by Martin and Townsend (1939). This is probably because the part of the 
counter slit outside the mica window was narrower than the part inside it, so that 


rays could suffer appreciable sideways deflection and still be counted. The 
conditions resembled those in Schonland’s experiments (1925) and by fitting 
his law of absorption for aluminium to the present measurements we estimate 


the fraction of the incident beam lost in the mica window to be 10%, 2:6%, 

-and 0-6% at Hp 1085, 1542 and 2186, respectively, These corrections have not 
been included in figure 4. They would increase the downward convexity of 
that end of the spectrum. 
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§5. PREPARATION AND CORRELATIOW OF THE SOURCES | 

The activating source was obtained by overnight exposure to thoron of a 
platinum disc P, 4-5 mm. in diameter, which was then screwed into the “ activator” 
shown in figure 1b. A high a-recoil efficiency was obtained by polishing le with | 
chromic oxide and exposing it to the thoron in air at 6 cm. pressure in order to 
reduce tarnishing. 

The hood H was pushed on to the activator to prevent the sideways escape, 
of recoil particles from P, which, if they fell near the glass window W of the | 
spectrograph, would give a disproportionately large y-count. The gap between. 
H and A was about 0-5 mm. 

In order to reduce the effect due to y-ray sources not on the leaf A, the 
y-counter was entirely surrounded by a lead shield 7-3 cm. thick, whose only open-J 
ing was the collimating channel E, 2:7 cm. wide and 0-8 cm. deep at its outer end.| | 

The usual working cycle was 14 minutes ; of this 6} minutes were used for. 
counting and the rest for activation. In order to obtain a reliable measure of the} 
source strength it was arranged that the y-ray count was more than 2} times the: | 
B-ray count, even at the peak of the continuous spectrum. 

During counting the activator was withdrawn about 1:3 metres down the} 
curved tube T by means of a sliding magnet acting on an iron rod I. When 
withdrawn an additional 24 cm. of lead separated the activator from the counters 
making the effect of its y-radiation negligible. 

Tests with a fluxmeter showed that the motion of the activator magnet altered 

‘the magnetization of the main magnet and that owing to magnetic viscosity theft 
latter took about 1 minute to reach a steady state. It was therefore necessa 
to go round the standard hysteresis cycle during each activation. 


SiG Wess, CORRECTION FOR AGGREGATE RECOIL 

Aggregate recoil is the only recoil process which may deposit thorium B o 
the source mounting A. It may thus be studied separately by measuring the 
strong F line which contains more than 0-25 electrons per disintegration and 
stands out about 30 times higher than the maximum of the thorium B continuou 
spectrum at the resolving power of our apparatus. ‘The growth in each day’g 
run was followed by making about four measurements of 8, the counting ratdy: 
at the peak of this line. 

The counting rate 8, due to aggregate recoil at a given setting of the magnetit 
field H can be calculated if By and B,/By are known. Using small greek letter} 
for rates and capitals for total counts, we have | 

B=B’+B,+Bp,+By 


= BY +[Bo+ Bs + Be(Pa/Bp) At, = 9 oes (1 a) 

where B is the total P-particle count in a counting period of length a (usual; f 

63 minutes) ; B’ is the count due to ThC’’D with half-life 3-1 min. ; , is thy) 

natural background rate in the f-counter, and fB, is the rate due ‘6 B- recoil} 

By decays with the 1-hour period of Th(C+C”) and thus may be consider 

constant during Af. 1 

The contribution [', of the aggregate recoil to the y-count was also pro 
portional to By. We have 

D=0’4+T94+T gtr, 

=I" + (yo +yp+0-038,)At, seeeee (1 5) 

where 0-03 is the experimental value found for y4/By. 


The B-ray spectrum of ThC'"D 471 


1000 


Counts per minute 


0 L I pie ear PLS yaw Oi | t 
0-309 0-313 0-317 0-347 0350 078 0:79 0-80 
Magnet Current 


Figure 2. f-ray lines. The F line of ThB.C by aggregate recoil and the G and P lines 
of ThC”. D by a-ray recoil. The broken lines indicate the continuous spectrum. 


From these equations B’/I’’, the standardized B-ray intensity of ThC’’D, 
was found for each value of H. ‘The measurement of 6, was made easier by 
| the flatness of the peak of the F line, as seen in figure 2. The flat part had a width 
| of about 0:5% of Hp and the field was increased through this zone in steps of 
about 0-2°%. [was taken as the highest reading lying between two lower ones, 
all three containing more than 2000 f-counts. This choice of the highest of 
three readings subject to fluctuation might give rise to asystematicerror. Analysis 
showed that the small overestimation of 6, from this cause was, by a happy chance, 
| nearly cancelled by the mean fall in the highest reading due to its distance from 
the exact top of the line. ‘This,mean fall increases with the spacing of the field 
| settings and the curvature of the peak. 
The aggregate recoil efficiency fell during the course of each day, in agreement 
' with the findings of Lawson (1919). It showed big variations from source to 
/ source. In order to 
| find B’/I’”’ by equa- 
tion (1), it is necessary 
to know f, during 
each counting period 
At. fy wasestimated 
“by the method of 
_ graphical interpola- 
_ tion shown in curve 
(1) of figure 3:, where Time of activation of source 
| B,, and A,,,, are values 
of By» found from the 
nth measurement of °° ene 4 6 
the F-line. In a group Bours 
_of activations running 


Time of measurement of F line 


Figure 3. Growth of contamination by aggregate and f-recoil. 
(1) Counting rate x 10% at the peak of the F line of ThB.C. 


_for le from A 
9 a : (2) nJne, the fraction of the equilibrium strength of the B- 
_to B,, a little more recoil contamination. (3) (n/ne)e~“st, proportional to its 


than half the time was absolute strength 7. 
31-2 
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spent in activation so that the mean rate of growth of By during activation was 
just less than twice the slope of A,B,. Since, however, this rate was always falling, 
a good compromise can be ieade ie taking the rate during the first activationy 
A,C, to be twice the slope of A,B, and then assuming a linear growth from 
C, to B,. The values of By are then read off along C,B,, C2B, etc. i 

The wales (8/8) in equation (1) were found in two steps. First, one value 
of this ratio was found by counting at the maximum of the aggregate recoil 
spectrum on days when it was strong, applying a correction for B-recoil. The 
second step was to fit this known value of the ratio at Hp=4541 to the known 
spectrum of thorium (B+ C+C”) in equilibrium on a thin aluminium leaf, from 
which other values of {,/8, could then be read. ‘The stepwise compariso 
was necessary because the relative height of the F-line above the continuous 
spectrum was not necessarily the same in the two spectra, the height and shape} 
of a f-ray line being sensitive to changes in the position and distribution of the 
source on A. No such changes occurred during the measurement of ThC"D 
because activation by recoil does not disturb A and seems to give a constaty 
surface distribution. 

The errors in the values of UB, given in table 1 arise (i) from the use of only 
2000 counts to find fy, (ii) from the graphical interpolation, (iii) from the valuéf 
of (B4/By) at Hp =4541 and (iv) from the spectrum of Th(B+C+C”) in whic | 
about 30 000 counts per point were used. It is thought that the errors are coverec di 
by assigning a standard deviation of 3°% to XBy. 


' 


Table 1. The continuous spectrum, the corrections for aggregate 
and f-recoil and values of the Fermi function F(Z, 7) 


He SB 2B, -2Ba Bas BB/EI” aN) gaa 
0 1,542 1017 0) 0 0-0078 — — | 
1085°5 22,455 2756 6254 306 0-0808 2-1 3-208} 
1541-6 32,261 2835 8548 547 0:1489 1197 4-906,| 
2185-8 21,228 1597 DUBS 468 0:2390 ley 76784) 
2770-1 17,906 IGS 2049 433 0-3082 1123} 10-667,| 
3306-7 20,839 1169 2847 635 0-3496 iles 13°78 
3843°5 19,674 1455 2443 595 0-3714 1°3 17-203 
4541-2 18,259 1176 2163 603 0-3313 13) 22°11 
5077-9 20,582 1409 2840 813 0-2663 1-3 26:235 
5594-7 25,767 2096 4579 1267 0:1787 ihoS 30°45, 
6128-9 23,834 2702 4230 1433 0-1104 1:8 35-085. 
6662-4 16,497 2847 2992 1151 0-0513 2°5 39°95 
6928-6 7,894 1776 1433 519 0:0319 S97) 42-47. | 
7194-4 2,578 1015 462 228 0-:0132 16:3 45-05, 
n=H p/1704-2 
B= Total number of f-counts. 
SB y= Natural background. 
2Ba=Counts due to aggregate recoil. 
»Bs=Counts due to f-ray recoil. 


2B” /=I’’=Standardized f-ray intensity. 
%o(N)= Percentage standard deviation of N (cf. eqn. (5)). 


$7) THE CORRE CRION THOR G=RE COL 
The counts assigned to this effect are shown as UB, in table 1. 


_ The counting rates are so small that in contrast to aggregate recoil no dire 
method is available for finding the growth of the contamination at intervals duri 


( 


f 
1 
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the day. For simplicity we assume that the ratio of the «- and f-recoil efficiencies 
is constant for deposition on A from the activator. The value of 6p in equation 
(1 a) can then be estimated for each counting period At. 

The first step is to assign to each At a measure of the relative strength of the 
f-recoil source at that time. This measure is given by n/m,, the ratio of the 
number n of B-recoil atoms on A to the equilibrium number », which would 
have been reached if the interval had been preceded by a very long period of 
activation. 

Let us assume that a constant fraction 8 of the ThB atoms decaying on the 
activator deposits recoil atoms of ThC on A. Then at time ¢ during an activation 
0 


n=neW tl) 4B vores 7 ge Cea), bia mates (2) 


where n; is the initial number of recoil atoms n on A at the start of the activation 
at t= £;, mp is the initial number of ThB atoms on the activator at ¢=0, and A, and A, 
are the decay constants of ThB and ThC. Between activations n decays with 
decay constant A,. The equilibrium value n, is given by 


Peace ee ee Ga cena (3) 


Using (2) and (3) average values of n/n, were estimated for each counting 
period At. An example of the serrated curve of m against ¢ is shown in figure 3 
for a typical day’s run. 

The second step is to find the ratio of the «- and B-recoil efficiencies. Defining 
the f-recoil efficiency f as in equation (2) and the «-recoil efficiency « as the fraction 
of ThC.C” disintegrations which deposits ThC’’ on A, let us consider the case 
of an activation so long that both y” and yp, have reached transient equilibrium 
and n/n, is unity. ‘Then every N atoms of ThB decaying on the activator will 
deposit V8 atoms of ThC on A by f-recoil and 0:34 Na atoms of ThC” by «-recoil 
(because only 34% of ThC atoms decay to ThC”). 

If we neglect the weak y rays of ThC.C’ and ThC.C” the y-activity on A 
will be due to ThC”’. D alone and will be the sum of the counting rates y”’ and yz 
due to the two modes of arrival by recoil. ‘The ratio of the numbers of ThC” 
atoms arriving by each mode (neglecting the 10°% difference between secular and 
transient equilibrium) is 


034Nx a mi) a 
NB x 034 p ve) caittiet |  See ee 


The ratio of y’’ to yz at equilibrium was found by measuring the decay of the 
activity of strong recoil sources made by long activations. ‘The decay of the 
f-activity at Hp =4541, and of the y rays, was followed for several hours and was 


analysed into three half-periods of 3-1 min., 60-5 min. and 10-6 hours due re- 


_ 


spectively to «, B and aggregate recoil. Extrapolating back gave yz and y”, 
but owing to the weakness of the effect large statistical errors in yp were present. 
Seven measurements of (y’’/77g)equi, 00 different days gave a mean of 51 withamean 
residual of 15°. Five earlier measurements of the decay of the y-activity, 
each at the end of a day’s run, gave the value 64 with 6°%, mean residual. As the 
value of 8 is known to be sensitive to changes in the emitting surface layer it is 
possible that the difference is significant and that it arose from a change which 
was made in the method of making the source. ‘The earlier value was therefore 
used for measurements taken before the change was made. 
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The value of yz was calculated for each counting period AZ, obtaining y | 
in (4) from the observed y-count I’, and using the assigned value of n/n,. ‘The | 
value of the f-recoil correction fy, in equation (1 a) was found by multiplying yz | 
by the value (8,/y,) appropriate to each value of Hp. ‘This ratio was obtained | 
by direct measurements of the B- and y-counting rates of 1-hour half- life at | 
Hp =4541. Its value at other fields was estimated from the known shape of the | 
B-spectrum of Th(C+C”) deduced from that of Th(B+C+C”) by subtracting 
from the four points of lowest energy the contribution due to thorium B. This | 
contribution was estimated from the curves of Flammersfeld (1939), making | 
a small correction for loss of ThC”’ by «-recoil from the aluminium. | 

In view of possible variations of f-recoil efficiency, a standard deviation. of | 
10°, is assigned to XB, in table 1 which makes it the major source of error at the | | 


four highest values of Hp. 


§8. THE p-RAY LINES | 

Flammersfeld (1939) has found intensities of the f-ray lines in the spectrum | 
of Th(B+C+C”) which are higher than those previously accepted. The | 
increase arises. in two ways. First, the absolute intensity of the F-line (1386 Hp), | 
used ordinarily as a standard, was found to be 0-284 electrons per disintegration | 
as against the former value of 0-250 (Gurney 1926). Second, the intensity of | 
the lines of ThC’’D was increased by about 30% by covering the source with | 
a thin film to prevent the escape of ThC” atoms by «recoil. This loss by recoil 
may be expected to occur from all uncovered sources which are. moderately | 
clean. 

It is interesting to check these intensities by a method which is independent 
both of the F-line strength and of loss by «-recoil. A second point of interest is to 
verify the assignment of some of the many lines of doubtful origin. A ThC’D |} 
line will be about three times as strong in our case as it is in the Th(B+C+C”) 
spectrum, whereas other lines will have reduced intensities set by the amount 
of the aggregate and B-recoil contamination. 

Approximate values for the intensities of the lines are given in table 2. The lf} 
absolute intensity of each line in electrons per disintegration is equal to the ratio |] 


Table 2. Line intensities in electrons per disintegration 


B-ray line Ec Fb G Jb Jb3 Ils, M N P.. Xait 
Hp (Ellis) 1251 1469 1594 1929 1980 2603 2887 3164 3924 10000 |] 
Abs. Int. in 
KC AB KO? OES WE = Dao) Oras O-0S Pale 1:42 0:38 0-26 0-18 | 
Int. in ; 
Th(B+C+C”) — 0°:356 2°73 0-416 — 2:26 1:75. 0:356 0-255 0-234 
(Flammersfeld) (Jb+-Jb1) 
Gamma ray Ec Fb G G G L L+M M 12 x 


of the area due to the line to that due to the disintegration electrons in the] 
momentum distribsion curve. This curve is obtained by plotting N against|] 


Hp where 
Be B” 
reo alr a -|>|) bre (5) 


[B’’/I’’], is the standardized B-ray intensity when H=0 and is due to the counts) 
caused by the y rays of thorium C” in. the B-counter. The area beneath each 
line is equal to the product of the height by the mean width. The latter wasif 
found to be 0:0137 Hp by measuring the area of the F-line shown in figure 2. 
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Flammersfeld’s values are given for comparison in the table, taking the 
thorium C branching ratio to be 33:7°% (Kovarik and Adams 1938). The in- 
tensities in row 2 are subject to error because it was necessary when measuring 
a ThC’’D line to reset the magnetic field after each activation. It was found 
impossible to reproduce H to better than about 1 in 300 after going round the 
standard magnetization cycle. It was thus difficult to build up a large number 
of counts on the flat peak of a line by repeated activations at the same H. A low 
B-count due to just missing the peak (about Hp/200 wide) is indistinguishable 


| from a statistical fluctuation. This uncertainty did not arise in the measurement 


of the F-line of thorium B because there was then no magnetic disturbance due 
to motion of the activator. 

In view of this limitation of the method the ratios of the intensities of the 
G, L, and M lines were found independently from the spectrum of Th(B+C+C’) 
in transient equilibrium. The lines now stand out less above the continuous 
spectrum, but it is possible to get about 20000 counts at each peak. 

The relative intensities of the lines are compared with the results of Flammers- 


feld (1939) and Oppenheimer (1936) in table 3. 


Table 3. f-ray line intensity ratios 


TSG: M:G XM eG Source 

0-826 0-641 Q-0857 Flammersfeld (1939) 
0-429 O=353 0-106 Oppenheimer (1936) 
0-709 0-597 0:0711 ThC”’D data in table 2 
0-692 0-619 — Th(B+C+C”) 


Special weight must be attached to the L:G ratio of 0-69 found using 
Th(B+C+C’) because it is the mean of three independent measurements. 
Tt seems likely that the correct ratio lies between the values of Flammersfeld and 
Oppenheimer. 

The uncertainty in setting the field can be seen from the typical plot of the 
G and P lines in figure 2. Each point is the value of B’/I” found from a single 
activation. B’ is the B-ray count not corrected for f-recoil nor for the counts 
due to the line itself in the aggregate recoil spectrum. ‘This correction, when 
made, improves the agreement between data taken on different days. The purely 
statistical standard deviation in B’/I” due to the finite numbers counted is about 
1-8°%, for the G line in which about 10000 f-counts and 12000 y-counts were 
recorded near the peak. In the case of the P line this error is about 10° and it 
‘is higher for the weaker lines. It is likely that the error due to just missing the 
peak was larger than the statistical error for the stronger lines. It would be 
expected to lead to an underestimate of the intensity because it would effectively 
flatten out the lines. Thus the fact that our intensities are rather lower than those 
of Flammersfeld may be due to a systematic error. ‘They confirm his results 
to the extent that they are appreciably higher than the values of Ellis (1932) and 
Oppenheimer (1936). 

The reduction in apparent height of a B-ray line due to error in setting the 
field H may be estimated roughly if we represent the line by a gaussian error 
curve 

B= Bye" + Bo, 
where f, is the counting rate due to the continuous background, B,=the peak 
height of the line at H=H,, and «=(H,—H)/H,. Let there now be added an 


476 D.G.E. Martin, H. O. W. Richardson and Yun-Kuei Hsii 


independent gaussian error in H with a precision index hy. ‘The broadened line | 
must have the same area so we get 


B’ = B,(h/hy)e"* + Bo | 

where h-2=hy?+hz?. If we take h,=1-23 x 10? and h,=3-54 x 10?, equivalent i 
to percentage standard deviations of 0-572 and 0:2 respectively, we get a decrease 
in the height f, of the line of 6%. | 
Small drifts in the field strength due to magnetic viscosity also lead to an | 
-underestimate of the height of the lines. | 
A search was made for the lines recorded by Meitner and Philipp (1933) as | 
having moderate intensity. Their line with Hp=541 is below the limit for | 
penetration of the counter window, while that with Hp=1379, presumably the | 
E, line, is too close to be separated from the F-line of thorium B in our case. | 
These authors did not survey the region containing the P-line which was | 
attributed to ThC’’D by Sze (1933). The P y-ray must carry an appreciable | 
fraction of the total y-ray energy. i 
A weak line was found near Hp=1990. This was probably the Jb3 line | 

of Ellis (1932) which he attributed to the internal conversion of the Ga y-ray 
of thorium CC” in the L, level. However, its energy is also close to that expected 
for the conversion of the strong G y-line of ThC’’D in the M, level. ‘The intensity 
we find for this line when using ThC’’D sources supports its assignment to the | 
G y-ray, but it may well be a case of a line with a double origin such as the M line. | 


Surugue (1937) found two close lines at Hp 1978 and 1981 which he attributed |} 


to the conversion of the Ga y-ray in the L, level and the G y-ray in the M, level |f 
respectively. 


$9: THE CONTINUOUS SPECTRUM 
The momentum distribution for an allowed f-transition (Fermi 1934) may | 
be written 
VAP(n F(Z, 9) =ANeo =<) © eee OR | 
where P(y)dy is the probability of emission of a B-ray with momentum lying | 
between 7 and 7+ dy and €, is the maximum value of ¢, the energy of the f-ray |} 
including its rest-energy, expressed as a multiple of mc?. We have used |} 
mc? = 510-89 kev and y= Hp/1704-2. 1 
For a given value of the atomic number Z, 4 is proportional to the square |] 
of the matrix element of the heavy particle transition in the nucleus | 
F(Z, ) = .e*" | T(S +tye/n)P 

where y = 2/137, and S=4/(1—y’). H)) 
F(Z, 7) has been computed with Z=82 for 13 values of Hp and is given in |} 
table 1. It is believed that the values of F are correct to four significant figures. |} 
P(y) is proportional to N, defined in equation (5). The graph of y/(N/F) |} 
against ¢ is given in figure 4 for values of ¢ chosen to avoid the f-ray lines. The |} 
errors to be expected in the values of N are expressed as percentage standard |] 

deviations in table 1. ‘These errors were found by assigning zero error to PD, and I'y3] 
and standard deviations of 1%, 3%, 10% and 5%, to By, By, Bg and T, re- | 
spectively. ‘The uncertainty rises rapidly near the end-point at ¢ =< . | 
The graph deviates from the straight line predicted by (6) for a single 6-tran- |}| 
sition, which suggests that the spectrum is the sum of several partial f-spectra 
ending on different energy levels of the product nucleus Pb2°’, All the stronger 
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partial spectra should be allowed and should have the definite shape predicted 
by the Fermi theory. At low energies the shape predicted may no longer be 
definite because it depends on the coefficients in the combination of terms which 
may be present in the expression for the interaction in B-decay. Evidence on 
this has been given by Lewis and Bohm (1946) for Cu®, but no clear indication 
of the effect appears in the present work. 

It can be seen from figure 4 that the end-point at «=<, is approached almost 
linearly, If we adopt the view discussed in the next section that the experimental 
end-point corresponds to the B-excitation of an energy level 3-20 Mev above the 
ground state of Pb2°8, we cannot expect another strongly excited level until a 


1-0 2-0 3-0 40 
S 


Figure 4. Comparison of the experimental spectrum with the spectra calculated from (a) the 
level scheme of Itoh and Watase and (b) from that of Oppenheimer. 


space equal to the energy of the G y-ray of 277 kev has been bridged. ‘Thus the 
last 277 kev of the spectrum should be linear, but below this there should be 
a small rise in slope. ‘The consequences of such a level scheme are discussed 
-elsewhere, but we may use it provisionally to estimate ey by the method of least 
squares by fitting the four points of highest energy « to a straight line. 

The deviations of the four points were weighted in the normal equations 
inversely as the standard deviations of N given in table 1. The value of (N/F)* 
at the point with « =3-73 mc? was reduced by 4% because it should receive a small 
contribution from the second partial spectrum. We get for the maximum 
energy €) =(4:507 + 0:014)mc?, or (1:792 + :007)Mev. ‘This agrees with the value 
of 1-795 mev found by Henderson (1934) using a method which avoided the 
deposition of aggregates on the source, presumably by carrying out the activation 


an air. 


478 D.G.E. Martin, H. O. W. Richardson and Yun-Kuet Hsti 


The probable error was calculated by means of the standard formula (Margenau 
and Murphy, p. 502). It depends only on the deviations of the four points and 
does not take account of a possible error in our scale of Hp and may thus be an 
underestimate. 


§10. COMPARISON OF ALTERNATIVE ENERGY LEVEL SCHEMES 

If the downward convexity of the spectrum in figure 4 is due to its being the 
sum of a series of m partial B-spectra all having the shape of an allowed transition, 
it should obey the equation 


ue 


VI NIE) =1/ (Age, 2) ee (8) 
n=1 
where A, =0 for «>e,. 

A,, is proportional to the square of the heavy particle matrix element for the 
partial B-transition to the mth energy level of the product nucleus, with end- 
point at e=e,. 

A,, can be connected with p,, the B-transition probability to the mth level, 
which is usually called the B-excitation. 

We have 


P= i pede = Ee (9) 


where p(e)de is the probability of emission of a f-ray in the mth partial spectrum 
with energy between « and «+de. ‘To evaluate the integral we choose as an 
approximation for F(Z, 7) in (6). 
F(Z, 7)=3-8540 + 0-413)5 2) eee (10) 
This gives, on integrating (9) 
Pn= KA, [isen" a 2En = en = i +0-418(€,? ag Lt (ata o ene ae is) 
+ 0-418 dc, In (e, +(€,2—1)83)]. (11) 


m 
K is a constant which is adjusted to make Xp, =1 in a given level scheme. 
1 


Equation (10) differs from the formula of Fermi (1934) in that 0-418 replaces 
0:355. We believe that the value 0-418 covers the energy range which makes 
the main contribution to p,, 1n our case. 

Among the more recent schemes proposed for the energies and f-excitations | 
P, of the levels of thorium D are those of Oppenheimer (1936), Arnoult (1939) | 
and Itoh and Watase (1941). ‘The second scheme differs from the first partly | 
in that it includes a number of weak y-rays some of which are of doubtful origin. | 
It would give much the same shape of spectrum. Using (8) and (11) we have 
calculated the shapes of composite spectrum predicted by the first and third | 


schemes, taking the largest end-point, €,, to be 4-50 mc?. The spectra shown in | 


figure 4 are fitted to the experimental curve at «=2-177 by adjusting K. The. | 
values of the constants are given in table 4. | 
Table 4. The energy level schemes and the resulting complex f-spectra 
Author Oppenheimer Itoh and Watase 
Energy above Bs 
ground level (Mev) { ALES 3-13 3°20 3:13 3°20 3-41 
y-ray transitions | 


to around levels oe es me LX MAX (Gi 
B-ray end point «, 4-50 4-04 3-90 4-50 4-36 3-96 
B-excitation p, 0-08 0-20 0-73 0-45 0-45 0-10 


ia 0-169 O22 3109 0-949 1:16 0-399 
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The discrepancies between the spectra seem to be beyond the limits of experi- 
mental error and are too great to be removed by any changes in the values of Pn 
which would be compatible with the intensities of the y-ray transitions. Upper 
and lower limits to these intensities can be deduced from the theory of internal 
conversion using the observed intensities of the B-ray lines. 

A difficulty in both these schemes is that they do not put the end-point of the 
B-spectrum on the energy level at 3-20 Mev. It was pointed out-by Henderson 
(1934) that in order that the total energy emission may be equal round both sides 
of the ThC branching it is necessary to place the B-ray end-point at 1-792 mev 
on or near the 3-20 mev level of ThD. Otherwise one must postulate different 
rest-masses for the neutrinos emitted in ThC’’D and in ThCC’. In order to 
give a correct balance of the branching energies Arnoult (1939) included in his 
level scheme a partial spectrum with an excitation p, of 10% extending to an 
end-point at 2:10 mev. No £-ray evidence has been found for such a strong 
high-energy tail to the spectrum. 

On adding the f-ray energy of 1-792 Mev to 3-204 Mev, the sum of the M 
and X y-ray energies, we get 4.996 mev for the difference in energy of the 
ground states of ThC” and Pb*°8. his differs by less than the estimated error 
from the value deduced from the «-decay energies Briggs (1936, and others), if 
the B-ray energy of ThC.C’ is taken to be 2:2505 Mev. 

Level schemes which include the P y-ray are discussed by Richardson (1948). 


$11.. THE POSSIBLE s-EXCITATION OF LEVELS 
OU SVEB ROYa iG EINGE RIG Ye 


The graphs of 1/(V/F) predicted by the level schemes can be seen to be nearly 
straight at their low energy ends. 

The experimental curve, in contrast, shows a marked downward convexity 
at energies as low as «= 1-4. 

The curvature, which is very small between «=3-2 and «=2:1 rises sharply 
near the latter point. ‘This excess of slow B-rays may be due to one or more 
partial spectra of low energy. Its magnitude can be estimated by subtracting 
from the observed values of N the values obtained by approaching «=1 along 
. astraight line forming a tangent to the graph at about e=2-2. We get surpluses 
of 16%, 11°%% and 5% in the values of N at the three points of lowest energy 
taken in order of increasing «. ‘These figures are compatible with a single partial 
spectrum with a transition probability p, of about 3% and an end-point «,, of 
about 2:1 mc? exciting a level of thorium D at about 4:45 Mev above the ground 
state. So large a value of p,, suggests that there may be several close levels and 
requires that the excitation is by allowed -transitions. ‘The selection rules 
for an allowed transition give the f-excited level values of spin and parity which 
permit quadripole y-ray transitions to all lower energy levels also excited by 
allowed transitions. As there are very probably three such lower levels we expect 
at least three quadripole y-rays of about 1 Mev starting from a level at about 
4-5 mev. Owing to the low intensities their internal conversion lines should be 
difficult to find, but the weak y-ray of (1350+ 20) kev found by Alichanow and 
Dzelepow (1938) may be one of these rays. Its intensity is stated by Latyshev 
and Kulchitsky (1941) to be 0-036 times that of the X y-ray. 

Analternative explanation of the surplus of slow f-rays is that they are generated 
in the aluminium leaf A by the processes of inelastic back-scattering and ejection 
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of secondary (knock-on) electrons in close collisions. ‘The experiments of 
Flammersfeld (1939), and others, show that these effects do not become appreciable 
at energies of about 100 kev until the thickness of the mounting material is about 
1 mg/cm2. In our case with 0:13 mg/cm? they should be negligible. 

A third possible explanation is that the surplus is due to an increased efficiency 
of the B-counter for slow rays due to their greater specific ionization. It seems 
unlikely that a rise in NV of more than 2% for the two lowest points can be accounted 
for in this way. 

$12. CONCLUSION 

The level schemes considered above all fail both to give agreement with the 
observed shape of the spectrum and to conform to the requirement of energy 
balance in the thorium C branching in a way which is compatible with the experi- 
mental value of the end-point at 1-792 Mev. It seems probable that no energy 
levels of thorium D lying below 3-20 Mev can receive more than 1 or 2% 
B-excitation. In order to account for this low excitation it is evident that the 
lower levels can only be accessible via forbidden f-transitions and that the transition 
to the ground level must be highly forbidden. 

Owing to the presence of five strong y-rays and at least three weak ones a 
great many level schemes are possible. A clue, which may, however, bea chance 
coincidence, is provided by the P, G, and M y-rays, which satisfy the combination 


(Av) q+ (Av) = (Av)p. 


The evidence is given in table 5. 


Table 5 
Author * Ellis (1932) Surugue (1937) 
Energy of G y-ray (kev.) 276:7 276°8 
Energy of M y-ray 582°5 583-4 
Sum 859-2 860-2 
P £-line 771-6 WHS 
Surplus 87-6 87-9 


The surplus agrees fairly well with the K-ionization energy of 88-10 kev | 
for Z=82. Each author’s values for the energies of the f-lines have been used, 
but the ionization energies are those of Ingelstam tabulated by Arnoult (1939). 

If we assume the spin of the ground state of ThC” to be four units of h/2z, 
the Gamow-T eller £-selection rules for an allowed transition permit the excitation 
of levels of ThD with spin quantum number J=3, 4 and 5. The distribution 
of convexity in the graph of 1/(N/F) shows that levels with these spins are present 
near the high-energy end-points between 3-1 Mev and 3-75 mev and that after 
an interval they appear again at higher energies above 4:4 Mev possibly in the 
form of a close group of levels. The ground state of ThD has very probably 
I=0 and the spacing of the levels would be expected to get smaller at high energies 
of excitation. 
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ABSTRACT. A method is presented for solving scattering and eigenvalue problems. 
This is of use mainly in cases where the forces are of short range. ‘The method is based 
on a Taylor expansion of the radial eigenfunctions about a point where the potential is 
already rather small. A procedure is given for determining the logarithmic derivative at 
this point of the solution which is “least singular” at the origin, this is in general the 
physically admissible solution. 


§1. INTRODUCTION 

tT is well known that for calculating scattering cross-sections one does not 
require the detailed solution of the relevant wave equation, but only its 
asymptotic form. ‘Thus for a system with forces of finite range a, the 
scattering problem is effectively solved when one has determined the logarithmic 
derivatives of the interior radial solutions at the boundary a; similar considera- 
tions apply also to the corresponding eigenvalue problem for the bound states. 
In more general cases where the forces are of short range one may solve the 
scattering or eigenvalue problem by subdividing the whole range (0), «) into an 
“interior region”? (0,@a) and an “exterior region” (a, 00). In the exterior 
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region, where the potential is already rather small, the solution may be obtained 
by a perturbation method together with the condition that this exterior solution 
fit smoothly on to the interior one. Here again the only information about 
the interior solution needed in the calculation is the value of its logarithmic 
derivative at the point of subdivision a. 

The object of this paper is to present a method of evaluating the logarithmic 
derivative of the interior solution at a, without solving the equation explicitly. 
Usually the origin is a singular point of the radial differential equation. The 
customary procedure is then to expand the solution in a power series about the 
origin so as to satisfy a condition of regularity or finiteness there. ‘The present 
method too is based on a Taylor expansion, not about the origin however, but 
about the boundary or subdivision point a. The condition of regularity at the 
origin is then equivalent to a condition on the radius of convergence of the series. 
It is this latter condition which permits the direct evaluation of the logarithmic 
derivative of the interior solution at the boundary. 


§ 2 [METHOD SORE CAW CUO AMBION 
In general the radial differential equation is linear, of second order and with 
the origin as a regular singular point : 


au du ; 
ae +z + p(rjaaOen > 1 5°) ae (24) 


The equation may also have singular points other thanr=0. We are concerned 
with that solution of (2.1) which is regular at r=0; more particularly we wish 
to determine its logarithmic derivative at the ordinary point r=a. 

Let r=0 be the singular point of (2.1) nearest.to a; and let d be the distance 
from r=a to the next nearest singularity. Since a is an ordinary point of (2.1), 
any solution can be expanded in a Taylor series about a : 


u(r)= © ¢,(r —a)" 
n=0 


n\\ dr® 


In general the radius of convergence of this series will be just a, and so by Cauchy’s 
convergence test : 


1 1 (Fs) 
where C¢, = —u,=— «5% ‘egy loetaemelete eee (22) 


lim 
N>D 


Cy |Pe LG. on ak a 2 Se ee (2.3) 


But when, and only when, u(r) is that special solution which is regular at r=0, 
the radius of convergence will exceed a and will in fact be d; i.e. 


lim |G, (Ug, or (2.4) 


It is immediately evident that c, can be represented in the form : 

on = F 60 a Gila, : 
where Fy=G,=1, F,=G,=0 and the coefficients F,, G,, (n=2, 3. . .) are 
uniquely determined by the differential equation. Let b=w'(a)/u(a) denote 
the logarithmic derivative at r=a of the regular solution. We then have 


b= —lim F,/G,,. eM BERS 


n> 
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This can be proved as follows: We are free to make either cy=0 or c,=0 in 
(2.5). ‘Then atleast one of the two corresponding functions u(r) must be singular 
at the origin and thus lead to (2.3). Hence at least one of the relations 
fina ei ela lim |G, "= 1a 
n>D n> D 
is valid. For definiteness suppose the second equation holds. Then when 
ic, /Co—0 we have : 
1 /d=lim | F,,¢9+ Gc |” 
n> 
=lim {|¢9)!"| GP" |(Fy/G.) +52} 
> 


and so it Gwe ol aldees | A a added (237) 
which establishes (2.6). 


' Equation (2.7) indicates that the convergence to the limit in (2.6), once it 
“sets in’’, is rapid provided a/d is not too close to unity. It also provides a 
useful control in some applications, enabling us to decide whether or not an 
empirically determined convergence is real. 

In general it is not feasible to attempt to find the exact limiting value b in 
(2.6). ‘The main use of (2.6), (2.7) in practice is as the basis for an approximation 
method; and it is just this feature which we wish to emphasize in the present 
paper. The procedure then consists in regarding b, = —F’,,/G,, for some value 
of m as representing 5, the particular choice of m being guided by the convergence 
criterion and the order of accuracy required. 

To illustrate the procedure we shall now discuss a number of examples : 
In these we shall write (cf. (2.2), (2.5)). 

SEL Ahab Selb Akad ey A Tinal-4 Oh Lh. Ani Pee tee (2.8) 
where A,, is some convenient numerical factor. ‘Then 
Wee = Ay, it I OF = Ay G,,/n 
and the mth approximation to d is given by 
Dee One Be)! Se aaa (2.9) 


§3. EXAMPLE 1 


a’ 
Consider the equation r = ST arte TTL DY AUS AER Se (3.1) 


which has as solution regular at y=0 the function y/7r1,(2/7), 1, being the Bessel 
function of imaginary argument. Differentiating (3.1) (7 —2) times and sub- 
s‘ituting r=a we obtain a recurrence formula for the w,, of (2.2) : 


au, = —(n —2)Uy_y + Un_o- se eeee (3.2) 
Writing au, = Pyiy + Oy 
it follows that the P, and the Q;, separately each satisfy recurrence relations 


of the same form (3.2) : 
ee —(n Se aes } 


O,.= —(n —2)Qn-1 + 4Qn-23 i sevens (3-3) 
Po=1;° Qo=0=P; Qi=4. 
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Repeated application of (3.3) then leads to : 


Vis: QO, =0 

P= =a Oa 

P,=2a+a@ Oj 2a 

P,= —(6a+ 4a’) O,=0e a 

P,=24a+ 18a? +a? O,= —(24a? + 6a?) 

P,= —(120a + 96a? + 9a?) QO, = 120a? + 36a? + a* 
P,=720a + 600a? + 72a? + a* O,= —(720a? + 240a* + 12a*) 


P,= —(5040a + 4320a? + 600a? + 16a*) QO, =5040a? + 18004? + 120a* + a°. 
The nth approximation, 6,, to u’(a)/u(a) is given by (2.9). To indicate the | 
order of accuracy attainable and the mode of convergence, successive values _ 
of b,, for three different values of a are listed in table 1. (The entries for. 7-703 
were obtained by use of tables of Bessel functions). 


Table 1 
n i (Gs) by (a=4) by (a@=9) 
3 1-000000 0-250000 0-111111 
4 1-500000 0:750000 0-611111 
5 1:-428571 0-550000 0°311111 
6 1-433333 0:583333 0-380342 
7 1:433121 0:578571 0-362540 
8 1:433127 0:579060 0:365300 
9 1:-433127 0:579022 0-365300 
oa) 1:433127 0:579024 0°365353 


Even for comparatively large values of the argument accurate values of 6 can be 
obtained rapidly. ‘To obtainscomparable accuracy with the use of a Taylor 
expansion about the origin requires the retention of a large number of terms and _ |f 
more laborious computations. | 

It is interesting to remark that the procedure may be extended to obtain 
representation of the function u(a) = 1/al,(2,/a) over a fairly wide range in terms | 
of elementary functions. Regarding b,= —P,/Q, as representing u’(a)/u(a) | 


with sufficient accuracy in a range (0, A), an elementary integration leads to the } 
result : ! 


ula) = aeat? (= + 20a zs =" | ia —24/10)a + 60 7125/12, V10, 
60 (10+ 2,4/10)a +60 
Here the constant of integration has been fixed by using the known behaviour 
of u(a) in the neighbourhood of the origin. The accuracy of the formula (3.4) 
decreases gradually as a increases. At a=9 it is in error by less than 0-01 %. 


| 
| 
| 
| 
§4. EXAMPLE 2 | 
| 


2 
Consider the equation r(r+1) ce U0 Sg at oe eee (4.1) 
r 


and suppose that the logarithmic derivative of the regular solution is required | 
for r=2. Here the direct method of series expansion about the origin breaks 
down on account of the singularity at y= —1; the series would have radius of | 
convergence unity. (Of course an expansion about the origin can still be used in | 
conjunction with a bilinear substitution, e.g. one which transforms —1 into w). 
It is here that the method of this paper can be used with particular advantage. 
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Writing (Gera) Upc) l= we Where (4.2) 
we obtain the recurrence formula : 
P,,= —(n —2)(2a+1)P,,_ , —[(n —2)(n —3) -1](@2 +4)P,,_» 
Qn= —(u —2)(2a +1), —[(n —2)(n —3) -1](@? +.4)Q,,-2 J 
avith P51, Oo al QO1,=(a +a). 
peccreding as in §3, successive approximations to m/% (with a=2) are given 
by bn = —PalQn : 
b,=0-8333 Dy = 07333 b; =0-7083 b, =0-6988 
b, =0-6945 b, = 0-6923 Dye= 0-69 11> 350-6903. 

é Here the differences between successive terms b,,, (k>4), agree quite well 
with the convergence criterion |b —b, |~K(2/3)" with K depending only slightly 
on. ‘This fact may be used to obtain a more accurate value for b : 

Oe = b,, —K(2/3)" = b,, =2(b 4 on): 
The corresponding values of b’, are b',=0-6879; b’,=0-6887; 6',)=0-6887. 


§5. EXAMPLE 3 
The differential equations 


ae at [e- i +B |u=0 ie) 


(J=0, 1, 2,...) occur in the theory of scattering with a central “‘ meson potential’. 
We consider the case /=0 and write 


CN Se RE ed ee a (5.2) 


V(x) must then be finite at x=0. By repeated differentiation of the differential 
equation for V(x) we obtain for V,=V (a) the recurrence formula : 


QV nso t+(Mt+2)Vniyth[aV, +nV, |= —Be“ Py (- 1y6(%) ers Garces (a3) 


Writing a”V,,=P,,Vo+QnV1 with Py=1, Qp=P,=0, OQ, =4, we obtain recur- 
rence formule Pe P., O,,separately : 


i +2 — —(n Te) Lappe Rai, rate) —aBe™ & ( or L)e @ tea 0 

k=0 
_ with an identical formula for the Q,,. Taking the case k=0-8, B=1-5 and 
1=0, for which the asymptotic phase has been obtained by numerical ee 
(Sachs and Goeppert-Mayer 1938), the following values of 6, = —P,,/Q,, 
obtained with a=2: 


b;= —1-161 bg = —1:284 [NOSE 8S) b= —1:375 

y= 1369, 9 0,5= 1396 0, = = 1-399: 
The differences * between successive b,, (n>4), behave regularly and provide 
-an improved extrapolated value for b=V'(2)/V(2): b= —1-401, and so 
u'(2)/u(2) = —0-901. 


* In most cases it is possible to obtain an improved value in this way by extrapolation using the 
differences between successive values of by. 


2 
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curate value for the phase shift 7 of the exact wave function 


Phe asymptotic phase », satisfies the exact equation * 
~X 
— Jae : z 
ksin(y)—7)=B| —wu(x)sin(kx+n)dx. ...... (5.4) 
; x 
Since the potential is already rather small for x>2, a fairly close approximation} 
nO yp Can be obtained by evaluating the integral in (5.4) with sin (kx +) substituted} 


Sor wx). A straightforward numerical evaluation of the integral then leads to} 
sin (y,—y)=00217 and so »,=0-8371. This is in good agreement with the} 
valne »,= (083706 obtained by Sachs and Goeppert-Mayer (1938). 


cedure converges more rapidly than for /=0; 
rapid is the convergence. This feature is not 


It can easily be understood in terms of thef 
relative degree of singularity of the regular and irregular solutions which increases} 


F (xyh{ Rx) +2'(x)x (kx) =0 
we obtain 2 pair of first order equations for jf, g < 


. a) 


j(x)=- E v(x) u(x) x, (Rx), ] 
> 


(7p. Pe? ae (5.5) 

ge (x)=+ % O(x) w(x) a, (Rx). ; | 
The condinen 7”(@)=0 then leads to an approximation for g(a)/f(a) and so td 
2h approximation for the phase shift at x=a. The repeated differentiationg 


are more Complicated than in the above direct method, but f, g are moray 
mg functions than w and the procedure would thus converge more 


$6. APPLICATION TO TENSOR FORCES 


ters the pair of simultaneous equations (cf. Rarita and Schwinger 1941), 
Y Sis? =— : = az ’ = 6 
Ze tt a@)eee)e; [St & ae +2(0) | wW=WU,3(x)u ......(6.1) 


. Tp, U, are Specified potential functions. The physical solutions are required 
‘Sty, amongst other conditions, that of vanishing at least like x at the origin 

or sumpheoity we shall consider the case of a rectangular potential; v,, v, 2 | 

én constant m (0, 1) and vanish for x>1. ‘ | 


5.4) as easily obtained e.g. by a slight modification of the argument leading to equation (6.6} 
m Bethe end Bacher (1936). 
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An examination of (6.1) in the neighbourhood of the origin shows that the 
physically admissible solutions are of the form : 


n=(0 n=0 


D v0) 
u(x)= 2 Aw *t + In x Cx t, i 


D D 
w(x)= 2 Boats in « Dx, | 
n=0 n=0 


‘The constants Ay, By can be assigned arbitrary values independently. There 


are thus exactly two independent solutions which vanish as required at x=0. 


On account of the linearity of (6.1) any linear combination of solutions is again 
asolution. A fundamental system of solutions 


= OAR) u= (x) : 
w= a ay 


can then be chosen so as to satisfy two arbitrary independent pairs of conditions 
atx=1. In particular, there exists one such system with the property 


aaa Sena (6.4) 


BU) =0 J p(1)=1 

Once the derivatives «’(1), B’(1), 6’(1) and (1) are known for the system (6.3) 
satisfying (6.4) (or any other set of conditions), the scattering problem is, in 
effect, solved. It only remains to calculate the asymptotic phases by fitting 
appropriate exterior solutions (expressible in terms of Bessel functions) on to 
the interior solutions at x=1. A knowledge of the above first derivatives of 
the interior solution is sufficient for this purpose. 

.The discussion in the preceding sections has been based on the translation 
of a regularity condition at the origin into a condition on the radius of convergence 


of a power series. ‘This point of view is, however, unnecessarily restricted. The 


selection of the regular solution is a special case of a more general property of the 
procedure. It is in fact applicable to problems in which the required solution 
is not regular at the origin. In this case the method selects the “least singular”’ 


solution consistent with the specified data*. A proof of this, in a form suff- 


ciently general for our purposes, is given in the Appendix. 
Now, the physical solutions of type (6.2), although not regular at the origin, 
have the property of being “‘least singular”’ in the sense of the Appendix. The 


_ remaining solutions become infinite at «=0. Our method can thus be extended 
_ to apply to this problem. 


We introduce new dependent variables U, W: 
Z=nU; WaweWes een a” OO eet (6.5) 


Then repeated differentiation of the differential equations for U, W, leads to 


recurrence relations for U,,= U1) and W,=W (1). ee 
X= UL, 41 As, | 
Y,,= W,,+2nW,_1+n(n —1)W,_2, 
aie MIDE. fez on (6.6) 
Zy,=(RP+U)X,, =U5y 5; | 
T= (R240) Vn —U34 n> J 


* We are indebted to Professor Peierls for this suggestion. 
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we have Uj is= —(1+2)0 7-2. \ 
Wris= —(2n+4)Wyiy —(n -1)(n+4)W,, —Tr. J 


(The relations, =a, 
of (6.6), (6.7)). 


We now write 


Y,.2=6W,,—T, are useful in application 


On = LP U, +0, OF iy Wot Sm W;; I. 
Wy=K,,U67 Lo Ie J 
Then 

( P)=0,=M,=N,=1. 


All other P,,,...; N,, with m=0, I zero. } 


It is evident that the pair of coefficients (P,,,, K,,), (w=0,1,2...) by themselves. | 
satisfy the recurrence relations (6.6), (6.7) with P,, replacing U,, and K,, | 
replacing W,,; the same is true for the remaining pairs (Qn, Ln), (Riny Min), | 
(S',N,)- This simplifies the calculation of successive U,,, Wn. 
The procedure is now to set both U,,=0 and W,,=0 in (6.8). Putting | 
U,=1, W,=0 then leads to the equations : 

0,0, +5SnMi=—Pas 14077N = oe (6.10) 

The values of U,, W, thus obtained provide the required first derivatives 
eat Oe B= | 
of one solution of the fundamental system (6.3, 6.4). 4’(1), &’(1) are obtained | 
in the same way by solution of : 
0,04 ae SinrW, ==. 
together with g(1)=U,; ¥'(1)=W,42. 
To illustrate the mode of convergence of the procedure, successive approxi- | 
mations to («’(1), 6’(1); (1), #/(1)) are listed in table 2. As our object in the 
present paper is exposition of method rather than actual calculation of cross- | 
sections the calculations for table 2 were performed with the simple values 


R+v,=3, k?+v,=—1, v3= —6. 


m > Ly,01+ NW, = =I, 


m 


Table 2 
n a’ (1) B’ (1) ¢’ (1) ob’ (1) 
3 SS PAGS —3-2308 +0-4615 +3-2308 
4 +-0-:5000 —1-3500 —2-8750 +2°5625 
5 —0-3421 —1:3953 —1-1245 +3:1031 
6 — 0-8654 —1-5320 —1-7084 + 2-8668 
W —1-0189 —1-°5711 —1-5093 +2-9178 
8 —(0:8721 —1-5316 —1-5381 +2-:9109 
9 —()-8634 —1-5287 —1-5260 +2-9141 
10 —0:8694 —1-5300 —1:5305 +2:9131 
14 —0-8690 —1-5297 —1-5292 +2-9134 
12 —0-8687 —1:5293 —1-5291 +2-9135 
APPENDIX 


Let f(z) and g(z) be analytic functions having no singular point other than |} 
z=0 in some circle |x —a|=6, (b>|a]). Without loss of generality, a may be |} 
assumed real and positive. If either f or g or both are many-valued with the 
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origin as a branch point, a cut is made along the negative real axis from 0 to 0; 
we then restrict our considerations to a single branch of fand g. 


Suppose, further, that in the cut plane: (a) there exists a real number o, 
(0 <o <1), such that 


Sime (Z)- Oras ere (9 Ne Pe Te, (7.1) 
(b) g(z) can be expressed as : 
£2) =e@) Fue) el ee 8 kis: (7e2) 


where 4(z) satisfies the same conditions as f(z); the “more singular” part :b(z) 
will be assumed to have either of the forms : 


eA; 
Y2= 2 Sey ae oe (1) 
C= 1.8 
ey, 
(2) =In ue ik pieces (II) 
ByNELGe Ay ee 
We consider the Taylor series 
[ sAx)= E flea, 
; : 5 
Ls = cers g(x Mac a)" = oy (bn + Pax ==)" 
and shall prove that ie nl ig, = 0, 


Let h<b—a and let @; ee the 
contour P+ Y;+ ys + L;5 in the cut plane. 
T is the circle with centre a and radius 
a+h, ys; the circle [lor L,and Lsare 
straight-line paths between —A and —6 
along the upper and lower lips of the cut 
respectively (see figure), ‘Then 


mei A(z) 
faa (z— Gane 


ail eae meas rane 


> Let ae = maximum of | f(z)| on I, 


m(h) = maximum of | z'~°f(z) |on the 
two lines —h<z<0, | 


M,(8) = maximum of | f(z) | on y;. 


‘Then 
1 ede ale oe gi f(z)dz 
Pale (2 —q)rti — L; 2-2 Sart 
L mh) (  \\rde 
S25 (a+5yr eee 
h°m(h) 


Qro(a+d)rt 
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Similarly for the integral along %;. Hence 
M,(h) 2h°m(h) 5M,(8) 
nl < (a+h)”" 2 2mo(a+8)"+1 — (a —8)h4}- 
By (7.1), 8M,(5)+0 as 80. Thus, proceeding to the limit 60 : 


lfrl< {Mh (1 4s a as ny 


oma 
Now, given «>0, however small, we can choose h so that 
h°mh) _ « 
70a De 


Then we can find a number Ng such that for all n >No 
M,(h)(1 +h/a)-" <e/2. 
Hence, for all n>Np, 


lf.|<elate 8 ee (7.3) || 
On the other hand, by direct differentiation of %b(z) in the forms (1), (II) we 


have at 
ne { (1) —=- (ae cae )+ a, ee (I) 


| 
Yu= ae Cede pee ae Ina— ae ae 
Qirn n AE k nr 
Lie 


d 
Hie hie ge een (11) 
where the dots denote terms which are of lower order as noo. In both caseS 
the coefficient of a+” is at least non-zero : in general it becomes infinite as |} 
n-> oo (except for A,=1 in case (1)). Thus, for both (1) and (II), there exists | 
a number N, such that, for all n>N, | 
en |>M,/a"+",. 
where M, is a constant independent of n. | 
Thus, combining (7.3), (7.4), A follows that for all n=>Max(No, N;), i 
lfrlbn | <a’e/Ms, and so bo n fn a0! | 


n-> 


Since ¢,,/¢,,>0 as n— co ree behaves ‘‘at worst”’ like f), we have lim fniSa= Mg 


If now f(x), g(x) form a fundamental system of solutions for a Sen second |f 
order differential equation, the general solution is of the form | 


as U,(x —a)" = Af(x) + Bg(x). - Then u, = Af, + Bg,. 
n=0 


Our procedure is to set u,=0. This gives B/A= —f,/g, which -0 as| 
moo. In this w ay the ‘least sinenler? solution f(x) is selected. | 
The proof can easily be extended to the case of a differential equation of | 
order k, with (Rk =) conditions specified at x=a; the method will again select |] 
the “‘least singular’’ solution consistent with the Heed data. 


| 
i 
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Characteristics of Radio Echoes from Meteor Trails: 
I. The Intensity of the Radio Reflections and Electron 
Density in the Trails 


ByTA. ©C.sBy EOVELL AnD Jo Ay CLEGG 


Physical Laboratories, University of Manchester 
Communicated by P. M. S. Blackett ; MS. received 26 Fanuary 1948 


ABSTRACT. Recent observations of radio echoes reflected from meteor trails have given 
a large amount of information about the electron densities in the trails and the physical 
processes occurring in the atmosphere after the passage of the meteor. In this paper, which 
is Part I of three papers dealing with the work, formulae are derived for the intensity of the 
radio wave scattered from a meteor trail on the assumption that the electrons are created in a 
long narrow column, of diameter small compared with the wavelength of the radio wave. 
Experimental work is described which shows that the predicted variation of received power 
with wavelength is correct for the wavelength range 4:2 m. to 8-3 m., and according to 
preliminary results down to 1-4m. The formulae can then be used to measure the electron 
density in the trails of meteors. Of particular interest are the measurements for meteors 
which are also observed visually. The results show that the density in the trail of a 5th 
magnitude meteor (on the limit of naked eye visibility) is approximately 2 x 101° electrons 
per cm. path. Brighter meteors (~ magnitude+1) produce 10 electrons per cm. path. 
These results are in very good agreement with contemporary theoretical calculations. The 
theory is compared with earlier calculations made by Eckersley and Pierce. 


si INTRODUCTION 

HEN a meteor enters the earth’s atmosphere it creates a trail of electrons 

at a height of about 100km. above the earth’s surface. ‘The trail is 

dense enough to reflect energy from an incident radio wave and under 
suitable conditions a short-lived response can be obtained on the receiver. ‘These 
transient radio responses have been studied since 1929 (see Lovell (1948) for a 
review of this work) and recently the technique has been developed for the study 
of the activity of meteor showers and the delineation of meteor radiants. 
(Appleton and Naismith 1947, Hey and Stewart 1946, 1947, Prentice, Lovell 
and Banwell 1947, Lovell, Banwell and Clegg 1947). ‘The present papers are 
concerned with the mechanism of reflection of the radio waves and the physical 
processes occurring in the atmosphere after the creation of the meteor trail. 
In Part I formulae are derived for the intensity of a radio echo received from a 
column of electrons in the atmosphere; from these can be calculated the density 
of electrons in meteor trails. Part II will deal with the distribution of electron 
densities in the trails of the shower meteors and with the distribution of meteoric 
mass. In Part III the factors which determine the duration of the transient 
echoes and the behaviour of the trails after formation will be discussed. 

S2PCALCULATION OFTHE INTENSITY OF THE RADIO 
REFLECTION FROM AN ELECTRON CLUSTER 


Blackett and Lovell (1941) have calculated the intensity of a radio wave 
scattered from a long narrow column of ionization in the atmosphere, assuming 
that the electrons scatter coherently. These calculations will now be extended 
to include the parameters of the radio apparatus and will be applied to the parti- 


cular case of meteor trails. 
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Consider a cluster of N electrons at a distance Rcm. from a radio transmitter 
and receiver working on a wavelength of \cm. If the size of the cluster is small 
compared with ) and if the electrons reach their full velocity under the influence 
of the impressed F.M.F. and are not impeded by collisions, then the cross section 
for scattering by the N electrons will be (87/3)(e2/mc?)2N*cm?, where e is the 
charge, m the mass of the electron and ¢ the velocity of light. 

If the peak power in the pulse of the radio transmitter is P watts, and if the 
aerial system of the transmitter has a power gain G’ over a source radiating 
uniformly through 47, then the power density at a distance R cm. will be 
PG'/4r7R? watts/cm?. Hence the amount of energy scattered by the electron | 


2 \2 PO’ 
SistCnnial ee (=) pee a i. 


3) \mc 47R? 

If the receiver is situated in the same place as the transmitter, the power 

density of this scattered radiation when it returns to the receiving aerial will be 
2 \2 ae 

Sn (<5) mie ae watts/cm?. (The factor 1-5 is introduced to allow for 

the fact that the electrons scatter as Hertzian dipoles and not uniformly 

through 47). 

The effective collecting area of the receiving aerial will be G)'A?/4a7cm? 
where G;’ is its power gain over a spherical source. Hence the power which it 
delivers to a matched load when situated in a field of mean power density o will 
be G)'A?o/87. Thus in the case under consideration the amount of energy 


delivered to the receiver input after scattering from the cluster of N electrons will 

be 

37 watts: ., < Deceee ee (1) 
For comparison with experiment it is convenient to take the gain of the 

transmitting and receiving aerials as the same, and to express the gain in terms 

of that of a half-wave dipole so that G’=G,'=1-64G, where G is the power 

gain of the aerial system over a half-wave dipole. ‘Then (1) becomes 


<= 37 me? 1672 Rt oe 


8 ( e s 7 OL G Gane 


2 \2 2\2 
= on (=) N? FO watts, ook. (2) 
In terms of the voltage amplitude V at the receiver input 
V=S2re, Se eee (3) 
where ¢ is the input resistance of the receiver. Or for constant P 
Vek XAN/R?, yee gee (4) 


2 PG? 


2078 aun 


2 
where the constant & is given by ke on (=) 


§3. APPLICATION TO THE CASE OF A METEOR TRAIL 

If it is assumed that the meteor trail consists of a long narrow column of 
ionization whose diameter is small compared with A, the appropriate number 
of electrons N to be included in (2) can be calculated by using ordinary optical 
diffraction theory. The voltage amplitude returned to the receiver from a 
small element of track at an angle 6 to the perpendicular from the receiver to the 
trail is given by dV =kiad6/R, where « is the number of electrons/cm. path in 
the trail and & is the constant occurring in equation (4) above. 
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The voltage amplitude from the portion of the trail @=0 to 6=0, is then 

‘given by Fresnel’s integral in the form 
3/2 v 3/2 - 
Y= alm) fs cos V8dV + 7 (R) [sing veav, 

where V = 264/(R/A). 

The voltage amplitude at the receiver due to the whole track is then 
V =(1/42)ka(A/R)??. Hence by comparison with (4) the appropriate number 
of electrons N is given by 

INEHOA URE) re ay oa (5) 

In this calculation the phase correction is applied twice because of the curvature 
-of the incident and reflected wave surfaces. The wavelength is then effectively 
reduced to A/2, giving equation (5) instead of the more familiar form N = a /(AR) 
used by Blackett and Lovell (1941). 

Hence from (2) the amount of energy at the receiver input due to scattering 
from a meteor trail with an electron density of « electrons/cm., will be 


= 8 e? \2 i AR PG?) ae PG?? / e2 \2 os 6 
e= 3 7 Fae a 7 32,3 R1) = a 472 R38 \me AWANEES, coo nce (6) 
; me?\ | eR? 
Alternativ ely “L= a/ 247 (=) G JF). Sabch Sak en (7) 


S45 WMSUS WIP IE(CME Ue IIRC AN KON) (SCID IL USINOUNKS 
The above calculations assume a cross section for scattering of 
(877/3)( e?/me?)2N2, 
This is correct only if the electrons are free, and it is necessary to examine the 
assumption, since an electron may collide with a neutral atom or molecule before 
it has attained its full velocity under the inflyence of the impressed E.M.F. This 
-effect was neglected in the calculations of Blackett and Lovell (1941) and it will 
be considered here for the particular case of meteor trails. 

Let v, be the collisional frequency of electrons with neutral atoms and vy, 
the radio wave frequency. ‘Then the cross section for scattering is effectively 
reduced to 

(Sas)\(etme*?)t NA (vijay) cee (8) 

Except at high altitudes this term has a considerable influence in any prac- 
ticable experiment. For example at 20 km. altitude in the atmosphere an electron 
will make approximately 10 collisions/sec. with atmospheric molecules. Now, 
-even with a wavelength as low as 10cm. (v)=3 x 109/sec.), (v,/7v)?—= 100; and 
for wavelengths in the region of 4m. it rises to about 2:5 x 10%. 

In the particular case of meteors the radio echoes are observed to originate 
_ataheight of approximately 100km. At this altitude v, 3 x 10°/sec. and except 
for long wavelengths (e.g. A>100m.) the effect of damping can be neglected. 
The assumption made in deducing equation (6), that the electrons are free, 
is therefore justified. 


§5. EXPERIMENTAL TEST OF THE THEORY 
An attempt to verify equation (6) experimentally as regards the predicted 
variation of « with ’ has been made by observing the same meteor trail simul- 
‘taneously on different wavelengths. The essential features of the apparatus 


494 A. C. B. Lovell and 7. A. Clegg 


have been described elsewhere (Prentice, Lovell and Banwell 1947; Lovell,. 
Banwell and Clegg 1947). The wavelengths used were 8-3 m., 6-5 m., 4-2 m., and 
1:4m., and the main characteristics of the apparatus are summarized in table 1. 
A preliminary account of the measurements on 8-3m., 6-5m. and 4:2m. has 
been given by Lovell (1947). 


Table 1 
A (m.) (a) (b) (c) (d) 
8:3 150 15 50 4-2 
6:5 50 NS 50 4-6 
4.2 50 8 50 De 
1-4 50 2 180 oe 


(a) peak transmitter power P (kw.); () pulse width (u sec.) ; (c) pulse recurrence frequency ;. 
(d) receiver sensitivity : voltage input for twice signal/noise (V.). 


For each equipment the aerial systems consisted of half-wave dipoles mounted 
at equivalent heights above the ground such that in equation (6) G can be treated 
as constant. Observations of the transient echoes were made visually on a 
cathode-ray tube display. The quantity measured is the amplitude of the echo 
in terms of the normal receiver noise level. This is converted into equivalent 
voltage input at the aerial by appropriate calibration of the receiver. Thus if 
simultaneous echoes are observed from the same meteor trail on wavelengths 
Ay, Ap of amplitude V,, V, equation (6) predicts that V,/V_=(«/€2)* =(Ay/As)*® 
since R, « are constants for the same trail, G is assumed to be the same and the 
equivalent voltage inputs are corrected for any variations in P by assuming 
V,/V2=(Py/Pe)t. 

In the actual reduction of observations the exponent m in the expression 


V;/ V, = (A,/A2)” 


has been calculated for each pair of simultaneous observations for comparison 
with the theoretical value of n=1°5. 


(1) Stmultaneous Observations on 72 Mc/s. (4:2 m.) and 36 Mc/s. (8-3 m.)* 


Simultaneous observations of the radio echoes from the same meteor trail. 
on wavelengths of 42m. and 8-3m. using apparatus with the characteristics 
specified in table 1 have been made in 1947 as follows: 12-27 August, 104 
observations on the Perseid radiant (21 observations rejected) ; 23-24 September, 


158 observations on an unknown radiant (17 rejected); 9-10 Octeber, 59 observa-. 
tions on an unknown radiant (2 rejected). 


For each of the 321 pairs of observations the exponent 1 has been calculated 
from equation (9) as described above. The distribution curve for these measure- 


ments is plotted in figure 1 and gives a mean value of 1-4 for nm with a standard 
deviation of 0-5. 


* Note added in proof. Equation (6) applies to the case of scattering from a meteor trail which 
crosses through the aerial beam in such a way that the foot of the perpendicular from the station 
to the trail lies within the beam. This is the only condition under which the main radio echo can 
be observed on wavelengths shorter than6'5m. On longer wavelengths some experimental evidence- 
exists that the trails lose this critical aspect sensitivity (Millman, McKinley and Burland, Nature, 
Lond., 1948, 161, 278, Lovell 1948). The measurements on 8:3 m. in this paper apply to the case 
of perpendicular reflection, since the only echoes selected are those appearing simultaneously on: 
the shorter wavelength, which implies that the aspect condition is satisfied. 
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Figure 1. Distribution curve for simultaneous observations on A=4:2 and 8:3 m. 
Mean value n=1-4; standard deviation=0°5. 


(11) Stmultaneous Observations on 72 Mc/s. (4:2 m.) and 46 Mc/s. (6-5 m.) 


A similar series of simultaneous observations on wavelengths of 4-2 m. and 
6:5m. have been made as follows : 12-27 August 1947, 57 observations on the 
Perseid radiant (11 rejected); 14 December 1947, 79 on the Geminid radiant 
(17 rejected); 3-4 January 1948, 30 on the Quadrantid radiant (2 rejected). 

The distribution in values of the exponent n calculated for each pair from 
equation (9) is plotted in figure 2 and gives a mean value of 1-7 for m with a standard 
deviation of 0-6. 
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Figure 2. Distribution curve for simultaneous observations on A=4:2 and 6:5 M. 
Mean value n=1-7; standard deviation=0:6. 


(iii) Triple Simultaneous Observations on 208 Mc/s. (1-4 m.) 72 Mc/s. (4-2 m.) 
and 46 Mc/s. (6:5 m.). 


Preliminary attempts have been made to observe the same meteor trail 
simultaneously on wavelengths of 1:4, 4:2 and 65m. Equations (6) and (7) 
indicate that the voltage amplitude at the receiver can be expected to decrease 
as.A3/2and hence only the.trails with very dense ionization can be expected to give 
echoes on the 1-4m. equipment specified in table 1. ‘This expectation has been 
realized and only four triply observed radio echoes have so far been measured. 
These yield 12 values for the exponent 7, calculated from equation (9) with 
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mean values as follows : 4 4:2m: and 6:5m., mean n=1-7; A 65m. and 1:4m., 
mean n=1-5; 1 4:2m. and 1-4m., mean n=1-4. 


§6. DISCUSSION OF RESULTS 


‘There are two principal uncertainties in these observations. (i) Although 
geometrically similar aerial systems were used in the comparison equipments, 
slight differences of terrain and surroundings inevitably introduce variation in 
the polar diagrams of the aerials. It can therefore be anticipated that a certain 
number of meteor trails will form in regions where the polar diagrams of the com- 
parison equipments are not exactly similar. (ii) The measurements of signal/ 
noise ratio were made by visual observation of the cathode-ray tube display. 
Two observers may not agree exactly over the amplitude of a single echo, and 
for a given observer there will be a certain inconstancy in his estimates of amplitude. 
Any systematic bias due to the first effect has been counterbalanced by alternating 
the observers between the equipments but the second effect could not be eliminated 
in this work. 

Both (i) and (ii) will give rise to a considerable scatter of the results about | 
the mean value. In the results described in §5(a) 40 of the 321 pairs deviated | 
from the mean by more than three times the standard deviation and were rejected. 
Similarly under § 5 (6) 30 of the 166 pairs were rejected. 

The final mean values for n are in very reasonable agreement with the predicted | 
value of 1-5 given by equations (6) and (7) over the wavelength range A=4-2 to | 
8:3m. and the preliminary results on A=1-4m. indicate that the law also holds 
down to this wavelength. 


~~ 


S7. ELEC ERON DENS TE YUN SEE RAIS 


By using equation (7) it is possible to calculate the electron density in the 
meteor trails, since G,P,A are constants of the apparatus, and « and R are 
measured for individual echoes. ‘The measured densities lie in the range 109 | 
to 10” electrons/cm. path. ‘The distribution of these densities in the major | 
showers will be dealt with in Part II. | 

Of particular interest are those cases where visual observations have been | 
made concurrently with the radio echo observations, since it is then possible to i 
measure the electron density of the trails of meteors which are of visual magnitude. | 
Coincidence measurements of this type between a visual meteor and the radio | 
echo on a single wavelength were made during the 1946 Perseid shower (Prentice, | 
Lovell and Banwell 1947), during the 1946 Giacobinid shower (Lovell, Banwell 
and Clegg 1947) and between a visual meteor and the radio echo on two wave- 
lengths during the 1947 Perseidshower. Detailed associations of visual magnitude 
with electron density have not yet been made but the results agree in giving | 
electron densities of 10!° to 10¥ electrons/cm. path for the visually observed 
meteors. ‘Ihe measurements on the 1946 Perseid and 1947 Perseid showers 
indicate that a 5th magnitude meteor which is on the limit of naked eye visibility 
produces approximately 2 x 10'° electrons/cm. path and that a very bright meteor _ 
{~magnitude+1) produces approximately 10” electrons/cm. path. These 
results are in very good agreement with the theoretical calculations of Herlofson 


(1948). 
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§8. COMPARISON WITH OTHER THEORIES 

Previous calculations of the densities of electron concentrations which give 
rise to the transient radio echoes have been made by Eckersley (1940) and Pierce 
(1938). 

Eckersley considered that the echoes were due to clouds of electrons and treated 
the scattering of radio waves as analogous to the scattering of «-particles by heavy 
atoms. On this theory the voltage amplitude returned to the receiver should 
be proportional to A*. Eckersley’s measurements in the wavelength range 
A=15m. to 35m. showed agreement with this \2 law. The measurements 
described in this paper do not cover this wavelength range but it is evident that 
the A? law does not apply over the wavelength range A=1-4 to 8-3 m. discussed 
here. Eckersley’s theory does not explain the critical orientation effect of the 
meteor trails with respect to the aerial beam (Pierce 1938, Hey and Stewart 
1946, 1947, Lovell, Banwell and Clegg 1947) and it also yields values for the electron 
density in visibly observed trails which are in disagreement with Herlofson’s 
calculations (1948). As discussed by Lovell (1948) the transient echo phenomena 
show a sudden change as the wavelength increases beyond 8m. Below this 
wavelength the radio echoes are due to the trails of meteors which are associated 
with the visible showers originating from definite radiants. In the higher 
wavelength range the echoes due to these shower meteors are generally submerged 
in a very high background rate showing a diurnal and seasonal variation (Eckersley 
1940, Appleton and Naismith 1947, Eastwood and Mercer 1948). ‘These 
diurnal and seasonal variations can be explained if the earth is passing through 
a general distribution of meteoric dust, too minute to produce visible meteors. 
Eckersley’s measurements refer to these latter type of transient echoes, and 
although the reasons are not yet understood it is possible that the change from 
the 2 to A*? law may be associated with the change in origin of the transient echoes. 

Pierce (1938) calculated the scattered energy from a meteor trail compared 

with that scattered from an ionized layer, by assuming that the electrons in the 
trail formed a long column with a diameter /arge compared with the wavelength 
and with such a high electron density that the radio wave did not penetrate to 
the inner parts of the trail. ‘These calculations predict the same variation of 
scattered energy with A and R as equation (6) and from this point of view the 
experimental results described here are not in disagreement with Pierce’s theory. 
There are, however, two facts which indicate that Pierce’s fundamental assump- 
tions are in error. a 

(i) Pierce assumes that the meteor trail has already diffused so that its diameter 
is large compared with A before the radio echo is observed. Other aspects 
of the work described in this paper (Part III, see also Ellyett and Davies 
1948) show that this assumption is incorrect, at least for observations in the 
wavelength range A=1-4 to 8-3m. The echo is obtained from the trail in process 
of formation as it crosses the beam of the aerial. The subsequent expansion 
of the electron column by diffusion, to the dimensions assumed by Pierce to be 
the initial state, can then be shown to be a necessary consequence of the change 
of duration of the transient echoes with the azimuth of the aerial beam described 
by Lovell (1948; see also Part III) and explained by Herlofson (1948). ein 

(ii) Pierce’s calculations lead to estimates of the electron density in visible 
meteor trails which differ by orders of magnitude from the densities calculated 
from equation (6) and are thus in disagreement with contemporary theory 
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(Herlofson 1948). Also in order to explain the intensity of a meteor echo 
compared with the intensity of the echo from the F region, Pierce, following 
Maris (1929) has to assume that the entire kinetic energy of the meteor is spent | 
in ionization, whereas Herlofson (1948) shows that only 10-° of the kinetic | 


energy of the meteor is spent in this manner. 
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LETTERS TO °THEsE DIOR 


Polarization of Second Order Raman Effect in Alkali Halides 


We have succeeded in obtaining Raman spectra of sodium chloride, potassium bromide 
and potassium chloride which show polarization effects. The spectra obtained for NaCl | 
correspond to the unpolarized spectrum obtained by Rasetti (1931) and by Krishnan (1946). | 
With a Hilger medium quartz spectrograph using 2537 a. as the exciting line, with a mercury | 
vapour filter to absorb the resonance line after scattering, and a Wollaston prism in the | 
scattered beam, useful spectra were obtained in two days with NaCl and KBr, and in four | 
days with KCl. Single crystals were grown for these experiments. 

In order to achievé this, a wide-angled cone of incident radiation was necessary. In | 
practice, no limitation was put on the beam other than that naturally afforded by the critical 
angle. For high values of depolarization little error results from this, while for low values. | 
for example for KBr, a depolarization of zero would be measured as 0:2. Thus it is clearl > | 
possible to distinguish between “ polarized” and “‘ depolarized ’’ scattered radiation 

The general appearance of the spectra is that of peaks on a continuous background 7 this | 
simply means that under the given conditions one could not distinguish between aril of | 
intensity in a continuum or lines imperfectly resolved. There is a considerable difference 
between the spectra of NaCl and KBr on the one hand and KCl on the other. The two 
former were similar in that there were a number of peaks which were partially polarized 

7. > 
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-and one outstanding peak depolarized with a depolarization in excess of unity. This peak 
was at 232 cm in NaCl and at 127 cm in KBr. 

KCl gave a much weaker spectrum, and there was no outstanding depolarized peak 
all those strong enough to be observed were well polarized. 

It is thought that an explanation can be sought in the approximate equality of mass of the 
two atoms in KCl and the disparity of mass of the constituent atoms in the other two crystals. 
On the cyclic lattice theory of Born, it is sufficient to consider the one-dimensional lattice 
of two atoms to see that peaks of intensity of the continuum present when the ratio of the 
atom masses differs considerably from unity will fade as the ratio approaches unity. This 
at once accounts for the relative absence of peaks of intensity in KCl present in NaCl 
and KBr. 

An attempt is now being made to grow a single crystal of RbBr which has a mass-ratio 
of 1-07 which, if the above explanation is correct, should give polarized spectra similar to 
those of KCl, which has a mass-ratio of 1-1. 

The type of spectrum obtained is seen in the figure, which shows microphotometer curves 
of the vertically and horizontally polarized spectra for rocksalt. In each case 232 cem- is 
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Raman spectrum of rocksalt. 


Note: The numbers 517 B.H. in the diagram are merely the recording paper code and are 
to be ignored. 


marked with an arrow. The dip in the high peak is the position of the exciting (resonance) 
line . The vertical component is privileged in the arrangement used, and this should be 
remembered in comparing the two curves. ‘The position of mercury lines is indicated by 
vertical dashed lines. 


Adam Hilger Ltd., London N.W. 1. A. C. MENZIES. 
30 January 1948. J. SKINNER. , 


KRISHNAN, R. S., 1946, Proc. Roy. Soc. A, 187, 188. 
Rasetti, F., 1931, Leipziger Vortradge, Molekiilstruktur, p. 59. 


REVIEWS OF BOOKS 


Theory and Applications of Electricity and Magnetism, by CHARLES A. CULVER. 
First Edition. Pp. 594, 407 figs. (London and New York: McGraw Hill 
Book.Co.elnc:, 1947.) 25s. net: 

This attractive text book was written to provide an introduction to electromagnetism 


for students at American Universities and ‘“‘ is intended to serve as a basis for advanced 
study in physics and chemistry and also to lay the foundation for courses in electrical 
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engineering ”’, but it should also meet the needs of second and third year students of physics: 
and electrical engineering at British Universities. 

In its thirty chapters the book surveys the elements of electricity and magnetism, 
applied electricity, particle physics and electronics. The theoretical approach to electro-. 
statics, magnetostatics and the magnetic effects of currents conforms to the conventional. 
pattern to be found in most text books, but the treatment is clear and comprehensive, with 
no mathematical demands on the reader beyond a knowledge of elementary algebra and 
calculus. The electrostatic, electromagnetic and practical systems of units are employed 
in the text but a welcome passing reference is made to the metre-kilogram-second system of 
units. 

The distinctive feature of the book which calls for special comment is the very large 
number of practical applications that are described both for their intrinsic importance and. 
co supplement theoretical discussions. To illustrate how the author has justified the word 
‘« applications ”’ in the title of the book, it suffices to state that there are excellent chapters on 
each of the following topics : D.c. measurements (bridges and potentiometers) ; applications 
of the thermal effects of currents (electric furnaces, welding, arcs and lighting) ; electrolysis 
and primary and secondary cells ; thermoelectricity and its uses ; the magnetic effect of | 
currents (ammeters and voltmeters). 

The book contains one chapter on A.c. theory which expounds the vector method and 
another short chapter explaining the j-notation with applications. Alternating current 
instruments claim a chapter and electrical units another. 

The last third of the book is a rapid survey of electronics and atomic physics with a 
stress on practical applications rather than on theoretical details. All this, within the 
author’s terms of reference, is very well done, although the chapter on electricity in gases. 
would seem to indicate that the author has not read sufficiently widely in this subject. 

The final two chapters deal respectively with thermionic tubes and their uses and with 
electromagnetic waves and their applications. 

The student’s interest is stimulated by frequent reference to the history of ideas and 
applications but in this the author proves himself by no means infallible. For instance, |} 
it comes as a surprise (p. 468) to learn that Moseley worked at the Cavendish Laboratory. |f] 
On p. 428 we learn, with reference to the nature of cathode rays, that ‘‘the German school 
of thought held, more or less, to the view that the rays consisted of some sort of wave 
disturbance ; while English physicists felt that they were dealing with some type of electri- 
fied particles’. This statement is inconsistent with the fact that in 1897 both Wiechert and | 
Kaufmann published values of e/m, and that the former was the first to draw the correct | 
conclusion concerning the mass of the cathode corpuscles. On p. 57 it is stated that sun- ' 
spots and terrestrial magnetism have been studied in relation to each other for over 1000 | 
years. ‘This is presumably a misprint for 100 years. 

These, however, are trivialities. More serious is a diagram (figure 338) of a cyclotron 
showing the pole pieces spanning scarcely one-third of the diameter of the dees, and a |} 
diagram illustrating the principle of radar (figure 399) in which six pulses are shown in | 
flight simultaneously between the ground equipment and the airborne target. It is also 
time that writers of text books abandoned the practice of attributing the external magnetic 
field of the earth to an enormous internal bar magnet (figure 24) ; a doublet source near the 
centre is a much closer representation. In this context the author remarks : “‘ The terres- 
trial poles are below the earth’s surface’. Whether they can be reached by a suitable | 
mining operation is not stated. 


Useful sets of test problems accompany most chapters and the book is well bound and 
excellently printed. 

The author ends his book on an optimistic note: ‘‘ The science of electricity and 
magnetism has brought countless blessings to mankind. In the years to come it will make 
the earth a still better place in which to live and have our being”. Be this as it may, the 
author has given practical effect to his belief in writing a book which can give to an imackam 
ative student a vision of what is possible. The book can be confidently recommended. 
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